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The refraction of shock waves at a gaseous interface 


By ROBERT G. JAHN 
Department of Physics, Lehigh University, Bethlehem, Pennsylvania 


(Received 24 April 1956) 


SUMMARY 

A programme of shock tube experiments has been conducted 
to study the refraction of plane shock waves at interfaces between 
two gases. Shocks of strength € = 0-85 (weak) and € = 0-30 
(fairly strong) were allowed to impinge, at various angles of 
incidence, on interfaces between air/CO, and air/CH,, and the 
resulting configurations were photographed through a Mach- 
Zehnder interferometer. From the interferograms, measurements 
were made of the strengths of the reflected waves, and of the angles 
of refraction, and the values were compared with the theoretical 
calculations of Polachek & Seeger (1951). Within the range of 
parameters for which the refraction model assumed by the theory 
is applicable—the so-called ‘regular refraction’ region—the 
observations were in excellent agreement with the theoretical 
predictions. 

When the study was extended to ranges of the parameters for 
which the theory is clearly inadequate, a succession of rather 
complex ‘irregular refraction’ patterns was observed. Although 
these configurations were highly interesting qualitatively, each of 
them involved curved shocks, non-uniform regions of flow, and 
other less simple processes which discouraged any formal 
theoretical analysis. On a less rigorous basis, however, it could 
be shown that these patterns were internally consistent, and that 
each represented a distortion of a regular refraction process which 
was reasonable under the prevailing aerodynamic conditions. 

Certain observations in these refraction experiments appear 
to be of some significance outside the specific problem. (i) The 
sensitivity of strong shock refractions to the values of the specific 
heat ratio y for the two gases suggests a possible technique for the 
measurement of y and its temperature dependence. (ii) ‘Two of 
the irregular refraction patterns display a transition process which 
would be equally appropriate to the onset of the Mach configuration 
in the shock reflection problem. (iii) Some irregular refractions 
can be considered as special cases in the problem of the interaction 
of a shock and a boundary layer. 


I. INTRODUCTION 


Of the many non-linear phenomena available for physical study, 


the 


behaviour of shock waves in gases is one of the neatest and least complex. 


F.M. 2H 








458 Robert G. Fahn 


Even so, shock wave problems are difficult to treat theoretically, and 
computations on even the most fundamental shock interactions tend to 
become involved and somewhat ambiguous. Normally it must fali to 
direct experiment to clarify the theoretical predictions, and to extend the 
understanding into regions too complex for mathematical solution. This 
paper discusses a programme of shock tube experiments which is intended 
to clarify and supplement the existing theoretical information on one of 
the most basic shock interactions—the refraction of plane shock waves at 
an interface between two gases. It is a condensation of a doctoral thesis 
submitted by the author to the Department of Physics, Palmer Physical 
Laboratory, Princeton University. ‘Throughout the paper, a dagger (+) 
is used to indicate specific points discussed in greater detail in the original 
work, which is available on request from the Princeton University Library. 


II. ‘THEORY 

Theoretical description of this problem has been formulated indepen- 
dently by Taub (1951) and Polachek & Seeger (1951). ‘They assume the 
same, somewhat idealized model for the refraction process, which is outlined 
briefly here to introduce the notation and terminology in this paper (see 
figure 1). 

A plane, flat-topped shock J of strength € (€ = pressure ahead/pressure 
behind) travelling in an ideal (inviscid, with constant specific heat and 
zero thermal conductivity) gas (1) is incident at some angle « on a plane 
interface O of another ideal gas (5). Assume that a plane, flat-topped 
shock 7 of strength &” is transmitted at an angle «”, and a similar shock RS 
of strength €’ is reflected at an angle «’. Assume also that the gas interface 
passed over by the intersection X is deflected through an angle A, and that 
each of the five angular regions (1) to (5) formed by this configuration is 
uniform in all its pertinent properties: pressure p, density p, temperature 7, 
velocity of sound a, entropy s, and gas flow velocity u. Such a process 
will be called a ‘regular refraction’. 

The problem is formulated from the point of view of an observer travelling 
with the intersection point XY. Such an observer sees a certain gas inflow 
velocity in regions (1) and (5), u,=u;. The gas flow u, is decelerated 
and deflected through an angle 6 by /, and again decelerated and deflected 
through 6’ by RS. ‘The flow ug, is deflected through 8” and decelerated by 7. 
The problem is solved by appplying the Rankine—Hugoniot shock equations, 
in the form relating change in flow velocity to shock strength, across J, 
RS, and 7, under the boundary requirements that the outflow velocities 
u, and u, be parallel (though not necessarily equal), and the exit pressures 
be the same, that is, under the requirements 6+6’ = 6” = A, and p, = py. 

Polachek & Seeger (1951) increased the range of applicability of their 
refraction model by admitting the possibility that, instead of a shock, the 
reflection could be a Prandtl—Meyer angular rarefaction wave, centred at _X. 
Such a reflection (hereafter denoted by RR) would accelerate the flow u, 
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and deflect it toward, rather than away from the interface. Furthermore, 
they recognized that any particular gas combination was not necessarily 
restricted to just one or the other type of reflection, but that reflected shocks 
might occur for small angles of incidence, and reflected rarefactions for 
more glancing angles, or vice versa. (The angle of incidence at which 
such a change takes place we call «,—the transition angle.) 


I a 








Figure 1. Assumed configuration for the shock refraction process. J = incident 
shock; RS = reflected shock; T = transmitted shock; «, «’, «” = angles of 
incidence, reflection, refraction; X = point of intersection of J, RS, and T 
with interface; OXO original interface; XD = deflected interface; 
A = angle of deflection of interface; u = gas flow velocity; 8, 5’, 5” = angles 
of deflection of flow by shocks J, R, T, respectively. Numeral subscripts 
refer to the 5 angular sectors of the pattern. 


Both Taub (1947) and Polachek & Seeger (1951) have carried out 
extensive numerical calculations, based on their own formulations, of 
the reflected and refracted wave strengths and angles for various gas 
combinations, shock strengths, and angles of incidence. In both cases, 
the automatic calculators can solve the problems only incompletely. For 
some values of the pertinent parameters they return several solutions ; 
for others, none at all. In the former cases it is necessary to choose one 
branch of the solutions on the basis of some boundary value criteria, such 
as the known behaviours in the limits 1 (acoustic case) and «— 0 
(normal incidence). In the latter situations it is presumed that this 

2H2 
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idealized model for the refraction process is inadequate to satisfy simul- 
taneously all of the aerodynamic requirements, and that some other, more 
complex processes must occur physically. 


III, DESIGN OF THE EXPERIMENTS 

The theoretical approach outlined above was given encouragement by 
some early experimental work of Stoner, Woodbridge & Davies (1952), 
Bitondo (1950), and the Princeton laboratory, which indicated that a 
process quite similar to that assumed by the theory did occur physically. 
This process was observed to occur for each gas combination and shock 
strength tested, over a range of angles of incidence from normal up to a 
certain limiting angle, which depended on € and the gas combination. 
In the range of « beyond this limiting angle, other more complex inter- 
actions, called ‘irregular refractions’, were observed, which were clearly 
incompatible with the simple process assumed by the theory. However, 
this early work was largely qualitative, involving shadow photographs and 
streak photographs, from which it was possible to measure wave angles, 
but not shock strengths or pressure fields directly. It was felt that an 
experiment which provided direct measurements of the wave strengths and 
pressure distributions throughout the entire refraction field was necessary 
if the assumptions and solutions of the theory were to be checked thoroughly 
and the complex irregular patterns were to be understood at all. 

Of the two sets of calculations that were available, those of Polachek & 
Seeger (1951), which included the cases both of reflected shocks and 
reflected rarefactions, seemed more appropriate for experimental compari- 
son. (Henceforth we refer to these exclusively.) ‘The results of those 
computations were presented in the form of graphs of 1/€’ vs « and «” vs « 
for seven shock strengths and twenty-four gas combinations. ‘These 
graphs, largely unpublished, were made available from the Naval Ordnance 
Laboratory, by the kind cooperation of Mr L. D. Krider. The experi- 
mental work described below was first directed toward a quantitative 
appraisal of these regular refraction solutions, and then extended to explore, 
qualitatively, the theoretically opaque regions of irregular refraction. 

The Princeton laboratory is fortunate to possess a 5-inch Mach-Zehnder 
interferometer, with which the gas density in the entire interaction field 
can be evaluated. With this equipment available, the experimental problem 
resolved itself into the creation of a suitable shock wave refraction process 
in the test section of the shock tube. This, in turn, depended primarily 
on the accomplishment of a nearly ideal separation of two gases to form a 
plane refracting interface there. To date, the only successful means of 
gas separation we have found for this experiment has been by some’type of 
film or membrane. Such a film must be sufficiently strong to withstand 
the hydrostatic pressure of the heavier gas, and must be impermeable to 
both gases. On the other hand, it must have a sufficiently small mass for 
its acceleration by the incident shock not to interfere appreciably with the 
refraction process. 
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Finding a film that satisfies these requirements is the key to the entire 
experiment. A natural first attempt was to try soap films, which have 
conventionally been used for similar purposes (cf. especially Stoner’s use 
of soap films for refraction experiments). It quickly appeared, however, 
that soap films had several inherent disadvantages in this application. 
(i) They must be drawn in place at the time that they are to be used. (ii) 
They evaporate, thereby contaminating with water vapour the gases they 
are separating. With some gases (for example, CO.) they react chemically. 
(iii) If used in any position other than horizontal, they drain toward the 
lower edge and become non-uniform in thickness. (iv) They stretch 
comparatively easily, and hence tend to bow out under the hydrostatic 
pressure of the heavier gas, thereby establishing a curved interface. (v) 
When broken accidentally, or by the incident shock, they splatter the optical 
flats with soap solution. 

After considerable investigation of this film problem, a more suitable gas 
separation technique was developed using a commercial lacquer film 
(‘Zapon—Aquanite A’, a product of the Atlas Powder Company, Wilming- 
ton, Delaware). Briefly}, the films are cast on a surface of distilled water, 
lifted off on stainless steel frames, and allowed to dry. ‘They weigh, when 
dry, approximately 5 micrograms/cm’, and in thickness are a small fraction 
of one wavelength of visible light. ‘These Zapon films overcome the above 
disadvantages of the soap films, yet are simpler and more convenient to 
prepare. ‘They can be made lighter, in total mass, than the soap films 
that would be needed to hold the same gases over the same size span. It is 
no exaggeration to say that it was the neatness and manageability of these 
thin films which made possible an accurate quantitative refraction experiment. 

The bulk of the apparatus had to be designed around these all-important 
films. For example, it was necessary to reduce the effective width of the 
shock tube test section from 4 in. to 1 in., even though this decreased the 
optical sensitivity of the experiment, since thereby we could gain a factor 
of about sixteen in the strength of a film of given thickness. Also, it was 
necessary to grind narrow slots in the glass windows of the observation 
section, in order that the edges of the stainless steel frames that held the 
films could be recessed therein and not obstruct the gas flow in the region 
of interest. 

The apparatus and technique of the experiment can best be visualized 
by reference to the diagram and photograph shown in figures 2 and 3 (plate 1). 
The principal piece of special equipment is a shock tube test section H, 
made of 1-in. dural plate and having internal dimensions 35 in. x 17} in. x Lin. 
Through the walls of this test section is mounted a cylindrical refraction 
cell K, which contains the two slotted observation windows R, and suitable 
valves and seals for introducing and holding the refracting gas. ‘This cell 
rotates, and with it the interface-forming film J, to permit incidence of the 
shock wave from any angle from normal to glancing. ‘The lower part of 
the test section contains a system of gears and shafts which permits a two- 
directional adjustment of an 18-in. leading plate P. The purpose of this 
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plate, and a similar but shorter back plate QO, is to present the incident shock 
to the gas interface in as neat and uncomplicated a manner as possible. 
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Schematic diagram of shock tube and auxiliary equipment for refraction 
experiments : (a) top view; (b) side view; (c) enlargement of test cell. A, B, 
C, and D are valves controlling the flow of gases to and from the rear half 
of the refraction cell; FE shock tube chamber; F shock tube channel; 


Figure 2. 


G = plastic diaphragm; H = test section; J interface-forming film; 
K refraction cell; L — exhaust tank; .W = mechanical pump; N = low- 
leading plate; Q = back plate; R = slotted glass 


shock front advancing on refraction cell. 


pressure reservoir; P 


windows; S 


The technique actually followed in these experiments might best be 
illustrated in the form of a step by step recipe. 

(i) ‘The films are cast on their frames and allowed to dry. 

(ii) ‘The entire shock tube is completely evacuated, and the appropriate 
pressure is set in the high pressure chamber E. The channel F and test 
section H are then returned to atmospheric pressure. 

(iii) ‘The Zapon film, held on the stainless steel frame on which it was 
cast, is set into place between the slotted windows J. The refraction cell K 
is inserted into the test section at the desired angle and sealed to it. (Note 
that it is necessary to carry out steps (ii) and (iii) in this order. ‘The thin 
Zapon films are so delicate that, if they were present inside the shock tube 
while the high pressure was being set in the chamber, they would be broken 
by the minute pressure pulses emitted by the plastic diaphragm G as it 
stretches into position.) 
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(iv) The leading plate P and back plate QO are adjusted until they nearly 
touch the film interface at its most forward point. 

(v) With valves A and B open, the shock tube channel and test section 
are evacuated again, quite slowly, until p, and p, are equal to the desired 
downstream pressure. (‘The connection through valve A equalizes the 
pressures on either side of the film during the pumping process more readily 
than the longer connection va valve B.) 

(vi) With valve A closed and B open, the refracting gas is now intro- 
duced very slowly into the rear of the cell via needle valve D, replacing the 
air there which exhausts into a large tank L mounted by the side of the 
shock tube. ‘The other end of this exhaust tank is directly connected to 
the shock tube channel, thereby automatically maintaining p, = p;. The 
tank is sufficiently large, and contains bafHes so arranged that practically 
none of the exhaust gas ever drifts into the shock tube. Since the flow 
path of this arrangement is long compared with the size of the various 
orifices involved in the plumbing, p; tends to remain slightly higher than p, 
for any reasonable rate of flow. ‘To compensate for this, needle valve C, 
connecting an evacuated reservoir (V) with the exhaust line, provides a 
further ine adjustment on p;. 

An external indicator of the pressure difference across the film, of the 
extreme delicacy needed to protect the films throughout the above process, 
might be difficult to arrange. But fortunately, since the interior of the 
refraction cell is visible through the windows, the films themselves can serve 
as pressure indicators. Normally, they respond to a pressure difference 
by a slight elastic bowing into the lower pressure region before they break. 
By careful observation of such distortions it is usually possible to control 
the gas flows sufficiently accurately to avoid a rupture. In practice, the 
operator, holding needle valve D in one hand and C in the other, watches 
the oblique reflection of a circular light source from the front face of the 
film. As the film bows slightly in or out under an incipient pressure 
difference, the reflected image becomes elliptical in one direction or the 
other. Keeping this reflected image as nearly circular as possible, the 
operator gradually opens the two valves until a reasonable flow of gas is 
established. 

(vii) After the gas has been flowing thus through the rear half of the 
cell for five or ten minutes, the shock tube is fired, sending a plane, flat- 
topped shock wave down the tube toward the test section. ‘I'wo sharp 
plates extending up the tube from the test section literally slice out of this 
shock front a 1-in. strip, which then travels over the leading plate and is 
presented to the 5 in. x 1 in. gas interface in the refraction cell. 

Spark photographs of parallel fringe interferograms of the resulting 
refraction process are taken through the Mach-Zehnder interferometer. 
These interferograms are analysed by the standard superposition method 
(see Bleakney & Griffith 1954) to determine the density field, from which 
the shock strengths and pressure distributions can be computed. 
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Measurements of the wave angles can also be taken directly from the 
interferograms. 

White light interferograms are obtained simultaneously, which aid in 
determining the integral fringe shift across the various density discon- 
tinuities. ‘I'wo light screens record the velocity of the incident shock wave 
as it travels down the tube and thereby offer a double check on its strength. 

There are three significant restrictions on the ideality of this experiment: 
the effect of the finite mass of the film, which will be discussed later; the 
inability to obtain completely pure gas in the refracting cell, which must be 
considered in the evaluation of the dataft; and one other, which arises as 
an unavoidable consequence of the physical boundaries at which the 
observed refraction process begins. Nowhere in the theoretical treatment 
is there considered the practical limiting condition that the outflow from 
the reflected and transmitted waves must be supersonic in their respective 
gases with respect to the intersection point if regions (3) and (4) are to 
retain portions undisturbed by foreign signals. ‘That is, regardless of the 
apparatus, the ‘corner’ at which the refraction process starts will radiate 
signals which travel out into the flow with the local sound velocity. For 
sufficiently glancing angles, such corner signals can ‘catch’ the intersection 
point, distort the reflected or transmitted waves, or both, and destroy the 
uniformity of all of regions (3) and (4), thereby violating the theoretical 
assumptions. ‘This effect imposes on the regular refraction process two 
limiting anglest, x,, and «,,, beyond which the observations cannot be 
expected to agree with the theoretical predictions. Fortunately these 
angles usually occur near the end of the theoretical regular refraction region, 
and therefore this situation does not require serious consideration until 
one wishes to study the irregular refraction configurations (see § V1). 

Five dependent quantities are available for measurement from the 
interferograms obtained by the technique outlined above: €’, &’, a’, a”, 
and A. Of these, &’ 2 dx” are the most significant to study. €’, the 
strength of the reflected wave, is the most sensitive index to the refraction 
process, changing by large fractions of itself for comparatively small changes 
in the experimental conditions ; and although the reflected waves are always 
quite weak, and the fringe shifts across them correspondingly small, the 
interferometric technique permits sufficiently accurate measurement of €’ to 
make it an ideal quantity to study. The angle of refraction ~”, though less 
sensitive to the parameters, is always sharply defined on the interferograms. 

The remaining three quantities cannot be obtained with comparable 
accuracy. ‘The evaluation of the refracted shock strength €” is rather 
sensitive to the optical properties, and hence to the purity of the refracting 
gas, which varies appreciably from shot to shot. The reflected wave is 
so weak that the angle «’ it makes with the interface is not measurable with 
adequate precision. (This measurement is particularly difficult when the 
reflection is a rarefaction fan.) Finally, the measurement of the deflection 
angle A is inaccurate because the deflected intertace is rather broad and 
indistinctly defined on the interferograms, and because this region of the 
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pattern is complicated by spurious signals from the corner. Fortunately, 
it is also true that the angles «’ and A are less sensitive indices to the refrac- 
tion pattern than the other variablest. 

Therefore, in these experiments we concentrated on the measurement 
of €’ and «”. Since one of these pertains to the reflection and one to the 
refraction, accurate evaluation of these two can check the entire regular 
refraction solution. In the analysis presented in the following sections, our 
approach is to superimpose experimental data of €’ and « corrected for 
conditional variations, on the theoretical curves of 1/€’ vs « for the appro- 
priate shock strength and gases. We also compare theoretical curves of 
a” vs x with the experimental values and show that these are less sensitive 
to the conditions of the problem. In a few cases, we include a few spot 
checks on «’, €”, and A to preclude any major deviation from the assumed 
configuration from passing unnoticed. 

Once the feasibility of the technique outlined above was established by 
a few preliminary tests, it became necessary to select the range of initial 
conditions that a systematic set of experiments should cover. In this 
experiment, as in almost any practical refraction situation, the five inde- 
pendent variables of the theory (€, «, y,, ys, @/a;) are embodied in three 
physical conditions, €, «, and the gas combination. Of the many combina- 
tions of these conditions which could be investigated in the shock tube, 
it was desirable to select those cases that could be handled most simply, 
yet still provide a truly representative check of the idealized theory. In the 
experiments described below, study was made of the refraction of shocks of 
strength € = 0-85 and & = 0-30 incident at angles from « = 0 to glancing, 
on interfaces between air and carbon dioxide, and air and methane (hereafter 
denoted by air/CO, and air/CH,). The reasoning behind these choices 
can only be outlined briefly heret. 

The theory predicts two distinct classes of solutions, those which 
involve a reflected shock at « = 0, and those which involve a reflected 
rarefaction at « = 0, depending on the condition 


rs = Svsl(vs+1)+7s— VE] | 2? 1 > 1 





P, luluitD+(i-De 4, 

More detailed study of the way in which the gas properties enter into 
the problem, and of the results of the computations for many varied gas 
combinations, indicates that, beyond this broad distinction, there is no 
further significant subdivision of the solution types with respect to gas 
combinations. 

In other words, we can reasonably expect that one representative gas 
combination from each of the two major classes will adequately simulate 
the effect of any possible arrangement of the conditions &, «, a@,/a;. 

An obvious simplification in the experiment is gained if we choose the 
gas (1), which in our arrangement must fill the entire shock tube channel, 
to be air. For the refracting gases (5), other than the requirement that one 
be ‘faster’ and one ‘slower’ than air, the only obvious properties to consider 
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are their manageability, index of refraction, and the ratio of their sound 
speed to that of air. In preliminary experiments, refracting gases which 
were considerably different from air, such as helium and freon, were tried 
on the theory that the magnitude of the refraction effects would be largest 
for the most dissimilar combinations. However, the large density differ- 
ences involved required rather heavy films to withstand the hydrostatic 
pressures. Also, the large differences in index of refraction made it 
impossible to obtain good optical fringes over the entire field of both gases. 
In the long run, it was found best to study two less severe gas combinations, 
air/CH, and air/CO,. Although the refraction effects are somewhat less 
pronounced for these cases, the increased precision with which the experi- 
ment can be performed more than compensates for this, and a more accurate 
check on the theory can be made. 

The shock strengths € = 0-85 and 0-30 were chosen, under considerations 
of operating convenience and pressure limitations of the shock tube, to 
provide the most reasonable check on the é-dependence of the refraction 
theoryt. 

The specific details of these experiments and the results obtained from 
them are described in the remaining three sections of the paper. § IV 
covers the weak shock experiments, and their results in the regular refraction 
regions. § V presents the strong shock, regular refraction experiments. 
In § VI we conclude with a description and brief analysis of the many 
curious patterns observed in the irregular refraction regions at both shock 
strengths and for both gas combinations. 


IV. EXPERIMENT 1: THE REGULAR REFRACTION OF WEAK SHOCKS 

The first experiment studied the refraction of shocks of strength € = 0-85 
at interfaces between air/CO, and air/CH, for the dual purpose of testing 
the technique and apparatus outlined above, and of checking the regular 
refraction solutions at this shock strength. ‘To accomplish this, refraction 
interferograms, samples of which are shown in figure 4 (plate 2), were taken 
in the manner described in $III. From these, measurements were made of 
reflected wave strengths, angles of refraction, and pressure distributions 
throughout the interaction field in sufficient quantitative detail to permit 
reasonable comparison with the appropriate theoretical computations. 

Unfortunately, the numerical solutions computed by Polachek & 
Seeger (1951) involve approximate values for the pertinent gas constants. 
Specifically, they have assumed that the ratio of specific heats, y, for both 
CH, and CO, is exactly 4, and have rounded off the air/CH, and air/CO, 
sound speed ratios to convenient decimals. Before attempting any experi- 
mental comparison with their results, therefore, it was necessary to adjust 
their solutions by using more realistic values for the gas properties. Since 
these are rather small adjustments, it is adequate to evaluate them by 
approximate methods. For example, we chose to find the corrections at 
a few convenient special angles, such as « = 0, x,, «,, where the computa- 
tions are simpler, and then to interpolate the remainder of the solutions 


graphicallyt. 
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The conditions of the experiments could not be held sufficiently close 
even to these more realistic theoretical assumptions to permit immediate 
comparison of the results. Each observation had first to be adjusted for 
the experimental deviations from ideality mentioned in § III. 

(i) The theoretical curve is based on € = 0-85. The shock strengths 
actually used differ by as much as 0-025 from this value. The normalization 
was handled by cross-plotting graphs of 1/€’ vs € for various values of « 
from the theoretical graphs of 1/€’ vs « for various vajues of €. Taking 
the slope of the former at the point € = 0-85 yields a function of « which, 
when multiplied by the deviation of the experimental value of € from 0-85, 
is just the correction desired. 

(ii) Any incompleteness in the gas replacement process resulted in an 
impurity of the refracting gas and a consequent discrepancy from the 
values of y; and a,/a; assumed by the theory. ‘The degree of this impurity 
was measured from a ‘ blank’ interferogram, taken just before the shock tube 
was fired. ‘The fringe difference across the boundary could be computed 
from the known index of refraction for the pure gas, and any discrepancy 
between this difference and that measured on the ‘blank’ could be used to 
determine the relative impurity. With this measure of the efficiency of 
the replacement process, and knowing the composition of the bottled gas, 
it was then possible to calculate a new y, and a; for the impure refracting 
gas; and from these, corrections to the experimental observations could 
be madef. 

(iii) The effect of the finite mass of the interface-forming film on the 
refraction configuration had to be checked experimentally. For this 
particular case, study was made of the density fields resulting from the 
normal incidence of shocks on films of various thicknesses, with air on both 
sides. The effect of films even five times more massive than those regularly 
separating the gases was found to be surprisingly small. Although sharp 
pulses were reflected from the films, their effective strength, that is, the 
ratio of the pressure before and after the pulse, was never larger than the 
0-001 uncertainty with which the pressure ratio could be measured. 

It is conceivable that some other complications could arise from this 
cause for gas combinations other than air/air and at angles other than « = 0. 
Isolated tests were made at other such conditions, the net result of which 
was to support the assumption that in neglecting the mass of the film the 
error introduced was not greater than the uncertainty in the interferometric 
measurement of 1/é’f. Actually, in the preparation of the data, account 
was taken of the almost trivial correction suggested by extrapolating the 
film test results to zero thickness. 

In considering the results of this experiment, it is first to be emphasized 
that the qualitative nature of the interferograms taken in the regular refrac- 
tion regions agrees closely with the assumptions of the theory. Namely, 
within the limits of the observation, each pattern is composed of an incident, 
reflected, and refracted wave, which are radially straight for a finite distance 
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from the point of intersection with the interface, and between which are 
regions of uniform pressure. In no case has any evidence been seen of 
other signals, either continuous or sharp (except, of course, the corner 
signals mentioned previously) or of any non-uniformity in the angular 
regions. 

The quantitative results may best be seen in figures 5 and 6. In figure 5 
are plotted the observed values of 1/€’ and «. The dotted curves of 1/&’ vs « 
represent the original solutions of Polachek & Seeger (1951). The solid 
curves are those same solutions, corrected by the use of more realistic gas 
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Figure 5. Experimental values of reflected wave strength &' superimposed on 
theoretical predictions, for air/CO, and air/CH, with € = 0°85. Dotted 
curves represent original Polachek & Seeger (1951) solutions. Solid curves 


represent the same solutions, adjusted for realistic gases. 


constants. ‘The experimental points themselves have been individually 
corrected for variation in incident shock strength, gas purity, and film mass, 
wherever appropriate. ‘The error boxes reflect the ability to evaluate 
1/é’ from the interferograms with an accuracy of from +0-0005 to +0-0020, 
dependent on the particular interferogram. 

Also available for comparison with the theory are the experimental 
values of «" and «. In figure 6 the dotted lines of «” vs « represent the 
original calculated solutions at € = 0-85; the solid lines are those solutions 
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corrected as above. The experimental points that are superimposed on 
them (error boxes are +0-5°) have again been corrected for gas impurity, 
but not for their variations in incident shock strength, since small changes in 
€ caused no observable differences in «”. As evidence of this insensitivity, 
these graphs also show (dashed lines) the computed solutions for a much 
stronger shock, € = 0-30. The curve for € = 1-00, when plotted, could 
barely be distinguished from that for € = 0-85. 
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air/CO, and air/CH, with € = 0-85. For explanation of dotted and solid 
curves, see figure 5. Dashed curves represent theoretical solutions for 
f'= 0-30 
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A few spot checks were made on the other three measurable quantities, 
é”,a’,and A. ‘These quantities were not studied so thoroughly as é’ and «’, 
and in no case was any attempt made to correct the data for variations in the 
experimental conditions, or to adjust the computations for more realistic 
gas combinations. Rather, these measurements, a few samples of which 
are shown in tables 1 and 2, served merely as a reassurance that no major 
deviations from the theory were being overlooked by restricting the study 
of these quantities. 
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V. EXPERIMENT 2: THE REGULAR REFRACTION OF STRONG SHOCKS 
The interferograms obtained of the stronger shock refraction patterns, 
a few samples of which appear in figure 7 (plate 3), appear much the same as 
those taken of the weaker shock case. 


quantitative changes in the various parameters, notably the increase in 
strength of each of the three signals, and the change in relative sizes of the 
angles of the configuration. 


Table 1. Rough spot checks of «’ and A. 
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€” theo. ( £&') | €” exp 
Air/CO, 
0-854 | 0-851 
0-858 | 0-859 
Sir/CH, 
| 
0-876 | 0-878 
0-850 | 0-851 








Table 2. 


Rough spot checks of &”. 


The deflection of the interface, A, for example, 


There have been, of course, certain 
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is uniformly larger, in accordance with the higher flow velocity set up by 
the stronger incident shock. Similarly, the compounding of this higher 
flow velocity with the relative velocity of the reflected wave has tended to 
increase «’ for those cases where a+’ > 90° and to decrease «’ for those 
cases where a+’ < 90° (i.e. the flow ‘blows back’ the reflected wave). 
Also, the angle of refraction x” now assumes a value closer to « for both the 
‘slow/fast’ and ‘ fast/slow’ cases, for all « < «,f. 

An interesting feature is exhibited in the white light interferograms, 
figures 7c and 7d (plate 3). The reflected shock from the air/CO, interface 
at « = 0 (figure 7c) is quite sharp and flat-topped. However, the reflected 
rarefaction from the air/CH, interface (figure 7d) has already spread out 
considerably. Its trailing edge is actually behind the original interface 
position. Also note that the leading edge of this rarefaction is preceded by 
a slight compression bump, due presumably to the finite mass of the film. 

All of the interferograms obtained in the regular region at this shock 
strength exhibit excellent radial straightness of the signals, and uniformity 
of density and pressure within each of the five angular segments, thereby 
displaying again a close correspondence between the experiment and the 
theoretical problem. 





Air/CO, : é 1/€’ 


39-0+ 1-0 0-305 1:0769 


Measured value 





Correction for 9°, impurity 








of CO, 0-4 = --0-0078 
Correction to make results 

relate to 0-300 — 0-005 + 0-0016 
Corrected value 39-4 0-300 1-0863 

















Table 3. Typical sample of corrections applied to measured values of 1/€’ and «”. 


Measurements of 1/é’ and «” were extracted from the interferograms 
just as before. Again, certain corrections needed to be made before the 
observed values could be compared with the theoretical calculations. The 
adjustment of the Polachek & Seeger (1951) solutions for the realistic gas 
constants was repeated as above. The corrections of the experimental values 
for the mass effect, the purity of the refracting gas and variation in incident 
shock strength also proceeded as described for the weak shock caset. As 
an illustration of the relative importance of the various corrections for a 
strong shock experiment, table 3 shows a breakdown of the computations 
for one typical pattern. 

In evaluating the results of these experiments it is well to bear in mind 
that the optical sensitivity of the measurements was considerably reduced 
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by the necessary use of a low p,. Assuming that it is possible to measure 
an interferogram to +0-03 fringes, 1/€ could be determined to +0-015 
and 1/é’ to + 0-006. Considering the other uncertainties of the experiment, 
the error limits for 1/€’ range from +0-006 to +0-009. The angle «’” 
could normally be determined to +1-0°, except in a few cases where the 
corrections were sufficiently large or uncertain to warrant +2-0° limits. 
Despite this decrease in precision, the results, plotted in figures 8, 9, 
and 10, show that the experiment is still sufficiently sensitive to trace out a 
curve of agreement with the corrected Polachek & Seeger computations. 
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Figure 8. Experimental values of 1/&° superimposed on theoretical predictions, for 


air/CO, with £ 0-30. For explanation of dotted and solid curves, see 


figure 5. 


In particular it is interesting to note that the experiment is again capable 
of choosing between the theoretical curves for the realistic and idealized 
gases, even though the values of y and the sound speeds were changed only 
by 2% or 3° Consideration of the sensitivity of the experiment to changes 


in y suggests the possibility of designing a more appropriate version of the 
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Figure 3. Refraction test section, with rotatable cell removed. 
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Figure +. Refraction interferograms for shock strength € = 0°85: (a) air/COs, | 
x 25 (reflected shock); (b) air/CO,, a 63 (reflected rarefaction) ; 
(c)- ait/CH,;. 27 (reflected rarefaction); (d) air/CH,, « = 45° (reflected 
shock). J, R, T = incident, reflected, and transmitted waves; O original 
interface; arrows indicate directions in which J, R, and T are travelling. 
Fringes shift toward lower right for increase in pressure. 
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Figure 7. Refraction interferograms for € = 0-30: (a) air/CO,., « = 53° (reflected 
shock); (b) air/CO., « 63 (reflected rarefaction); (c) air/COs, 
x 0° (reflected shock); (d) air/CH,, « 0° (reflected rarefaction). 
I, R, T incident, reflected, and transmitted waves; O = original interface ; 
D deflected interface; arrows indicate directions in which waves are 
travelling. Fringes shift down for increase in pressure. 
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Figure 14 (a). Air/CHy,, regular refraction. (J, RS, T = incident, reflected, trans- 
mitted shocks; O, D original, deflected interfaces; C corner signal.) 
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Figure 14 (b). Air'CH,, irregular refraction,a >«,. (I, RS, T = incident, reflected, 
transmitted shocks ; O, D = original, deflected interfaces; r = sonic rare- 


faction front.) 
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Figure 14 (c). Air 'CHy,, irregular refraction,a >a 4. (J, RS, T = incident, reflected, 
transmitted shocks; O, D = original, defiected interfaces; S = stem shock; 
SS = slipstream; B = back plate.) 
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Figure 14 (d). Air/CH,, irregular refraction, « x;, (U, RS, T = incident, reflected, 
transmitted shocks; O, D original, deflected interfaces; c, r sonic 
compression, rarefaction fronts.) 
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Figure 14(e). Air CHy, irregular refraction, a Me ie, eS, incident, 
reflected, transmitted shocks; O, D original, deflected interfaces; 


B back plate; c, r — sonic compression, rarefaction fronts.) 
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Figure 14(f). Air CO,, regular refraction (J, 7 — incident, transmitted shocks; 
RR reflected rarefaction; Qo DBD original, deflected interfaces; 


B — back plate; TR = reflection of T from B.) 
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Figure 14 (g). Air COsg, irregular refraction, « x, (CU, T = incident, transmitted 


shocks: RR reflected rarefaction; O, D original, deflected interfaces.) 
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Figure 14(h). Air CO,, irregular refraction, « x,. (1, T — incident, trans- 
mitted shocks: RR reflected rarefaction; O, D original, deflected 
interfaces; B back plate; é £ compression, rarefaction fronts; 


TR reflection of T from B.) 
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Figure 18. Irregular refraction interferogram for air CH, with & 0-30, a X2, 
1 


and back plate B set parallel to the original interface O—a simulation of a 


thermal boundary laver. 
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normal refraction case, specifically to measure y for either the refracting 
gas or air. For example, in going to progressively stronger incident shocks, 
it would be possible to study the dependence of y on temperature. (These 
ideas were first suggested to the author by Professor W. Bleakney.) 

The agreement between the theoretical and experimental values again 
permits endorsement of the Polachek & Seeger computations, and their 
choice of the physically realistic branch of the numerical solutions, for the 
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Figure 9. Experimental values of 1/€’ superimposed on theoretical predictions, for 
air'CH, with € —- 0-30. For explanation of dotted and solid curves, see 


figure 5 


specific cases studied here. More generally, by the reasoning offered in 
$111, it can now be inferred that the theory of regular refraction from which 
these computations were made is valid for any shock strength from the 
acoustic limit, € = 1-00, to € = 0-30, and for any gas combination having 
constant values of y,, y;, and a,/a; in this range of values of €. 
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Figure 10. Experimental values of x” superimposed on theoretical predictions, for 
air CQO, and air CH, with & 0-30. Solid curves are adjusted Polachek & 
Seeger (1951) solutions. Dotted lines follow experimental points into 
irregular refraction regions. 
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VI. EXPERIMENT 3: IRREGULAR REFRACTION 

In the mathematical formulation of the theory of shock wave refraction 
there arise certain upper limits on the angle of incidence x, beyond which 
no real solutions can be obtained for the assumed configuration. The 
region from x = 0 to the smallest of these limiting angles for a particular 
retraction problem is defined as the region of regular refraction for that 
problem. So far our experiments have been concerned exclusively with 
this region and have indicated the validity of the refraction theory within it. 
We now consider the refraction problem for angles of incidence beyond the 
limiting angles, where the idealized theory becomes inadequate. 

The failure of the theory to produce realistic solutions beyond the first 
limiting angle implies physically that the assumed regular refraction process 
has been replaced by some other type of interaction. Actually, we find 
experimentally that several different configurations occur in this region, 
each of them more complex than the regular refractions. ‘These irregular 
interactions are a significant part of the overall refraction problem. In no 
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case do the regular refraction solutions persist through to « = 90 ; in fact, 
for many reasonable gas combinations, the first limiting angle is quite 
small, so that the irregular regions for these problems are larger than the 
regular ones. For this reason alone, some survey of the irregular phase is 
needed to complete the experimental refraction study. But even beyond 
their significance in the specific problem, irregular refractions are interesting 
phenomena in their own right. 

In presenting the results of this group of experiments, the approach 
will be simply to display examples of the various types of configurations 
that were observed in the irregular refraction regions, and describe the 
principal features of each pattern. ‘Then, in the absence of any quantitative 
theory with which to compare the observations, it is interesting to consider, 
at least qualitatively, the correspondence between the theoretical limiting 
conditions of regular refraction theory and the irregular patterns that result 
when those conditions are exceeded. For example, it is possible to show 
that the various patterns are related to each other in a continuous sequence, 
and that they may be considered to evolve from one to another, beginning 
with the regular pattern, under the influence of those aerodynamic require- 
ments which bring about the limiting conditions. In order to appreciate 
this relation of the irregular patterns to the limiting conditions, it is best 
first to discuss the theoretical critical angles and the physical restrictions 
they represent. 

Recalling the formulation of the regular refraction problem outlined 
in §II, it is possible to anticipate several situations which will make one or 
more of the assumptions of the idealized model untenable. ‘hese may 
be discussed most conveniently in terms of the smallest angles of incidenc« 
at which, for fixed values of the other parameters, the complicating effects 


first occur. ‘These limiting angles are as follows. 


(i) a, and x, 

The theory assumes that the reflected wave will be either a radially 
straight shock, or a Prandtl-Meyer expansion fan, both centred on the 
point of intersection of the incident shock with the gas interface. For 
sufficiently large z, however, region (2) becomes subsonic with respect to 
the point Y. Hence, the angle at which this situation begins is a limit on 
the assumed configuration. Note that this angle, denoted by ~,, depend: 
only on the value of y for gas (1) andon the incident shock strength. Actually, 
if the reflection is a shock, the limit x, is not strictly appropriate. Rather 
than require the outflow from the incident shock to be just sonic with 
respect to the reflection, we should require the flow to be sufficiently super- 
sonic to permit the existence of the finite shock demanded by the pressure 


and deflection conditions. For example, for € = 0-30 in air, « 6i:S° : 
however, for air/CH,, a reflected shock with 1/€’ = 1:20 occurs at this 
angle. Sucha shock requires the inflow Mach number to be approximately 


4 


1-09, and hence the true limit, which we denote by « ,. is some 4 or 5 degr 


ie@L 


lower than ~,. 
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(2) «, and x; 

In those cases where a, < a; (‘slow/fast’), x” normally exceeds «, and 
thus reaches 90° at some « -~ 90°. Expressed more rigorously, it is a 
consequence of the assumed pattern that Ma,/sina = M’"a,/sin«”, where 
M”, the Mach number of the transmitted shock, is fixed, independently 
of this relation, by the pressure and deflection requirements of the problem. 
Clearly, if the value of WM” computed from these requirements exceeds 
Ma,/a;sinz, the value of «” predicted for the assumed pattern becomes 
imaginary. Physically this must correspond to a limit on the regular 
refraction configuration. ‘The angle of incidence at which this effect 
begins is denoted by «,, and is given by 

sin?a. = (re) = (3)(22) Lina + 1) + (ya ~ DE) ; 
M'a;, Qs) \yi/[(ys+ IVE) + (y5— DEI 
(in calculating this angle for tabulation, Polachek & Seeger (1951) have 
introduced the approximation 1/f’= 1. In some cases, particularly as 
€-» 0, this leads to sizeable error, and even to a misplacing of this limit 
with respect to the others.) 

In one sense «, corresponds to the angle of total internai reflection in 
optics, since, at this angle, the transmitted signal is travelling parallel to the 
interface. ‘lhe analogy is incomplete, however, since in no case short of 
an infinite impedance ratio across the interface will the incident shock be 
totally reflected. 

It is expedient to define another critical angle x; which is related to, but 





always larger than, ~, by 
sina; = ——. 


When » = x;, the point of intersection of the incident shock, X, travels 
along the interface with a velocity just equal to a;. For « > a, therefore, 
region (5) becomes subsonic with respect to this intersection point, and 
consequently the transmitted wave can no longer be a shock. Also, some 
of the region ahead of X, both in (5) and (1), can now be influenced by things 
happening at and behind LX. 

(111) a, 

There exists one other limit, which is somewhat more intrinsic than 
the others, in that it concerns the extent of compatibility of the pressure 
and deflection requirements. Stated simply, as « becomes larger, an angle 
of incidence is reached beyond which it is no longer possible to satisfy 
both the pressure and deflection requirements by any arrangement of 
strengths and positions of the reflected and transmitted waves. The 
situation is sufficiently analogous to the extreme angle occurring in the 
theory of regular reflection to justify the same notation, ~,, for this regular 
refraction limit. Beyond this angle, the equations produce no real solutions. 
In most cases, the theoretical curves of 1/£’ vs x double back to presumably 


unrealistic alternative values of €’ at ~,. 
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‘The significance of x, is increased by the fact that it invariably precedes 
x, for any gaseous refraction. ‘The proof of this can be stated very simply. 
At x = a, x” = 90°, hence 8”, the deflection of the flow by the lower path, 
is zero. ‘Then, in order to have +8’ = 5”, 8’ must be equal to —4, and 
this is obviously impossible with €’ an order of magnitude weaker than &. 
Therefore the x,-condition must have occurred before «=a, In other 
words, a regular refraction pattern with x” = 90° cannot be achieved. 

‘Thus it must be x, or x, (or « ,, whichever is appropriate) that establishes 
the extent of the regular refraction regime and introduces the first irregular 
pattern. ‘The order in which these limits occur depends on both the gas 
combination and the incident shock strength. For most combinations in 
which a; is substantially larger than a,, the deflection condition becomes 
critical before region (2) becomes subsonic. In these cases, x, occurs 
< a,, however, the 


first. kor most combinations in which a; = 4a,, or a; 
limiting angles x, and «, (or x,) are quite close and the order of their occur- 
rence depends on the incident shock strengthf. 

(At this point it might be argued that, for any given refraction problem, 
all limiting angles except the first have only academic significance, for they 
are derived as limits to a regular refraction process which terminates before 
they are reached. While it is true that the numerical values of these angles 
computed on this basis are not entirely realistic, the requirements on the 
refraction process that they represent are still pertinent, even for the irregular 
configurations, and the effects resulting from violating these requirements 
are observed experimentally at sufficiently large x. In fact, the effects 
associated with exceeding x, and %;, which can never be the first limiting 
angles, are the most pronounced of all the irregular interactions.) 


(tv) x,, and x,, 

In any physicaily realistic refraction situation, it is also necessary to 
consider two other practical limits to the existence of a regular pattern. 
As mentioned in § III, the ideal regular configuration can be realized 
experimentally only so long as the outflow in regions (3) and (4) is fast 
enough to prevent signals from the ‘corner’ from overtaking the inter- 
section point. Beyond some limiting angle of incidence, this condition 
cannot be met and various signals, associated with the physical corner at 
which the refraction process began, are able to overtake and contaminate 
the region of interest in the pattern. ‘his complication is not, of course, 
peculiar to the refraction problem. Almost any two- or three-dimensional 
shock interaction experiment must reckon with this situation in one form 
or another (in particular the Mach reflection problem, cf. Fletcher, ‘Taub 
& Bleakney 1951). ‘lhe persistent occurrence of corner interference in 
such problems has often inspired philosophical discussion about the futility 
of considering them without complete specification of the physical corners 
at which the processes will begin. Rather than digress into that arena, 
however, we will merely state here the nature of the corner signals in the 
refraction problem, and indicate methods of treating this complicationf. 
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[In a refraction problem, the corner signal may be either a compression 
or a rarefaction front. Its phase and strength are determined by (1) the 
phase and strength of the reflected wave, (i1) the interaction of the transmitted 
shock with the back boundary, and (iti) the nature of the mechanical corner 
formed by the leading plate and back plate, around which the flow behind 
the transmitted shock must deflect (see figure 11). In a sense, each of 





| bi. Y in a retraction experiment. ‘Lhe rerlection 7'R of the 
t lt leflection of the flow f around the corner formed by 

t nd t ick plate B, and the finite length of the reflected 

Omposite corner signal ( which radiate s tron 

peed of sound in the two gases. In general, these 

i secondary signal c’ is transmitted across the deflected 

the { to t slower medium, ahead of the original cornet 


these three effects contributes a component to a composite corner signal 
which radiates from the point of origin ot the process with the local velocity 
of sound, compounded vectorially with the flow velocity. 

it is the ‘catching up’ of this composite signal with the intersection 
point in the incident and refracting media that constitutes the experimental 
limits x,, and lf this corner signal is a net rarefaction, the effect of its 
overtaking a shock will be to weaken the shock and curve it backward (making 
it less normal to the flow) ; if itis a net compression, it will tend to strengthen 
the shock, thereby curving it forward (more normal to the flow). 


Note that the strength and phase of the first component is fixed auto- 


matically for a given refraction situation, but that the latter two are under 
the control of the experimenter who, by changing the angle of the back 
plate, can vary the strength of these components and hence the strength of 
the composite signal. In practice, use was made of this valuable feature via 


two complementary procedures which served to separate the effects associated 


entirely with the idealized irregular refraction process from those which are 


dep ndent on the cornet situation. 
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(i) A given irregular refraction was performed with the back plate set 
to form a corner signal which was a net compression. ‘The same refraction 
was performed with the back plate reset to form a net rarefaction corner 
signal. ‘Those effects which appeared the same for both cases were assumed 
to result solely from the refraction process. 

(ii) ‘The back plate was adjusted so that the various contributions to 
the corner signal compounded to zero in the region of interest. This could 
be done by setting the back plate to some intermediate angle such that the 
net contribution from the mechanical corner and the transmitted wave 
interactions just balanced the uncontrollable component associated with the 
reflection process. In this way, the most interesting part of the refraction 
pattern could be observed, free from corner interference. 
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Figure 12. The critical angles for air/CHy, vs incident shock strength. Circled letters 
refer to corresponding interferograms in figure 14 (plates 4 to 8). 

The considerations outlined above suggest the following physical 
behaviour. As the angle at which the shock impinges on a specific gas/gas 
interface is made more and more glancing, an angle of incidence is reached— 
the first limiting angle for the case involved—beyond which the regular 
refraction pattern cannot occur, but is replaced by a new configuration. 
We expect that this new configuration will be related to the regular pattern 
it has replaced, in fact that it will be a distortion of it, and that this distortion 
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will be such that the flow requirements that brought about the demise of 
the regular pattern are now reconciled. As the angle of incidence is increased 
further, a limiting condition for this irregular pattern may be reached, at 
which point a new distortion of the configuration must occur to accommodate 
the second critical requirement. The pattern will continue to evolve in 
this way, adjusting itself to each new requirement by some distortion, 
until « reaches 90°. 

In figure 14 (plates 4 to 11) are displayed samples of the super-critical 
patterns that were observed for the air/CO, and air/CH, gas combinations. 
Figures 14a to 14e show an essentially complete series of irregular configura- 
tions for air/CH,, arranged in order of increasing angle of incidence. 
Figures 14f to 14h show a similar series for air/CO,. Each interferogram 
is accompanied by a line diagram to clarify the position and nature of the 
various signals involved in the interaction. 
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Figure 13. The critical angles for air/CO, vs incident shock strength. Circled 
letters refer to corresponding interferograms in figure 14 (plates 9 to 11). 


The relation between a given configuration and the various critical 
angles may be visualized by reference to figures 12 and 13, in which are 
plotted the pertinent critical angles as functions of the shock strength for 
the two sets of gases. On these figures, the regions characterized by the 
subsequent interferograms are indicated by circled letters corresponding to 
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the appropriate patterns in figure 1+. Note especially that the critical angles 
have substantially different values at the two shock strengths used, namely 
€ = 0-85 and € = 0-30. This feature, and the previously mentioned extra 
degree of freedom available in the back boundary, were utilized extensively 
to separate the various irregular refraction effects from one another and from 
the corner effectst. In a few cases, two critical angles fall so closely together 
that it is not possible to separate out their two effects by a single picture 
taken between them, and less direct methods have had to be used to study 
that region. Nevertheless, it has been possible to construct, from inter- 
ferograms such as these, a rather complete catalogue of the irregular 
refraction effects and their relation to specific critical conditions. ‘This 
may best be illustrated by brief separate discussions of each of the 


patterns. 


The air,CH, irreguiar refraction patterns 

Figure 14a (plate +) shows the regular refraction pattern that exists for 
the air/CH, combination just below the first critical angle. Note that all 
three shocks are straight, and that the angular regions between them are 
uniform in the neighbourhood of the intersection point. In the regions 
behind the reflected and transmitted waves, it is possible to observe a weak 
rarefaction front. ‘This is the corner signal for this particular arrangement 
of the boundaries; at this angle of incidence it has not been able to overtake 
the region of interest. 

As the angle of incidence is increased, the first limiting angle to be 
exceeded is «,. Figure 14b (plate 5) shows the pattern that is observed 
for this situation at € = 0:30. Since x,, has also been exceeded, it becomes 
necessary to investigate the extent to which the corner signal is contaminating 
the refraction pattern. In practice this was done by comparing many 
patterns taken in this same irregular region, at different shock strengths 
and with different settings of the back boundary, to separate out those 
effects which were boundary-dependent and hence not directly related to 
the exceeding of the x,-conditiont. This done, the following features of 
figure 14b remain attributable to the refraction process. 

The reflected wave RS is curved all the way into the intersection, and 
its strength increases very rapidly near the intersection. In its region of 
greatest strength, RS is followed by a sharp expansion region, which has 
the effect of ‘ peaking’ the reflected shock, the front of which is now com-, 
parable in strength to the incident shock! Curiously enough, none of the 
other regions or signals in the pattern are appreciably disturbed. Even the 
deflected interface, which is adjacent to the outflow from this anomalous 
reflection, is still radially straight from the intersection point. 

It is best next to consider the effect of exceeding «,, which is shown in 
figure 14c (plate 6) for € = 0-30. When adequate compensation for the 
corner effects has been made, in a similar manner to that described above, 
the following features of the pattern remain pertinent to the refraction 
process for this case. 
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A configuration much like a Mach reflection exists on the air side of 
the interface, complete with a Mach stem, a curved, varying strength reflec- 
tion, and a slipstream from the triple point. ‘The reflected signal and the 
stem are both peaked near the triple point. ‘he deflected interface is 
now curved along its entire length and approaches quite close to the slip- 
stream as the two near the back plate. ‘Thereafter they both exhibit an 
almost turbulent curling. 

‘The etfect of exceeding x, is not very noticeable at the stronger shock 
strengths, and z; is imaginary there, so to study these we returned to the 
case € = 0-85. ‘The onset of the process is shown in figure 14d (plate 7). 
Note that several significant alterations have taken place. ‘The transmitted 
shock has degenerated, near the interface, into a continuous compression, 
the leading edge of which is overrunning the intersection point of the 
incident shock slightly. A compression pulse has been retransmitted into 
region (1), ahead of the incident shock J. ‘This, in turn, deforms / near the 
interface. ‘lhe reflected shock and the sharp rarefaction region following 
it (cf. 14c, plate 6) have been distorted into a complex pulse, the leading 
edge of which has advanced up the back of / to the point where / intersects 
the retransmitted compression. 

As the angle of incidence increases still further, the effects just described 
become more distinct. For example, figure 14e (plate 8) displays the 
situation for% — x;. In this case, the inflow in region (5) on the transmitted 
wave would not be even sonic if 7’ did not run out ahead. ‘Two differences 
may be noted between this pattern and that for « > z;,. ‘he transmitted 
shock has now degenerated completely into a compression front at the 
interface. Also, the reflection is a cleaner shock, and the rarefaction 
following it is much less steep, than those appearing in the previous case. 
The air/COy, irregular patterns 

The irregular refraction region for the air/CO, problem is considerably 
less complex. ‘The only pertinent limiting angles are x, and x, (and x, 
for € 0-37), and, as may be seen in figure 13, the curves for these, along 
with that for x,,, are crowded together in a narrow band over most of the 
range of values of €. Nowhere is the angular interval between «, and «, 
(or «,) large enough to permit an observable development of the irregular 
ettect related solely to the former, before it is obscured by the rather prom- 
inent ettect arising from the latter. Consequently, it is possible to observe 
only one type of irregular pattern for this problem—one that exhibits 
primarily the effect of exceeding the x,-condition. ‘This isa rather interesting 
case, however, since it displays the transition of the reflected rarefaction 
from a supersonic to subsonic phenomenon. 

An example of the regular refraction pattern that exists at values of 
just below «, is shown in figure 14f (plate 9). ‘he irregular pattern that 
appears when ~, is exceeded follows in figure 14g (plate 10). ‘The principal 
changes to be noted concern the reflected rarefaction which, by the definition 


of x, now finds itself in a subsonic region. RR is no longer a centred wave, 
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but its leading edge has advanced into region (2), having climbed a short 
distance up the back of the incident shock /. (‘This intersection of the 
leading edge of RR on J subtends an angle at the corner which is approxi- 
mately equal to z,.) “That portion of / that is overtaken by this rarefaction 
edge is weakened by it, and consequentiy curved backward slightly, causing 
it to make a smaller angle with the interface. ‘The back edge of the rare- 
faction remains fixed at the intersection point. With these exceptions, the 
remainder of the fields (2) and (3) are uniform, as in regular refraction. 

As « is increased further into the irregular range, the situation described 
above changes in magnitude only, with no other noticeable ettects appearing. 

The pattern shown in figure 14g was obtained with the back plate set 
normal to the transmitted shock. ‘Te illustrate how the corner effect from 
a less favourable boundary arrangement can confuse the refraction process, 
consider figure 14h (plate 11), which was obtained with the back plate set 
normal to the incident shock. In addition to the effects observed in the 
preceding pattern, there now appears a compression region behind the 
reflection—a feature which was totally absent in the other cases. Apparently 
this results from the draining of the excess pressure behind 7°R across the 
interface into the upper region (3), in which the compression can travel fast 
enough to reach the reflected wave. Some of this effect even extends 
back across the interface into the region (4) behind 7, making it also non- 
uniform. ‘The relief of the pressure behind 7°R is evidenced by the strong 
rarefaction following it. From another point of view, the shock TR can be 
considered to be refracting, irregularly, at the CO,/air interface, producing 
a configuration not unlike that seen for air/CH, at glancing incidence (cf. 
figure 14e, plate 8). 


The interpretation of the irregular patterns 

From the experiments we know that each irregular refraction pattern ts 
a pseudo-stationary equilibrium configuration; that is, the pattern grows 
linearly with time+. As such, the process can be treated as a steady state 
when expressed in distance/time coordinates. A formal theory of this 
interaction would have to derive the steady-state configuration from basic 
principles and initial conditions, just as was done for regular refraction. 
In this case, though, the assumed models would need to be quite complex, 
and the resulting computations correspondingly complicated. However, 
since we are dealing with a pseudo-stationary problem, some qualitative 
predictions can be made from perturbation considerations. From the 
experiments we have seen that the various configurations that develop as 
the angle of incidence is increased are related to one another in a continuous 
fashion, forming, as it were, a sequence, or evolution of patterns, rather 
than exhibiting any abrupt changes in detail. Hence, it should be possible 
to infer the nature of any one of them from a knowledge of an adjacent 
pattern, and consideration of the flow conditions that necessitate the change. 

As an illustration, consider first the one observable irregular pattern 
for air/CO,. ‘The regular pattern observed for « just below «, is shown in 
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figure 14f (plate 9). Now, perturb this stable configuration slightly by 
increasing « until it is larger than x, and then readjust the pattern accordingly 
until a new equilibrium condition is achieved, as follows. 

When ~ exceeds ~,, the gas flow in region (2) becomes subsonic relative 
to the reflected rarefaction RR. ‘Therefore, the Prandtl-Meyer wave, 
which requires sonic flow normal to all its radial characteristics, can no 
longer exist as such. Instead, the leading edge of the rarefaction is free to 
run ahead into region (2). Consequently, this front, remaining straight 
but no longer centred at the intersection point -Y, will overrun part of the 
incident shock /. If « is not too much greater than x,, the expansion will 
still be sufficient to return the flow to sonic velocity behind RR, and the 
trailing edge will retain its centred position. 

‘That portion of / near the interface that has been overrun by the rare- 
faction will now be weaker and hence more inclined to the incoming flow 
there. ‘lhe ettect, thus, is to curve / from the interface up to the inter- 
section with the rarefaction front, causing it to be incident at a smaller 
angle on the interface. ‘he pattern just constructed is essentially that 
observed experimentally (figure 14g, plate 10). 

A similar line of reasoning can be applied to the series of air, CH, irregular 
refraction patterns. ‘lhe first critical condition is associated with the 
z,-limit, the violation of which results in the most subtle, yet probably the 
most significant, of all the irregular etfects. If we make ~ slightly larger 
than x,, We impose upon the reflected shock RS and the transmitted shock 
7’ incompatible pressure and deflection conditions. ‘he essence of the 
situation is that, if we required RS and 7 to meet the pressure requirements 
(p; = p,), the outflows in regions (3) and (4) would be convergent. 

A clue to the resolution of this paradox can be taken from a distantly 
related situation. In the problem of the supersonic aerofoil, for a given 
flight velocity, there is a limiting body angle beyond which a straight bow 
shock cannot deflect the flow parallel to the surfaces. For cases above this 
limit, the observed behaviour is that the bow shock becomes stronger near 
the leading edge, and less inclined to the flow there. ‘This stronger shock 
produces a subsonic rather than supersonic region behind it, in which 
subsonic region the necessary flow deflection is taken care of continuously. 
(Strictly speaking, it cannot be continuous right in the corner between the 
shock and the boundary. Here the ‘continuous’ variation must be 
condensed into an infinitesimal region—a common type of subsonic 
singularity, cf. subsonic flow over a sharp corner.) ‘This subsonic region is 
clearly an expansion, and is most severe where the shock is strongest. If 
the body angle is increased yet further, the bow wave detaches from the 
leading edge, and advances into the oncoming flow. 

The reflected shock in the supercritical refraction pattern apparently 
cures its deflection dithculty by essentially the same process, that is, it 
becomes stronger and less inclined to the inflow, and in so doing produces 
a subsonic region behind it which handles the deflection and presssure 
adjustment in a continuous fashion (see figure 15). 
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This mechanism can only exist as shown over a small range of a As « 
is increased further, the gas flow in region (2) relative to the reflection 
becomes less supersonic, and this reflected strong shock must become more 
normal to it to maintain its strength. Beyond the situation where it is just 
normal to the inflow, it has no other course than to leave the intersection 
point and advance into region (2), much as the bow wave detaches from the 
aerofoil. ‘This, of course, is just the « ,-limit, beyond which a new irregular 


pattern occurs. 
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Fig re 15. The x -adjustment. The reflected shock R Decomes stronger near the 


interface, and is followed immediately by a steep subsonic rarefaction zone, 
which completes the necessary flow deflection and pressure adjustment 


That the process described above is far from insignificant is clear from 
the observed patterns of figure 14b (plate 5) which show the great increase 
in strength of RS near the interface, and the sharp rarefaction region behind 
it. It is also impressive to compare the measured magnitude of this 
shock front near the interface with the reflected wave strengths in regular 
refractions (see figure 16). 

It would seem only reasonable that some process like that developed 
above would likewise be expected in connection with the transmitted shock— 
to help out, as it were, with the flow deflection problem. Yet it remains 
a Curious point in this case that there is no observable participation of the 
transmitted shock, and the flow behind it, in this readjustment process. 
As far as can be seen from the interferograms, 7 remains straight and 
unchanged in strength, and the field behind it remains very nearly uniform 
even though a rather drastic process is taking place in the adjacent region. 
Even the deflected interface seems totally unperturbed. 

The x,-adjustment for the air'CH, problem follows directly from the 
x,-configuration just discussed. When « exceeds «,, the gas flow in region 
(2) is no longer fast enough with respect to RS, even when normal to it, to 
maintain it at the required strength. As in the corresponding case for 
air/CQ,, RS will then advance into region (2) at a speed sufficient to maintain 
the necessary gas inflow on itself. In so doing, RS overtakes a portion of J, 
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thereby strengthening it, and establishing a pattern quite similar to a Mach 
reflection (cf. figure 14c, plate 6). ‘he difference in entropy change for 
the gas crossing / and then RS, and that crossing the stem, necessitates a 
slipstream from the triple point. ‘This slipstream is observed experimentally 
and is valuable in the interpretation of the patterns, since it assumes the 
direction of the local flow velocity. 
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Figure 16. Strength of reflected wave front (1 €) vs angles of incidence (a) for 
air CH, with € 0°30. Solid line follows regular refraction solutions 
Crosses mark experimentally observed strengths. Note the extreme maximum 
immediately following «,. 


The postulated pattern is quite eficient in relieving the critical conditions 
that necessitated the alterations. ‘lhe stem shock is normal to the interface 
and strong enough to leave a subsonic region behind it, the combination of 
which handles the deflection of the flow near the interface quite adequately. 
‘The reflected wave is now somewhat removed from the interface, and no 
longer need deflect the flow behind it to be exactly parallel to the outflow 
trom 7’, since there is a subsonic transition between regions (3) and (4). 
That this is indeed the case may be seen in figure 14c (plate 6) where the 
slipstream from the triple point, which follows the local flow direction, is 
clearly not parallel to the deflected interface. ‘The less strenuous deflection 
requirements on RS are reflected in a lessening of its strength, quite sharply, 
as x increases (cf. figure 16). 

For € = 0-30, the limit x, is closely followed by z,, which then is the 
last pertinent critical angle. However, for the weak shocks (& = 0-85), 
x, precedes x, by 13° and therefore produces a much more clearly defined 


effect. For this reason we will construct the x, pattern at € = 0-85, with the 
understanding that a similar situation must occur for the stronger shocks 
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also. At ¢ = 0-85, the previous pattern was that developed for « > «, 
(cf. figure 14b, plate 5). As & is increased past x,, the transmitted shock 
(which by definition of x, now makes an angle of x” = 90° with the interface) 
finds itself in a region where the inflow velocity is too small for its shock 
strength. ‘lo remedy this, the shock runs ahead in region (5), no longer 
maintaining contact with the intersection point (cf. figure 17a). In so 
doing, however, a pressure difference is set up across the interface between 
the intersections of J and 7. ‘This of course requires a compression signal 
to be sent upward into (1) ahead of /, and, simultaneously, a rarefaction 
downward behind 7 (figure 17b). ‘The effect of the compressed region 
ahead of J is to curve / forward, causing it to make a larger angle of incidence 
with the interface (figure 17c). ‘The region behind this curved portion of J 





Figure 17. The x,-adjustments. (a) The transmitted shock T leaves the intersection 
point and runs ahead into region (5). (b) This causes a compression front to 
propagate into region (1), ahead of the incident shock /, and simultaneously, 
a rarefaction front to advance down into region (4) behind 7. (c) Conse- 
quently, 7 is weakened near the interface, while / travels more rapidly into 
the gas passed over by the compression. (d) The region behind the curved 
portion of J is now subsonic, hence the reflected shock RS runs ahead into 
region (2), until it reaches supersonic flow again, that is, behind the straight 


portion of J. 


is no longer uniform. In particular, it is no longer supersonic with respect 
to the reflection, since the angle of incidence of this stem is past the « ,-limit. 
‘The process adopted in a case like this, as we have already seen in other 
cases, is for the reflection to advance into (2), climbing the back of / to the 
point where an adequate flow normal to itself is available, which in this case 
should be near the start of the curvature of / (figure 17 d). 

The rarefaction that advances into (4) will weaken 7 somewhat, thereby 
reducing the inflow requirements slightly. ‘This feature tends to stabilize 
the situation, in a sense restraining 7 from getting too far ahead, ‘The 
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equilibrium situation finds 7—a finite shock, followed by a region of 
compression to p,—running just slightly ahead of the intersection point 
of J. The pattern that actually occurs for «~- «, (figure 14d, plate 7) 
corresponds to this equilibrium configuration. 

If x is increased further, the inflow velocity in region (5) becomes 
subsonic (x ~ «;), in which case there is nothing to prevent 7' from degenera- 
ting completely into a continuous compression, the front of which will run 
ahead with velocity a;. The resulting other adjustments will be qualitatively 
the same as those described above for the «,; case, although the interaction 
region will now be much larger. ‘This is the configuration displayed in 
figure 14e (plate 8). 

The process cannot proceed to this extreme in the stronger shock case, 
since the velocity of the configuration along the interface is more rapid, 
and hence the inflow deficit on 7 is relatively less significantf. 

Figure 18 (plate 12) shows the z, effect for € = 0-30. Note that even at 
this large « (65°), 7 is only slightly ahead of the incident stem. One of the 
most interesting features of this picture is that the stem has been deformed 
by the precursor compression in such a way that it is actually incident at 
a > 90°! As a consequence of this, the flow behind it temporarily has an 
upward component, and the interface is seen to be deflected up above its 
original position for a short distance. 

It is clear from even such a sketchy analysis as that outlined above that 
the irregular patterns have a significance beyond their specific application 
to the refraction problem. In addition to providing some basis for empirical 
treatment of the irregular refraction regime, they are an abundant source 
of information on various shock interactions and flow field configurations 
which, in many cases, would be difficult to observe by more ad hoc experi- 
ments. In a sense, each irregular refraction process acts as a matrix in 
which are embedded, for our convenient observation, interactions between 
shock and shock, shock and rarefaction, shock and contact surface, shock 
and subsonic region, and many others, all existing in their ‘natural’ state. 
Such interactions are the fundamental elements of most complex shock 
processes, and whatever understanding of them can be acquired from this 
source should aid in the comprehension of other involved processes in 
which they also participate. 

In this survey we have included at least one example of every type of 
irregular pattern we have observed, and it now appears that extension of the 
study to other gas combinations would not produce any configurations 
fundamentally different from those shown here. Such an extension 
would not be without value, however, since the order of occurrence of the 
effects and the quantitative nature of the various separate interactions can 
be significantly varied by the proper choice of gases. In fact, it is this 
possibility of using different gas combinations that promises the necessary 
breadth of range for the use of irregular refraction patterns as display cases 


for aerodynamic interactions, 
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Supersonic flow past a bluff body with a detached shock 
Part II Axisymmetrical body 


By W. CHESTER 
Department of Mathematics, University of Bristol 


(Received 20 August 1956) 


SUMMARY 


This paper considers the flow at high Mach number behind the 
curved shock formed when a supersonic stream impinges on an 
axisymmetrical body with a rounded nose. 

The solution is obtained as a double expansion in 

6 = (y—I)/(y+)), 
where y is the adiabatic index, and M~; the expansion is developed 
to within terms of order (6+ M~?). 

Expressions are obtained for the distance between the body 
and the shock, the radius of curvature of the shock compared with 
that of the body, and the pressure distribution on the body. 


INTRODUCTION AND FUNDAMENTAL EQUATIONS 

The technique evolved in Part I (Chester 1956) is used here to investigate 
supersonic flow past a blunt-nosed axisymmetrical body. ‘The exposition 
is less detailed, for the ideas involved are essentially those of the previous 
paper. 

The common axis of the body and of the shock wave formed ahead of it 
is in the direction of the uniform flow upstream of the shock. ‘The 
x-coordinate is measured along this axis from an origin at the vertex of 
the shock. ‘The radial coordinate is denoted by y. 

Upstream of the shock the velocity, pressure and density are denoted 
by (V,0), Po, po respectively. Corresponding quantities in the region of 
disturbed flow behind the shock are (Vu, Vz), pp Vp and py p(y + 1)/(y— 1); 
as before, y is the adiabatic index. : 

The equations of conservation of mass, momentum and entropy are, 


with (y—1)/(y +1) = 8, 
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The first of these equations implies the existence of a stream function % 
such that dy, = puy, d, = —pvy. In the uniform flow ahead of the 
shock, % = }y? and is continuous across the shock. On the body, # = 0. 

With the independent variables transformed to (%, y), equations (1) 
become 
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pe” = fl), 
where f is an arbitrary function of #. 
We consider the particular case of a paraboloidal shock with unit radius 
of curvature at the nose. ‘The equation of the shock profile is then x = }y?, 
and the shock transition relations give the following boundary conditions 
to be satisfied on ys = $y?: 
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From these equations it follows that 
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From the third of equations (2) we also have 
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and, hence, the equation of the body is 
~hy? 1 7 
x = }y2+8 | (=) dp. (7) 
Jo \pvy 
As in Part I, it is not difficult to show that equation (6) includes 
automatically the stand-off distance between the body and the shock. 
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BASIC SOLUTION 
When 6 = M~ = 0, equation (4) becomes 
1 
” ~1)-2 
f(t) sai 1+ 2 +6 M ’ 
and equations (2) are easily solved using the boundary conditions (3). 
The results are 
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where 
A = [(2%)'?(1 + 2)!*— sinh-(2p)!? + y(1 + y?)!? + sinh-!y]. (9) 
The next step is to improve this approximation so as to obtain expressions 
for the velocity components which are uniformly valid as Y—>0. The 
expressions for u/v, p and p do not require modification; and so we have, 
correct to the first order of small quantities, 


£ = fp’! = fld){l +25 log p} 
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This expression for p/p is substituted in the second of equations 
(2) to give 
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where d = 6+ M~, 

Equation (11) is now solved in conjunction with u = vy, and all terms 
which are uniformly O(d) are discarded. ‘The following relations are then 
obtained : 
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FIRST-ORDER APPROXIMATION 
From equations (8) and (12), we can now obtain the following expression 
for 5/pv, which is uniformly correct as far as terms which are O(d): 
5 —-2y(1 + y?)2(d + 2M). 
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We can then obtain a first approximation to the equation of the body 
by the substitution of (14) in (7). The result is 
ety? dds 
x = dy?+2d(1+y?)? 
v= gy ted ty) | (apy + 2p A 
2y(1+y?)1?2 2d3?2Bu2(1 + a 
~ y(1+y2)¥2 4 sinh-y 








(15) 





+ M-*(1+y2)logd 


\v(1 +9?)!? + sinh“y 
Note that equation (15) contains terms which are O(d)??, and is in fact 
correct to within O(d?). 
To obtain the velocity components we must first calculate u/v from 
(6) which gives, correct to within O(d?), 
u/v = y(1+g) + 2dypy-*(1 + y?)"2(1 + 2h)"?(26 + dB)??, (16) 
where 


g = —2dyy-2(1 + y2)¥2(1 + 2p)! Nap alte AZ 4 
7, \ ~. 
M- 





(Qi + 2p) * 
H+ y2)Plogi2y(1-+y)!24) - 
9) 4 4y2)2 
- y(1 aN . et ly las sian (24) | (17) 
The form of the expression for u/v anticipates the fact that, near y = 0, 
u~ 2d¢/py? ~ 2dib/y? 
(the argument is similar to the corresponding one in Part I). Hence we 
may write 





u= U+2dp/y", (18) 
where U is O(y?) near y= 0. Equation (16) then gives 
U/v = y(1 +8); (19) 


and the singular term in (16) has now disappeared. 
When these relations are substituted in (11) and the resulting equation 
solved for U, we find that 


aa + y(1 + 2f)-22(1 + ayaa( 1 4a)x 
~~ re “T+y 


x | 26+ 2¢d-+ 2g ‘log = aaeny | (20) 


from which u and wv are easily obtained using equations (14) and (15) 
respectively. 
Finally the pressure is obtained by integration of the equation 

op leu 

op y oy 
using the previously calculated expression for u and the boundary condition 
on = 4y* obtained from the last of equations (3). Although the variation 
of p with both % and y is required before the next approximation can be 
found, the complete expression is somewhat involved and we quote here 
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only the expansion for small values of y and %Y = 0. ‘The basic terms are 
given exactly, but the rest of the terms have errors which are O(y®): 


ma ie sm, ee oy | FH ,.. Sef 
p = (Lestye] a+ yies SP +d{- 5-J559°+ ism") 


7 13741 8d\3% 
1 rr a8 PW devs 
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SECOND-ORDER APPROXIMATION 





The precautions referred to in Part I also apply here; otherwise the 
process is straightforward, and only the final results are given. Apart 
from the basic terms, fourth and higher powers of y have been ignored. 

The equation of the body is 
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and the pressure on the body is given by 
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DISCUSSION OF RESULTS 

Although, for reasons given in Part I, the analysis is carried out for the 
case of a shock having unit radius of curvature at the nose, the following 
results have been corrected to apply to a body having unit radius of curvature 
at the nose. 

Figures 1 and 2 show respectively the radius of curvature of the shock 
and the distance between the body and the shock along the axis of symmetry 
for various values of d. ‘The isolated point in each figure is an experimental 
result given by Oliver (1956). The curves are calculated from equation (22), 
except for the broken parts which are simply extrapolations through values 
of d for which the convergence of (22) may not be satisfactory. 

A comparison of corresponding expressions in Part I and II shows that 
the convergence is poorer in the axisymmetrical case. Nevertheless, one 
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or two general conclusions may be drawn. As in the two-dimensional 
problem, the stand-off distance seems to be virtually independent of M~, 
except in so far as this parameter influences d, and the single curve in 
figure 2 covers all cases. 

As would be expected, under similar conditions both the radius of the 
shock and the stand-off distance are smaller in the axisymmetrical problem; 
a rough guide for the distance between the body and the shock is 2d in the 
two-dimensional case, and 3d in the axisymmetrical case. 





























Figure 1. The radius of curvature of the shock on the axis of symmetry for various 
values of d(= (y—1)/(y+1)+M~?). The unit of length is the radius of 
curvature of the body on the axis of symmetry. The cross denotes an 


experimental result. 
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Figure 2. The distance between the body and the shock along the axis of symmetry 
for various values of d( (y—1)/(y +1) M-?). The unit of length is the 
radius of curvature of the body on the axis of symmetry. The cross denotes 


an experimental result. 


Because of the large factor multiplying d* (compared with the multiplier 
of d) in the coefficient of y? in equation (23), the expression for p cannot 
be regarded as reliable save for very small values of d in the neighbourhood 
of 0-01. 

When d= M~? = 0, the pressure distribution on the body is known 
exactly, and the ordinate of the sonic point is 0-680 compared with the 
value 0-629 in the two-dimensional case. According to calculations based 
on equation (23), the initial tendency is for the ordinate to decrease as d 
increases (and M-? = 0) in contrast to an increase in the two-dimensional 
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case (see Part I, figure 3), the ordinates agreeing in the two cases for a value 
of d around 0-02. But up to this value the decrease is only 0-3°%, and the 
most that can be said is that initially there is no appreciable variation in the 
position of the sonic point. 

Since the pressures on the body at the stagnation point and the sonic 
point are functions only of d and M-? (see equations (29) and (30) of Part I), 
it follows that, with similar conditions, the pressures are the same at these 
corresponding points in the two-dimensional and axisymmetrical cases. 
Hence we may deduce that for sufficiently small values of d, the pressure 
distribution in the two cases do not differ appreciably in the neighbourhood 
of the nose. 
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An inequality concerning the production of vorticity 
in isotropic turbulence 
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SUMMARY 

An inequality is demonstrated involving the rate of production 
S of mean-square vorticity in isotropic turbulence and a factor y 
which may be said to allow for intermittency or for the non- 
vanishing of fourth-order cumulants. An extreme state, corres- 
ponding to equality of this relationship, occurs if S = 0-64 and 
y = 2:14. The experimental values are S=0-4 and y= 4. 
Another kinematical relation shows that the mechanism of vorticity 
production resembles collision between fluid particles rather than 
the swirling of contracting jets. 


1. INTRODUCTION 


When the motion of a fluid is turbulent, the fluctuations of the vorticity 
are of major importance. For isotropic conditions they are directly related 
to energy dissipation by viscosity. They aiso show a sufficient degree of 
isotropy and homogeneity in general to justify the corresponding simplific- 
ations of the theory. 

These fluctuations obey a system of equations in which the pressure 
does not appear and thereby the degree of complication is reduced. Accord- 
ing to these equations, viscosity tends to eliminate vorticity, by a well- 
understood process. ‘These equations also show that the mean-square 
vorticity can be increased by inertial effects. This aspect of the problem is 
still somewhat mysterious. What determines this production of vorticity 
and how does it occur? 

The problem has been stated in terms of correlation functions, but this 
procedure leads to a number of unknowns always larger than the number of 
equations. 

The Fourier transforms of the correlations can be introduced, amounting 
to a spectral analysis of a sample of turbulent flow. They lead to a relation 
between vorticity production and the rate of energy transfer from large 
eddies to small eddies. Each hypothesis concerning energy transfer has 
thus been followed by an evaluation of the rate of production of vorticity. 

In this paper, we investigate the largest rate of production of vorticity 
compatible with the requirements of isotropy, homogeneity and incom- 
pressibility. We find an upper bound for the rate of vorticity production, 
although it depends upon a fourth-order mean value and therefore does 
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not constitute an absolute maximum. Comparison with old and new 
experimental data indicates that a turbulent flow produces vorticity at 
about 50°, of this limiting rate. It is quite possible that dynamical effects 
are responsible for this factor and that we are confronted by an extreme 
state of affairs, characteristic of turbulent flows. 


2. FUNDAMENTAL INEQUALITY 
We intend to establish an inequality valid at every point of an incom- 
pressible flow and at every instant. We shall use cartesian coordinates x,, 
denote the velocity vector by u,, and write u,; for cu,/cx;. The condition 
of incompressibility, with the usual convention about summation of 
repeated indices, is 
We can write u,; = 5,;+7,;, where s,; is the symmetric rate of strain 
tensor and r,; the antisymmetric vorticity tensor. At any point and instant, 
we can find the principal axes of s,; and, on referring to these axes, we have: 
| a A —B 
Mey = | —C b A : (2) 


B —A c | 

where a, b, c are the eigenvalues of s,, and A, B, C are half the vorticity 
components in the directions of the principal axes. In the sequel we shall 
designate the directions associated with a, b, c as the a-axis, b-axis and c-axis. 
Furthermore, we shall adopt the convention 

a| > |b| > Ic. (3) 
By (1) we have 

a+b+c = 0, (4) 
and consequently 6 and c must be of the same sign and a and 6 of opposite 
sign. 

We shall make use of the general Cauchy inequality (Hardy et al. 1952) 

(a* 4 53 me c3)? ; (a? + b2 a c*)(a* + }4 +c*) (5) 
with equality if a=b=c. We shall also use the following identities, 
obtained by raising both sides of (4) to powers 2, 3 and 4, 

a®+b?+c? = —2(ab+bce+ca), 
a+63+c? = 3abc, (6) 
aé+bhi+e 1(a? + b? + c?)?. 

Values of a, b and c for which the two sides of (5) are nearly equal are 
excluded by the incompressibility condition. Indeed, when a, 6 and c 
satisty (4), the ratio of the left-hand side to the right-hand side of (5) has a 
maximum value of } at —a = 2b = 2c. This allows a factor } to be inserted 


II 


on the right-hand side of (5), and so, with the use of (6), we have 


abc ae (a? + b? +. c?)3?, (7) 
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It is known that the set of eigenvalues of a tensor is invariant to trans- 
formations of the coordinates; therefore each term of (7) is invariant. 
Indeed, each term can be expressed by invariants such as §,;5;; OF 5); 85). 8,;- 
Consequently (7) remains valid if we pass to a fixed system of coordinates, 
independent of the orientation of the a-, b-, and c-axes. We now consider 
an ensemble of points of the fluid, each with its own set of eigenvalues 
with respect to these fixed axes, and average both sides of (7). We can use 
a space, a time or an ensemble average. ‘Then, with ¢ ) indicating the 


average, we have 





| <abc)| ((a? + b? + c?)8?), (8) 


The right-hand side of (8) cannot be measured by conventional methods 
and comparison with experimental results is not possible unless we introduce 
a modified inequality. We can write a single inequality involving N sets 
of eigenvalues, viz. 
N 2 N N 
(5 at+o3+ c) (3 a2-+63+ <)( S af +b5+c4) (9) 
a=1 1 1 
By taking (4) and (6) into account, and introducing the average as a limit 
for large N, the inequality (9) becomes 
1 

3y3 

The right-hand side of (8) is always less than or equal to the right-hand 
side of (10), with equality if a, 6 and c do not fluctuate. If, in addition, 
a= —2b = —2c, the two sides of (10) are equal. In view of (6), the ratio 
of right-hand sides of (8) and (10) depends only upon the probability 
distribution of a?+6?+c?. In the case of a Gaussian probability distri- 
bution of a? + 6?+c?, the right-hand side of (8) is only 8%, smaller than that 
of (10). In general, this ratio is close to unity, unless the probability 
distribution has a pronounced peak at the origin. Thus, the inequality (10) 
is not much weaker than (8). 


36 





{abc )| (a? + b?+c?)'2(at+b4+ ct), (10) 


3. RELATIONS WITH MEASURABLE QUANTITIES 

The various invariants occurring in (10) can be measured with a single 
hot-wire anemometer. ‘This has been done in the turbulence produced 
by a grid obstructing a parallel flow of air. The hot-wire is located 
sufficiently far down-stream to give some guarantee of isotropic conditions, 
and its orientation is normal to the mean flow. ‘The turbulence passes by 
the wire with a mean velocity about one hundred times larger than the 
velocity fluctuations, and consequently the time derivative of the hot-wire 
signal corresponds reasonably well to the space derivative of the velocity. 
With index 1 denoting the direction of mean flow, this means that the 
differentiated signal is proportional to wu, ,. 

Let us find the relation between u, , and the invariants of (10). At any 
particular instant, the hot-wire responds to a fluid particle whose principal 
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axes of rate of strain can be specified. Let ¢ be the angle between the 
l-axis and the a-axis (latitude) and ¢ be the angle between the c-axis and a 
direction normal to the a-axis and 1-axis (longitude). It then follows from 
tensor calculus that 

Uy, = acos*d +bcos* sin’ +c sin*y sin?¢. (11) 
The vorticity does not contribute to this particular signal, but this is peculiar 
to the component uw, ,, 

In the course of time the hot-wire is coincident with different fluid 
particles and the angles ¢ and @ vary as well as the eigenvalues a, b,c. The 
eigenvalues are absolute quantities. Consequently, if the turbulence is 
isotropic, the angles must be statistically independent of the eigenvalues. 
To clarify this point let us consider all fluid particles for which a = 1, 
b= —0-7, c= —0-3. If, in the average over these particles, we find a 
definite orientation of the principal axis, we can use the numbers 1, — 0-7, 
—0-3 to specify a direction. Such a direction would remain invariant to 
a rotation of the coordinates and this would reveal a lack of isotropy. 
Therefore, we must admit that, for any set of values of a, 5, c, the orientation 
of the principal axis are random; that is, we have statistical independence 
of eigenvalues and orientations. 

With (1/47)sin ¢ d¢ dis as the probability of finding ¢ and u inside a small 
solid angle, and statistical independence, we find 


‘ 
(ui) = is (a? +b? +c?), (12) 
The same method can be used to demonstrate 
* 24 
(ui) = i05 <abc), (13) 
8 
{ui,) = ins (a*+b*+c*). (14) 


Some of these results have been derived by others from correlation functions. 
When the following non-dimensional parameters are introduced, 
3 4 
Sis (uj) y = (Uin > 
os =— —s .ae? a ee mS 
(uz )°* «uy )° 
the inequality (10) can be expressed as: 
2 


=? 


(15) 


S ”, (16) 
For normally distributed velocities we have y= 3, and in that case 
S| < 0-756. 

In the extreme case —a = 2b = 2c = constant, we find by (12), (13) 
and (14), y = 2-14, |S|= 0-638. Proudman & Reid (1954) calculated 
the value of S with the assumptions of no viscosity, zero initial triple 
correlations and a relation between three-point quadruple velocity correla- 
tions and double correlations appropriate for normal probability distribution. 
They found S = 0-78, and by (16) this implies y > 3-2. This seems to 
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mean that the four-point quadruple correlations could not satisfy exactly a 
similar relation with double correlations, since this would imply y = 3. 


4, EXPERIMENTAL RESULTS 

The numerous measurements of S and y by Townsend (1947, 1948, 
1951) gave S=0-38, y=4. Later Stewart (1951) found values of S$ 
decreasing with increasing Reynolds number of the turbulence, in the 
range S = 0-5 to0-3. From the measurements of (0?u,/éx?) by Liepmann 
et al. (1951) one can find S, by making use of the vorticity equation. The 
resulting values of S increase from 0-4 to 0-5 with increasing Reynolds 
numbers. 

The author measured S and y for a variety of meshes and with improved 
electronic devices in the low-turbulence wind tunnel of the University of 
Maryland. The hot-wire signal was carefully compensated for thermal 
inertia. A simple low-pass filter was used to avoid excessive electronic 
noise, and it was verified that the filter did not affect the measurements of 
mean cubes. 

To measure mean squares and mean cubes, I used chains of biased 
diodes. For mean fourth power, I first obtained the instantaneous square 
with a quarter-square multiplier using four chains of diodes. 

The first grids used were of 25cm mesh, with 38 and 66% open 
area and air velocities of 2 to 20 m/sec. ‘The wire was located at distances 
from 30 to 100 mesh lengths from the grid. ‘The results consistently 
indicated S=0-4+0-1, without systematic variation. 

With a grid of 0-6cm mesh, 21% open area and a wind of 4m/sec I 
found S=0-4+0-05 at distances of 20 to 150 mesh lengths. With a first 
grid of 25cm and 54% open area, followed at 125cm downstream by a 
second grid of 1:25cm and 76% open area, I found S=0-4 at a distance 
of 200cm from the first grid. When several strips of paper tape were 
added to the first grid the turbulence was considerably altered thereby, 
but S was found to be unchanged. 

In a crude pipe flow and in the wake behind a cylinder, I again found 
S=0:-4. The central part of a turbulent boundary layer gave the same 
result. 

Only a few measurements of y were made, and they indicate y = 4+0-5. 

The constancy of S is remarkable, and holds at Reynolds numbers 
lower than those for which Kolmogoroff’s theory might be expected to be 
valid. Comparison with the inequality (16) shows that the experimental 
value of S/y!? corresponds roughly to 50% of the maximum permissible 
value. With y = 4, the inequality requires |S| < 0-872. 


5. PRODUCTION OF VORTICITY 
We shall now introduce the dynamical equations. With additional 
indices to indicate successive derivatives, we write the Navier-Stokes 


equations as 
Ou,/dt + Uz U;,— Un, = —Pi/p, (17) 
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where v is the kinematical viscosity, P the pressure and p the density. 
After differentiating with respect to x; multiplication by u, ;, Summation 
in z and j, and averaging we obtain 
( 
> CU; 5 Uy 5 +2 Cy Uy 5 Uy, ej) + 2 (Uy 5 Uy UZ) — 


I 
—2v{u; 54; jc.) = — 5 kis Pig) (18) 


For the first term we can write, with reference to (2), 
Pa re 2, pra ( 
(u,;U,;) = (a? +6? +c?) +2(A?+ B?+C?). (19) 
We now consider the invariant (u,;u,;;). In homogeneous and incom- 
pressible flow, we have 


‘ Oo 
(u; ;Ujy;) = = {u,u;;) = 0, (20) 
a” . 

giving the following relation between the mean-square vorticity and the 
mean-square rate of strain: 

(a? +6?+c?) = 2(A?+B?+C?). (21) 
The invariant (u,u,;U;;;) also vanishes, either by homogeneity or by 
isotropy. Also, for the same reasons, ¢P,;u,;) = 0. The viscous term 
can be related to the vorticity or the rate of strain, and can also be written as 


CU, Uigke) = — (Uign Mig) (22) 
The invariant (u,;u,,u,;) can be written as (u,;u;,U,,), Where the 
sequence of indices must be noted. With the use of (2), it becomes 
(u,; Uj .U;,,) = (a2 +63 +03) — (aA? + bB*+cC?). (23) 
Let us consider the invariant (u,;u;,u,;), where the index sequence is 
now different. ‘This quantity vanishes because it can be expressed as the 
divergence of a vector; 


C \ 
(u; Uj UR) = om CU, Uj Uy — 3U; Uy, Uy.) = 0. (24) 
a 


Equation (24) is a relation imposed by homogeneity, and with the use of 
(2) it becomes 
(+63 +c2) = —3(aA?+bB?+cC?). (25) 
This relation is comparable with (21) in that it is purely kinematical. The 
right-hand side describes the stretching of vortex lines, as discussed by 
Taylor (1938). 
We can now use relations (19), (21), (22), (23), and (25) and rewrite the 
dynamical equation (18) as 


(A?+ B?+C?) = — abc) - fVCU; 5p U; jn) (26) 


The viscous term is always positive and indicates dissipation of vorticity. 
Production of mean-square vorticity occurs if <abc) < 0, which is equi- 
valent to S > 0. Our inequalities therefore specify a maximum rate of 
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production of vorticity with a limit determined by mean square and mean 
fourth powers of the eigenvalues. 

Let us examine the effect of intermittency of the flow. We will say 
that the flow is intermittent if u,; — 0 in a variety of large and well defined 
regions, occupying a fraction « of the space filled by the fluid. The distri- 
bution of these regions can be isotropic and homogeneous and they can 
eventually have a size comparable with the integral scale of the turbulence. 
The quantities S and y can be measured as previously or similar quantities 
S’ and y’ can be obtained by averaging only over those regions where 
u,; fluctuates. It follows from (15) that S’ = 1/(1—x)S and y’ = (1—a)y. 
Hence S’/yy’ is independent of «, and this means that the presence of y 
in (16) allows for possible intermittency of the flow. 

The factor y’ is related to quadruple four-point correlations, and one 
could perhaps assume that, within the regions of finite u, ;, the fourth order 
cumulants vanish. ‘This would give y’ = 3. Observation of turbulent 
flows suggests that y and a fortiori y’ are not very different from 3, and 
that u,; is mildly intermittent. ‘This means that S could not exceed a 
value of about 0-9 without serious intermittency or without non-zero 
fourth order cumulants. 


6. A REMARK ON VORTICITY PRODUCTION 
From (6) and (25), we have the important relation 
—<abc) = (aA? +bB*+cC?). 
In a different form, this relation was first written by ‘Townsend (1951). 
We must remember that 4 and ¢ are of the same sign and that a is the largest 
eigenvalue. ‘hus abc has the sign of a and, for production of vorticity, 
we must have a predominance of points at which a < 0. 

If a < 0, the flow behaves locally as a jet parallel to the a-axis and 
impinging on a wall represented by the 4- and c-axes. ‘The dynamic 
equations show that A (the vorticity component normal to the wall) is 
attenuated, and that B and C are amplified. If this situation lasts long 
enough, we shall find, for any initial vorticity, that a4?+bB?+cC? > 0. 
This situation is therefore such that both sides of (27) are positive. If 
a > 0, the flow behaves locally as if it were entering a contraction cone. 
A small fluid sphere is deformed into a cigar-shaped ellipsoid. Only A is 
amplified, and if this situation lasts long enough, it will contribute positively 
to the right-hand side of (27) and negatively to the left-hand side. 

It is clear, therefore, that production of vorticity is associated essentially 
with a < 0 and amplification of B and C. This suggests that the most 
important processes associated with production of vorticity and energy 
transfer resemble a jet collision and not the swirling of a contracting jet. 


7. REMARK ON CONTINUITY 


So far, we have not mentioned any requirement of continuity, and the 
existence of second derivatives of the velocity has been assumed implicitly. 
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This question is of major interest, especially in the case of very small 
viscosity. Let us consider the extreme case of our inequality, that is, 
a= —2b= —2c=constant. Equations (21) and (25) impose some 
restrictions on the vorticity, and it turns out that they involve only <A?) 
and (B?+C*). This leaves a great freedom in the choice of the vorticity 
and the orientation of the principal axes. However, this does not guarantee 
that we can connect smoothly the values of u,; at a point with its values at 
some other point. When the viscosity is very small the momentum equation 
can conceivably include large values of w,; ;.,, but there may be a limit to the 
type and magnitude of the viscous terms. ‘This leads to the following 
question: Is it possible to construct a velocity field with maximum (abc), 
isotropic and homogeneous properties and such that it satisfies the momen- 
tum equations? In the limit of vanishing viscosity discontinuities of 
u, ; occur, and we can ask whether the condition a = —2b = —2c = constant 
makes the flow everywhere discontinuous or whether we can choose the 
other components of u,; in such a way that the discontinuities occur only 
along isolated surfaces (shear layers). ‘The answer to such questions could 
very well lead to a lower maximum value of S and eventually close the gap 
between experiments and theory. 
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viscous fluid between concentric spheres 
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SUMMARY 

‘Two concentric spheres are supposed to rotate about the same 
axis with almost the same angular velocity, so that the viscous 
stresses over the surfaces of the spheres induce a flow which may 
be represented by a small perturbation superimposed upon a rigid 
body rotation of the fluid as a whole. ‘The governing equations 
are therefore linearized in the magnitude of the perturbation, and 
it appears that the validity of this linearization is independent of the 
Reynolds number of the primary rotation. Attention is then 
restricted to the case in which the Reynolds number is large, the 
principal object of the note being to exemplify some of the pro- 
perties of rotating systems at large Reynolds numbers in terms of a 
particularly simple mathematical model. 

t is found that the cylindrical surface that touches the inner 
sphere (the axis being the axis of rotation) is a singular surface 
in which velocity gradients are very large. Everywhere outside 
this cylinder, the fluid rotates as a rigid body with the same angular 
velocity as the outer sphere. Inside the cylinder, the velocity 
distribution in the central (inviscid) core of the motion is shown to 
be determined by the velocity distribution in the boundary layers 
over the spheres, and explicit solutions are obtained for all these 
velocity distributions. ‘The mechanics of the cylindrical shear 
layer itself is also discussed, though no explicit solution is obtained 


in this case. 


1. INTRODUCTION 

In the literature relating to rotating systems in fluid dynamics, very 
little attention has been paid to problems in which the configuration of the 
boundaries is such that the flow may be represented by a small perturbation 
superimposed upon a rigid body rotation of the (incompressible) fluid as a 
whole. Indeed, the only relevant work of which I am aware is that of 
Squire (1953) and Stewartson (1953). Squire considered the flow due to 
a rotating disc when the fluid at a great distance from the disc is rotating as 
a rigid body about the same axis with almost the same angular velocity, and 
Stewartson considered the similar problem in which two discs rotate about 
the same axis with almost the same angular velocity. In both these 
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problems, the potential value of the basic approximation was, perhaps, some- 
what obscured by the fact that the general form of the exact solutions of the 
Navier-Stokes equations had previously been found by Batchelor (1951), 
The peculiar property of the general class of flows under discussion is, in 
fact, that the governing equations become linear in a manner that is in- 
dependent of the Reynolds number, and the much greater tractability of 
linear equations leads to the expectation that substantial analytical progress 
can be made in such problems. A discussion of many of the interesting 
properties of rotating systems is undoubtedly prohibited by the lineari- 
zation, but this is obviously not true of all such properties. In particular, 
the phenomenon of flow at large Reynolds numbers being governed by 
linear equations is of real interest, and may be expected to yield useful in- 
formation concerning the behaviour of boundary layers in rotating systems. 

The present note is concerned with the particular problem in which 
two concentric spheres rotate rapidly about the same axis with slightly 
different angular velocities. In such a system, it is clear that the viscous 
stresses over the boundaries must induce a secondary flow, and the prin- 
cipal purpose of the note is to examine the properties of this secondary flow 
when the Reynolds number is large. 


2. "THE GOVERNING EQUATIONS 

Let the radii of the inner and outer spheres be a and za, respectively, 
and let the corresponding angular velocities be 2 and Q2(1+€), where e is 
very small. Let (r’, 6, ¢) denote spherical polar coordinates in which the 
line 6 = 0 coincides with the axis of rotation, and let uv’, v’, w’, be the corres- 
ponding components of velocity (see figure 1). 

By symmetry, all dynamical variables must be independent of ¢, so 
that the incompressibility condition may be integrated in terms of a stream 
function w&’, where 





: 1 cw’ 7 1 cy’ 
i = Se2- on ae? a eae = (2.1) 
r’*sin@ cé r’sin@ er 
Then, writing 
, 
, x > 
vy = —. 22 
r’ sin 6’ (2.4) 


the Navier-Stokes equations of steady motion become (see Goldstein 1938, 


ch. 3) 
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1 o(d’, x’) 
ee eT RD y) 
and 7¥sinb a(r',6) vD*y’, (2.4) 
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The boundary conditions for equations (2.3) and (2.4) are : 
; 

, ow , , 7 . ¢ | 

yr’ = 4, —=y =(Q, x’ = Qa? sin??@ ; 
er 

‘ 

Aus’ (2.6) 
cy ’ , 9.9 .°_¢ } 

y’ = ad, —=y' =), xX = Q(1+6€)e7a* sin?6. | 
A J 


CID 2 (1+ €) 











Figure 1. Notation. 


r When « = 0, the exact solution of the problem is obviously the rigid 
body rotation represented by 
pf’ = 0, xX = Qr’ sin?é. 
Hence, when « is non-zero but small, it is reasonable to write 
Yb’ = eB Qu, xX = Qr’? sin?@ + ea? Qy, (2.7) 
and retain only terms of the first order in «. Thus, introducing the 


dimensionless variable ; b. 
r=r /@, (2.8) 


the dynamical equations (2-3) and (2-4) reduce to 


2 %.. _ lox i ~ 4 9) 
2(2 cos 4 ag sin a) RP yb, (2.9) 
-2( 3 cos 6— 7388) = 2D (2.10) 
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where D? is the same operator as (2.5), but with r’ replaced by r, and R is the 
Reynold’s number given by 


aQ 
R=—. ZA 
~ (2.11) 
‘The boundary conditions become 

i oe y=0, yx=0, 

or 

a , (2.12) 

ous re 
r= 4, —=7%=0, x = a’ s1n76. 

or 


‘The equations (2.9) and (2.10) form the basis of the present note. 

It is worth noting here that the approximation leading to (2.9) and (2.10) 
is somewhat analogous to the approximation introduced by Oseen (1910) in 
problems of uniform streaming past obstacles. For, in both cases, a linear 
approximation to the inertia terms in the Navier-Stokes equation is made. 
In Oseen’s case, however, it is an immediate (kinematic) consequence of 
the boundary conditions that the approximate form of the inertia terms is 
invalid in the neighbourhood of the obstacle, so that the approximation is 
spatially uniform only if the Reynolds number is very small. In the present 
case no such troubles seem to arise, since the boundary conditions are 
kinematically consistent with the linearization over the whole field of flow, 
and the Reynolds number of the primary motion remains a significant para- 
meter on which there is no immediately obvious restriction. Nevertheless, 
when the Reynolds number is very large, the possibility that intense velocity 
gradients 1n the secondary flow might invalidate the linearization must not 
be ignored. ‘This last point will be a matter for subsequent examination. 

Although the spherical polar coordinates are particularly suited to the 
geometry of the boundaries, many aspects of the natural dynamical be- 
haviour of the system at large Reynolds numbers are of a much more two- 
dimensional nature than the boundary conditions suggest. It will there- 
fore be useful to set out here the governing equations in terms of the cylin- 
drical coordinates (see figure 1) 


p=rsin§, C=rcosé. (2.13) 
Thus, (2.9) and (2.10) become 
2X _ = py, (2.14) 
ey 1 “ 
and ~2% px, (2.15) 
2 1a0( 
where DP oe = +e. (2.16) 


3. THE BOUNDARY LAYERS OVER THE SPHERES 


When the Reynolds number is very large, viscous forces must be 
negligible over the whole field of flow, except, possibly, in the neighbour- 
hood of certain singular surfaces. The cylindrical polar equations (2.14) 
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and (2.15) then show that, apart from such surfaces, the solutions must 
take the simple form 

b= (p), x = Xx(p); (3.1) 
and since it is clearly impossible to satisfy the viscous boundary conditions 
(2.12) with such functions, it follows that the singular surfaces must indeed 
exist in this problem. 

Now, it is natural in the first place to suppose that the singular surfaces 
take the form of conventional boundary layers over the two spherical bound- 
aries. ‘Thus, taking the case of the inner sphere first, the usual boundary 
layer approximations in the dynamical equations (2.9) and (2.10) lead to the 
equations 





Oo / 1 0 4); Ous 1 o2 / 
pe a oe —; — 2—cos@ = — cx (3.2) 
or R ort ( R or? ; 


In view of the large velocity gradients implied by the approximate equations 
(3.2), it is then proper to enquire whether these equations are consistent 
with the original linearization in «. Returning, therefore, to the exact 
dynamical equations (2.3) and (2.4), the neglected inertia terms in the first 
of these equations are 


| 1 0(s, D2) 


2y (oe hey... | 
e2al)? E= | 2x cos @— - ox sin@ +} — 


ale 


r?sin?@ | or r 00 r?sin?0  d(r,0) 
2D7s ( dus lows . 
———- — cos#— -— siné > |, 
r“sin°@ | or r o6 J 


which, on the assumption that a boundary layer exists at r = 1, reduce to 


ac a 4}. 93 Af a8 a) 98 
2aQ? 2cosd dx 1 (oped ap OM \ — 2 opar 
or © rc g 


sin2@ * Or sin \ 0@ dr? sin? or or? 











Thus, since we may assume 
b = O(R-), 4m 008), 
e/ér = O(R"?), 0/00 = O(1), 
the order of magnitude of the neglected terms is 
eaQ?R?, 
The order of magnitude of the inertia term retained in the first of equations 


(3.2), on the other hand, is 


cal? cos 8 or €aQ?R!?, 

Hence the neglected terms are smaller by a factor « than those retained, and 
the presence of large velocity gradients does not affect the basic approxi- 
mation. The argument may not, of course, apply near particular values of 
6, notably near the equator. Similar remarks apply to the other dynamical 
equation. 

The property of the equations (3.2) that is chiefly responsible for the 
success of this investigation is that they are both ordinary differential 
equations. This remarkable circumstance possibly embodies a point of 
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where D? is the same operator as (2.5), but with r’ replaced by r, and R is the 
Reynold’s number given by 





aQ 
R= . (2.11) 
Vv 
The boundary conditions become 
y= i, O chad x = 0, | 
or 
24 (2.12) 
Cw . r 
r= q, x =%=0, x = «sin?6. 
or 


‘The equations (2.9) and (2.10) form the basis of the present note. 

It is worth noting here that the approximation leading to (2.9) and (2.10) 
is somewhat analogous to the approximation introduced by Oseen (1910) in 
problems of uniform streaming past obstacles. For, in both cases, a linear 
approximation to the inertia terms in the Navier-Stokes equation is made. 
In Oseen’s case, however, it is an immediate (kinematic) consequence of 
the boundary conditions that the approximate form of the inertia terms is 
invalid in the neighbourhood of the obstacle, so that the approximation is 
spatially uniform only if the Reynolds number is very small. In the present 
case no such troubles seem to arise, since the boundary conditions are 
kinematically consistent with the linearization over the whole field of flow, 
and the Reynolds number of the primary motion remains a significant para- 
meter on which there is no immediately obvious restriction. Nevertheless, 
when the Reynolds number is very large, the possibility that intense velocity 
gradients in the secondary flow might invalidate the linearization must not 
be ignored. ‘This last point will be a matter for subsequent examination. 

Although the spherical polar coordinates are particularly suited to the 
geometry of the boundaries, many aspects of the natural dynamical be- 
haviour of the system at large Reynolds numbers are of a much more two- 
dimensional nature than the boundary conditions suggest. It will there- 
fore be useful to set out here the governing equations in terms of the cylin- 
drical coordinates (see ‘figure 1) 

p=rsin§, C=rcos@. (2.13) 


Thus, (2.9) and (2.10) become 





ox 1 
2—=— us 
21 RP ds, (2.14) 
' 5 1 D? 215 
anc “+o Bee (2.15) 
where D? = a : + = (2.16) 
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When the Reynolds number is very large, viscous forces must be 
negligible over the whole field of flow, except, possibly, in the neighbour- 
hood of certain singular surfaces. ‘The cylindrical polar equations (2.14) 
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and (2.15) then show that, apart from such surfaces, the solutions must 
take the simple form 

b= de), x = Xx(p); (3.1) 
and since it is clearly impossible to satisfy the viscous boundary conditions 
(2.12) with such functions, it follows that the singular surfaces must indeed 
exist in this problem. 

Now, it is natural in the first place to suppose that the singular surfaces 
take the form of conventional boundary layers over the two spherical bound- 
aries. ‘Thus, taking the case of the inner sphere first, the usual boundary 
layer approximations in the dynamical equations (2.9) and (2.10) lead to the 
equations 

4 p 7 a2 

22x cos@ = rey —2 - cos@ = E. at. (3.2) 

or R or'’ or R or? ‘ 
In view of the large velocity gradients implied by the approximate equations 
(3.2), it is then proper to enquire whether these equations are consistent 
with the original linearization in «. Returning, therefore, to the exact 
dynamical equations (2.3) and (2.4), the neglected inertia terms in the first 
of these equations are 

ace | 2x _ {x cos # — Oh @>-—< - yD") + 

r r 00 r’sin?@ o(7r,6) 


QDs { du Lap. 
et a ey cos 6— == SE OG > s 
r“sin’U | or r o0 j 





¢ 
t § 


r?sin?@ | ¢ 








which, on the assumption that a boundary layer exists at r = 1, reduce to 


2.02 E cos 0 ex 1 ( ous O3s ous oFds ) Z, ous =a | 
eal? p+ — ——| = = ~ , 


sin2@ * Or sin@\ c@ er? ~— er er?0@ sin? or or? 














Thus, since we may assume 
¥% = O(R-™), x = Of4), 
d/ér = O(R"?), 0/08 = O(1), 
the order of magnitude of the neglected terms is 
eal? R}?2, 
The order of magnitude of the inertia term retained in the first of equations 
(3.2), on the other hand, is 


Oy ; 
ean? Xcos@ or caQ?R?. 
or 


Hence the neglected terms are smaller by a factor e than those retained, and 
the presence of large velocity gradients does not affect the basic approxi- 
mation. ‘The argument may not, of course, apply near particular values of 
6, notably near the equator. Similar remarks apply to the other dynamical 
equation. 

The property of the equations (3.2) that is chiefly responsible for the 
success of this investigation is that they are both ordinary differential 
equations. This remarkable circumstance possibly embodies a point of 
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principle which might have been anticipated on general grounds. It is, of 
course, a consequence of the linearization that the entire secondary flow 
must be reversible. In particular, this reversibility must be true of the 
flow in the boundary layers. Now, the boundary layer flow in an axial 
plane presumably involves a general mass flux from the poles to the equator, 
or vice versa, depending on whether « is negative or positive (for the inner 
sphere). Hence, if one regards this flow as being determined by conditions 
just outside the boundary layer, together with an initial condition rep- 
resented by the fact that the flow is of some generalized ‘stagnation’ type 
at the poles (say), then the corresponding solution at the equator must be 
such that, when used as an initial condition for the reversed flow, it yields 
an identical flow pattern. ‘This would appear to impose a severe integral 
condition on the flow in the boundary layers, and by far the most attractive 
method of satisfying this condition is to suppose that the velocity profile at 
each value of @ is determined independently. ‘That such is the case is 
shown by (3.2), and the explanation given seems at least a possible one. 
One integration of (3.2) gives 


1a%% | 
2 cos 6(x — Xo) = Rap? | 
4 a 
1 ex (3 ) 
—2 cos Ab — py) = Ray’ | 
: 


where yy and yy are functions of 6 only, which are determined by conditions 
at the outer edge of the boundary layer. From (3.3) it follows that 
n4 


9 
we (b—w) = —4R* cos*6(ys— yp), (3.4) 
5 
so that a suitable measure of distance from the boundary is 
n = (r—1)(Rceos8)!2 (0 <@ < 4z). (3.5) 
The four independent solutions of (3.4) are then of the form 
exp{(+1+7)y}, 
and the bounded solution for % that has a double zero at » = 0 is 
Yb = po{l—e "(cosy +sin7)}. (3.6) 
The corresponding solution for % follows immediately from the first of 
equations (3.3), and is 
X = xXo—2(Reos 6)! 75, e~” cos n (3.7) 
Then, since x must also vanish at 7 = 0, we must have 
Xo = 2(Rcos 8)! hy (3.8) 
The formulae (3.5) to (3.8), which represent the explicit solution for the 
flow in a three-dimensional boundary layer, were first obtained by Ekman 
(1902). Actually, Ekman considered the flow near a rotating disc when 
the relative flow at a great distance from the disc is a uniform stream, but, 
to the order of the boundary layer approximation, this flow is clearly identical 
with the flow near the spherical boundaries of the present problem. Since 
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the boundary layer is very thin, each element of the surface of the sphere 
may be regarded as a large disc rotating with the normal component {2 cos @ 
of the angular velocity of the sphere, and the relative velocity outside the 
boundary layer is locally a uniform stream. 

The most striking feature of the solution is the oscillatory approach to 
conditions in the ‘inviscid’ core of the motion (the ‘Ekman spiral’). The 
same feature appeared, of course, in Squire’s and Stewartson’s solutions 
for the flow near a rotating disc. Another point of interest is the explicit 
formula (3.5) for the thickening of the boundary layer towards the equator. 
Very close to the equator this thickening becomes very rapid, and when 
}—@=O(R") the boundary layer approximation clearly breaks down 
altogether. 

It should also be noted at this stage that the relation (3.8), which is 
demanded by the viscous mechanics of the boundary layer on the inner 
sphere, implies a relation between the functions % and y throughout the 
whole of that part of the core of the motion which lies within the cylinder 
p = 1 (i.e. the cylinder that touches the inner sphere). ‘The reason for 
this is the form (3.1) that the solution must take in the absence of viscous 
forces. Actually, it is not strictly correct to apply the first of the results 
(3.1) without further explanation because (3.8) ensures that y = O(R-?) 
which vanishes in the limit of infinite Reynolds number. ‘The strictly 
‘inviscid’ solution is therefore 

~=0, x = (p). 
But, of course, when y% = O(R?) and y = O(1), it is still true that the 
viscous terms in the dynamical equation (2.15) are negligible in the core of 
the motion, so that the first of the results (3.1) is still valid. In this restricted 
sense, the motion in the core is ‘inviscid’, and we may write 


t= dlp), x = Xol(p)- (3.9) 


Since p = sin 6 on the inner sphere, it then follows that (sin @) and y,(sin 4) 
are the functions of @ only that appear in (3.3), so that the condition (3.8) 


becomes Xo(p) = 2R12(1—p?)'4y(p) (p < 1). (3.10) 


The analysis of the motion in the boundary layer on the outer sphere is 
substantially the same. Here, the integration of (3.2) gives 
ae 
2 cos 8{ x — x9(a sin @)} = Ro 
(3.11) 


& | 


lox | 
Ror’ | 
where ys, and x, have the same meanings as in (3.9), since «sin @ is the value 
of p on the outer sphere. Then, introducing the variable 
n’ = (a—1)(Reos 6)}2, (3.12) 
the relevant solution for & is 
ob = y(xsin9){1—e-"(cos y’+sin7’)}, (3.13) 


— 2 cos O{ib — ¥(« sin @)} = 
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and the corresponding solution for x is 

X = Xo(%sin @) + 2(R cos 6)! *yh,(x sin A)e-”’ cos 7’. (3.14) 
Moreover, the boundary conditions (2.12) require that y = a?sin*@ over 
the outer sphere 7’ = 0, so that (3.14) gives 


a” sin?@ — y,(xsin #) = 2(Rcos 6)! *(« sin 4). (3.15) 


Now, the result (3.15) provides a relation between the functions ys, and 
X» Which must be satisfied throughout the entire core of the motion. Thus, 


; 2 \ 14 ve 
p?— x0(p) = 2R!2( Ae 5) tp), (3.16) 


without restriction on the value of p. In the region p < 1, therefore, the 
two relations (3.10) and (3.16) uniquely determine the functions ys and yp, 
and the relevant formulae are readily found to be 


‘ pe j p” 14 ° 1 _ 
Ye) = peray(1- Sa) +0 (3.17) 
; p* 14 ) 1 
and Xo(p) = p*( l —p*)} 4. (1 ia f) +(1 —p*)i4 esd (3.18) 


The result (3.18) is particularly interesting inasmuch as it is an explicit 
example of a velocity distribution of order R° being determined by viscous 
mechanics in a case where Ris large. “The manner in which this takes place 
in the present problem is, of course, very simple. In the core of the motion, 
the streamlines in an axial plane are parallel to the axis of rotation, so that 
fluid leaving the boundary layer on one sphere round the circle p = p, must 
enter the boundary layer on the other sphere round a circle of the same 
radius. ‘This is one essential connection between the boundary layers on 
the two spheres. ‘The other essential connection arises from the fact that 
the azimuthal velocity in the core is also a function only of the radial distance 
from the axis of rotation, and must be such that the associated centri- 
fugal and coriolis forces are just sufficient to require the same mass flux 
from beth boundary layers at any particular value of p. 

It follows from (3.18) that the cylindrical surface p = 1 must be another 
singular surface of the motion. It is true that the formula (3.18) is not 
valid for values of p close to unity, but this cannot affect the deduction that 
velocity gradients become very large at such values of p. Thus, it seems 
likely that viscosity must again be taken into account near this surface. 

The problem of finding the distribution of % and xy») for values of p 
greater than unity is slightly different, because the result (3.10) is no longer 
relevant. It may be replaced, however, by the symmetry condition that 
y = 0 over the equatorial plane €= 0. Hence, either the equatorial plane 
is also a singular surface of the motion (so that ¢ at € = Vis not equal to 4,(p)), 
or else 


%o(p)=0 for all p > 1. (3.19) 
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As regards the former possibility, the appropriate boundary layer approxi- 
mations in the cylindrical polar equations (2.14) and (2.15) yield the equations 


Mee, Mera ae, 
O¢ Cc * og oC* 
so that 
a! ; 
aril —Yo(p)} = —4R*s—do()}- 


But there are clearly no solutions of this last equation that remain bounded 
for all values of { (both positive and negative), so that (3.19) appears to be 
the relevant deduction. It follows immediately from (3.16) that 
Xo(p) =p” for all p > 1, (3.20) 
which shows that the whole of the fluid outside the cylinder containing the 
inner sphere rotates as a rigid body with the same angular velocity as the 
outer sphere. ‘There is, of course, no boundary layer over the portion of 
the outer sphere that bounds this rigid body rotation. 
There now remains the problem of finding the motion in the free shear 
layer near the surface p= 1. ‘This problem is discussed in the following 


section. 


4. THE CYLINDRICAL SHEAR LAYER 
In the immediate neighbourhood of the surface p = 1, velocity gradients 
with respect to p must surely be very much greater than those with respect 
to ¢, so that it is again permissible to make a boundary layer type of approxi- 
mation in the governing equations. ‘Thus, the equations (2.14) and (2.15) 


may be written, approximately, 





CY tds ar 2 
Ws = ~2RY - eX (4.1) 


aL apt? aL Op? 


wx 


However, the approximate equations (4.1), unlike those governing the 
motion in the boundary layers over the spheres, are not always consistent 
with the linearization in e. For, if the exact equations of motion are 
written in cylindrical polar coordinates and a shear layer is supposed to 
exist at p = 1, it is a simple matter to show that the neglected terms on the 
left hand sides of the first and second of equations (4.1) are 

R( 2x4 + _ wa) ; 


AY 


of = =epap*af = aL Op 


opal = al ap 

respectively. Then, if 5 is the thickness of the layer, and if y is assumed to 
be O(1) throughout the layer (y is actually equal to the azimuthal velocity 
when p = 1), the orders of magnitude of these terms are 


2 | 
eR (1 + 4 and «RY , 
53 


and eR 


(x Ox dv *) 
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respectively. Hence, comparing the second of these estimates with the 
second of equations (4.1), the validity of the linearization is certainly depen- 
dent upon the condition that «/5 should besmall. ‘Thus, if ’—- 0 as R-> a, 
the equations (4.1) cannot refer to the limit R —> o for fixed small e«. How- 
ever, if, as seems likely, the orders of magnitude of #, x and 6, are determined 
only by R when « is small, the equations should apply to the limit* « — 0 for 
fixed large R. In this sense, the linearized equations (4.1) will be assumed 
to be valid. 

Now, it appears that the boundary layer type of approximation leading 
to (4.1) is the only feature that the shear layer has in common with conven- 
tional boundary layers. It does not seem to be possible, for instance, to 
have a shear layer of this kind in which dynamical variables are of the same 
order of magnitude as in ordinary boundary layers. For, if that were the 
case, we should have 


b= O(R2), x = O(1), 


together with the information that the thickness of the layer is O(R~?). 
But, with these orders of magnitude, the first of equations (4.1) reduces to 
any 
— = 0, 
cp 
and, quite apart from the dynamical unlikelihood of a layer with negligible 
acceleration in an axial plane, this equation does not possess solutions that 
are bounded as (p—1)R!? —> + ow. 

Again, it might be thought that the orders of magnitude must be such 
that viscous and inertia forces are comparable in both equations of motion. 
This would at least correspond in principle, if not in detail, to the mechanics 
of an ordinary three-dimensional boundary layer. Under these circum- 
stances, elimination of % from the equations (4.1) gives 





from which it follows immediately that the thickness of the layer must be 
O(R-'%), and hence, from (4.1), that 


/y = O(R™). (4.2) 


Now, it must surely be the case that y = O(1) throughout the layer, so that 
(4.2) gives y = O(R-™%), which implies a general mass flux by the velocity 
components in an axial plane of the same order of magnitude. But an 
important raison d’etre of the shear layer is to carry fluid from the boundary 
layer on one sphere to that on the other, and this mass flux is only O(R-?). 
The conclusion would have to be, therefore, that, to the order of the boundary 
layer approximation, the total mass flux parallel to the axis of rotation was 
zero, and this would imply that there was a substantial ‘return flow’ in the 
shear layer. I have not been able to show that such a situation is impossible, 
but the assumption of comparable viscous and inertia forces does seem to 


* A somewhat easier limiting process to realize in any practical arrangement. 
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be a most unnatural method of attempting to satisfy the conditions of the 
problem. 

The only acceptable alternatives appear to be those in which the orders 
of magnitude of #% and x in the layer are those suggested by the boundary 
conditions, namely 

p= OR"), — x = O(1), (4.3) 
and the thickness of the layer is sufficiently large to ensure that viscous 
forces are comparable with inertia forces in either only one, or possibly 
neither, of the dynamical equations (4.1). Taking first the case in which the 
mechanics is not entirely inviscid, it is a simple matter to see that the equation 
in which viscous forces are appreciable must be the second of equations (4.1) ; 
from which it follows that the thickness of the layer is O(R-'*). The 
governing equations may therefore be written 





Ox 0 _ 23" _ ay (4.4) 
og’ of p’?’ 
where p’ = (p—1)R"* (4.5) 
and y’ = R'?, (4.6) 


If this interpretation of the mechanics is correct, it appears that viscous 
layers in rotating systems can have an entirely different structure from 
ordinary boundary layers. ‘This structure is a highly anisotropic one in 
which viscous forces oppose only the rotary motion, so that the motion in 
planes containing the axis of rotation remains ‘inviscid’. In the special 
problem under discussion, the inviscid nature of the flow in an axial plane 
provides a neat explanation of how it is that fluid travels from one sphere 
to the other within a layer which remains thin. ‘The essential point is that 
the layer is sufficiently thick, and the velocity sufficiently small, to make 
viscous forces negligible. 

The general solution of (4.4) is easily seen to be 


p= -3LF'(p')+ Gp’), | 
where F and G are arbitrary functions of p’, and F’”’ is the second derivative 
of F. The form of these arbitrary functions is presumably determined by 
conditions at the points ¢ = 0 and ¢ = («?—1)"* where the shear layer joins 
the boundary layer on the spheres. Unfortunately, the solution for the flow 
in the neighbourhood of these points lies outside the scope of the boundary 
layer approximations, and would require a much more comprehensive 
analysis of the motion than I have been able to make. However, it can at 
least be shown that the structure of the solution (4.7) is consistent with the 
boundary conditions as p) ++. Asp’ ->+0o, the motion must tend to 
the rigid body rotation represented by (3.19) and (3.20), so that we must 
have 

F(oo) = 1, G(o«) = 0. 
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As p’ -»—o, on the other hand, the boundary conditions are found by 
writing p = 1 in the formulae (3.17) and (3.18). In this way we get 
1 1\-1/4 
F(—o) = 0, G(- 2) =3(1- 5) 
It may then be noted that F’’(p’) + 0 as p’ -» +, so that the dependence 
of the solution (4.7) on ¢ disappears as the inviscid regions are approached. 
‘The remaining possibility is that the mechanics of the shear layer is 
entirely inviscid. ‘This clearly corresponds to the case in which the thick- 
ness of the layer is of greater order of magnitude than R-!*. In sucha case, 
both the thickness of the layer and the velocity distribution within it (now 
of the simple form (3.9)) would presumably be determined by conditions 
near the joins of the shear layer with the boundary layers over the spheres. 
It would seem that the only satisfactory way of deciding which, if 
either, of the preceding two mechanical structures is correct is to determine 
the asymptotic form, tor large Reynolds numbers, of the exact solution of 
the linearized equations (2.9) and (2.10). 


I am indebted to Mr W. W. Wood for several helpful comments on an 
earlier draft of this paper. 
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The speed of drifting bodies in a stream 


By J. R. D. FRANCIS 
Imperial College, London 


(Received 12 May 1956) 


SUMMARY 
The speed of free floating bodies on the surface of a water 
stream in a sloping channel has been found to be sensibly the same 
as the mean speed of the layer of water in which the bodies are 
floating, contrary to some recorded opinions. 


It is sometimes believed that a ship drifting in a river without power 
or sail travels faster than the mean speed of the layer of water in which 
it is floating; it is said that the relative speed of the water past the rudder 
enables the ship to be steered. Prandtl (1952, p. 179) for example says 
‘“‘ The boat, in fact, hurries ahead of the water, and gets enough way on her 
to be steered ”’. 

The basis of the belief is that the weight force W of the ship is directed 
vertically downwards, but that the hydrostatic buoyancy force on the ship 
is directed at right angles to planes of equal pressure, that is, nornial to the 
free surface. Since the surface in a river in uniform flow is inclined by 
friction at an angle 7 to the horizontal, it follows that a component force Wi 
(¢ = sinz for usual slopes) acts on the ship in a direction down the surface 
slope. ‘This component force is balanced in steady motion by the hydro- 
dynamic drag force due to the relative motion between ship and water. 
Alternatively, it is believed that a ship travels faster than the equivalent 
volume of water because there is an exchange of momentum by turbulent 
motions across the boundaries of the water. Since the mean speed of the 
stream decreases with depth, the exchange of momentum causes a resistive 
force to be exerted on the water volume. ‘The solid boundary of the ship 
prevents this exchange, so that the only forces on the ship restraining 
the down-gradient motion are shear stresses caused by the motion of the 
ship relative to its surroundings, it being tacitly assumed that the flow and 
pressures in the surroundings are not affected by the replacement of the 
water volume by the ship. If, in fact, the solid boundary modifies the flow 
in the neighbourhood, then the forces on it may be quite unlike those at 
the boundaries of the same volume of water. ‘The buoyancy force on the 
ship may not then be normal to the free surface, so that the down-hill force 
can no longer be so easily found. 

It is difficult to show if ships do indeed travel faster than their 
surroundings; for the effect must be small in slow flowing rivers (with 
small 7), and it is clearly imprudent to allow a ship to drift unrestrainedly 
in a fast and turbulent river (which has a larger 7). It is more practical 
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to test the existence of this relative velocity on a cylinder floating with its 
axis vertical in a stream. Such floats are sometimes used by engineers to 
determine the speed of a river, and it may be desirable to ascertain if there 
is any systematic error in using them for this purpose. 

The force W on such a cylinder of diameter d and length / when floating 
in a fluid of density pis W = }md?lpg. If the cylinder is moving at a velocity 
u,, relative to its surroundings, then the drag force F on it is approximately 
that which the same cylinder would experience if it was travelling at the 
same speed relative to still water. Onthis approximation, F = C,dl(}pu,;,); 
and, equating F to the downstream component of W, we have 

C,dl(ipuz,) = 47d? lpgi, 


rel 
where C, is the coefficient of drag of the cylinder. It is convenient to 
substitute the frictional properties of the river channel for the variable 2, 
and one empirical formula used by engineers is 
u = Cy/(mi), 

where C is a coefficient (Chezy’s coefficient) predominantly determined 
by the roughness of the channel, wu is the mean velocity, and m is the hydraulic 
mean depth, i.e. the cross-sectional area divided by the wetted perimeter 
of the channel (equal to the depth if the channel is wide). Substituting 
for 7 in the previous equation, we obtain 


Ur I md \ 
uC 2mC 4 


It will be seen from this equation that a high value of u,,, should be 
obtained in rough-walled channels (having a low value of C) and with 
cylinders having a large d/m ratio. Engineer’s floats may be cylinders of 
about 3 in. diameter, and might well be used in a typical river of depth 
m = 10ft., C = 100 ft.!2 sec-1. The equation shows that the theoretical 
value of u,,,/u is about 0-01; it is probable that this small increase of velocity 
(if it exists at all) is hidden by experimental error and by the large-scale 
turbulence of the stream. 

The point may be better investigated in a laboratory channel which 
can be made to have a great roughness, and experiments have been so made 
in the Civil Engineering Laboratories of Imperial College. Systematic 
roughness elements, of a type proposed by Denil and tested by White & 
Nemenyi (1942), were arranged on the bottom of a glass sided channel 
11 mlong and 30 cm wide. ‘These roughnesses are zigzag walls 1-6 cm high 
running across the channel. From alternate 90° angles, longitudinal walls 
of the same height connect one wall with the next, 7-5 cm upstream or 
downstream. ‘The bed of the channel may be tilted so that it is parallel to the 
surface of the water. Four floats were made, all circular cylinders ballasted 
to float to a draft of 2-5cm with their axes vertical. A quantity of 
0-00652 m* sec"! was made to flow down the channel; and it was found that 
if the bed was inclined at 7 = 1/271, the depth was constant at 7-5 cm above 
the top of the roughness. Under these conditions, C = 21 m!? sec”? for 
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u = 0-286 m/sec. 
in table }. 


In this table the coefficients of drag have been taken 
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The theoretical value of u,,,4 for each cylinder is given 


as those 


applicable to cylinders whose length-diameter ratio is double that of the 
cylinders tested, so as to allow for the three-dimensional flow round the 


lower end (Goldstein 1938, p. 439). 

















Float diameter (cm) : | 6°3 1-9 0-7 0-2 
Calculations | 
| Length diameter | O-+ 1-3 3-6 12°5 
| Drag coefficient | 0-6 0-7 0-8 0-95 
| Upey iu | 0-270 0-137 0-078 0-038 | 
| 467 (cm/sec) Pee 5-9 pi. os ee 
pees | — | 
Observations | 
Number of successful drifts | 20 21 22 16 
Mean time to traverse 2°88 m (sec) | 6°83 6-90 6°80 6-89 | 
Mean deviation (sec) | + 0-24 + 0-31 + (0-25 + 0-24 
Mean speed (m/sec) | 0-421 0-417 0-423 0-418 
| | 








Table 1. Calculated and measured velocities of cylinders 2-5 cm long drifting in a 
stream of mean velocity 0-286 m/sec. 

The observed speeds of the floats are also shown in table 1. They 
were timed to the nearest 1/10 sec with a stopwatch calibrated in 1/100 sec, 
over a 2:88 m length of the channel, starting 5 m downstream of the channel 
inlet and finishing 3m upstream of the outlet. The floats were released 
1m upstream of the starting line in mid-stream; and, if the float later 
drifted so that its axis came outside of the middle third of the channel, 
its time was disregarded. Only about 1 float in 6 drifted successfully within 
the middle third. 

It will be seen that the calculation predicts that the largest cylinder 
should travel some 6:6 cm/sec faster than the smallest, but that the observa- 
tions do not disclose any significant difference in speed for the large range 
of d m used. 

It is, however, possible that a// the cylinders drifted at the same speed 
relative to the upper 2:5cm of water, in which they all were immersed. 
Two tests were made to explore this possibility. In the first, drops of dye 
were put into the stream just ahead and around the largest cylinder. Onno 
occasion did the cylinder appear to overtake the dye and to float into clear 
water. If the relative velocity exists, then it should have overtaken at 
7-7 cm/sec, a speed easily observed. 

In the second test, carried out concurrently with the timing of the 
floats over the 2:88 m test distance, the mean speed of the upper 2:5 cm of 
water was measured by a current meter. ‘This was a paddlewheel, 30 cm 
radius, rotating on a horizontal axis, with light aluminium blades immersed 
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to a depth of 25cm. ‘The part of the wheel in the air was protected from 
draughts. ‘The meter was calibrated before and after the test by towing 
it through still water. ‘The mean speed of the water was found five times 
during the timing of the floats both upstream and downstream of the test 
length. ‘The mean water speed so obtained was 0-420 m/sec; this was 
close to the mean speed of all the floats (0-420 m sec + 0-017 m, sec mean 
deviation of 79 observations). 

It therefore appears that the floats were in fact travelling, as nearly as 
could be measured, at the same speed as the water; and that the ‘hurrying 
ahead’ of the float was either not present at all, or was much smaller than 
predicted. Perhaps the presence of shear in a stream rearranges the 
hydrostatic pressure on the float, so that the resultant upthrust exactly 
balances the weight force; in this respect, a drifting body therefore appears 
to affect the flow in the surrounding fluid. This is fortunate, for many 
experiments in fluid mechanics are carried out using particles of near 
neutral buoyancy as tracers. If the gradient effect had been real, then 
systematic errors would occur if pressure gradients occurred. In 
meteorology, pilot balloons would not travel with the wind; since the 
atmospheric pressure gradient is at right angles to the wind, a sideways 
component would be given to the balloon, giving a false direction. ‘The 
above experiments should not be taken to apply to a flow which has a surface 
gradient but no shear, in which case a floating object might possibly move 
relative to the water down the gradient. 

If, as was mentioned at the beginning, ships can indeed be steered while 
drifting down a river, then a possible explanation independent of the 
‘hurrying ahead ’ effect may be that the ship’s draft is nearly the same as 
the depth of the river. ‘The lower part of the rudder is therefore in water 
which is considerably retarded by the boundary layer of the river bed. The 
ship travels at the mean speed of the upper layers, which is faster, so that 
there is relative motion over the rudder and the ship is steered. ‘The effect 
will be enhanced if the ship is of small draft, with a rudder projecting well 
below its bottom. ‘Tests carried out in the laboratory channel show in 
fact that a model boat can be steered when drifting if it has such a deep 
rudder. It cannot be steered if the same rudder is arranged so that it does 
not project below the boat. 


I am indebted to Professor C. M. White, who pointed out that Prandtl’s 


statement was untested. 
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as 
ing SUMMARY 
we In this paper the spatial variations and spectral structure of 
the steady-state turbulent shear flow in channels are investigated 
hy without the introduction of empirical parameters. ‘This is made 
one possible by the assumption that the non-linear momentum trans- 
ne port has only stabilizing effects on the mean field of flow. ‘Two 
ear constraints on the possible momentum transport are drawn from 
en this assumption: first, that the mean flow will be statistically 
In stable if an Orr-Sommerfeld type equation is satisfied by fluctua- 
he tions of the mean; second, that the smallest scale of motion that 
_ can be present in the spectrum of the momentum transport is the 
- scale of the marginally stable fluctuations of the mean. Within 
ce these two constraints, and for a given mass transport, an upper 
ne limit is sought for the rate of dissipation of potential energy into 

heat. Solutions of the stability equation depend upon the shape 
le of the mean velocity profile. In turn, the mean velocity profile 
™ depends upon the spatial spectrum of the momentum transport. 
ae A variational technique is used to determine that momentum trans- 
“ port spectrum which is both marginally stable and produces a 
” maximum dissipation rate. The resulting spectrum determines 
t the velocity profile and its dependence on the boundary conditions. 
t Past experimental work has disclosed laminar, ‘transitional’, 
1] logarithmic and parabolic regions of the velocity profile. Several 


experimental laws and their accompanying constants relate the 
: extent of these regions to the boundary conditions. ‘The theore- 
tical profile contains each feature and law that is observed. First 
approximations to the constants are found, and give, in particular, 
a value for the logarithmic slope (von Karman’s constant) which 
is within the experimental error. However, the theoretical 
boundary constant is smaller than the observed value. ‘Turbulent 
channel flow seems to achieve the extreme state found here, but a 
more decisive quantitative comparison of theory and experiment 
requires improvement in the solutions of the classical laminar 
stability problem. 
INTRODUCTION 
The significant difficulty encountered in the theoretical study of fluid 
turbulence is the non-linearity of the equations of motion. The classical 
inquiries (e.g. Goldstein 1938), which are concerned with gross spatial 
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wa 


features of the flow, avoid dealing with non-linearity by deducing 
macroscopic laws from the observations. ‘The more recent statistical 
inquiries (e.g. Batchelor 1953), which are mainly restricted to homogeneous 
fields of motion, deal with non-linearity by various hypotheses concerning 
the transfer of energy from one scale of motion to another. ‘The emphasis 
in these latter studies is on the local random aspects of turbulence and 
on those regions of the spectrum of turbulence far removed from external 
energy sources. Hence the spatial variations and spectral structure of 
momentum and heat transport, which are intimately coupled to the energy 
sources, receive little theoretical attention. 

In this paper the emphasis is on the statistical organization of turbulent 
fields rather than on their statistical randomness. Hence the momentum 
and heat transport of a steady-state turbulent field and the dependence of 
these transports on boundary conditions are of primary concern. However, 
this will not be a detailed mechanistic study. Indeed, there is little hope 
that exact time-dependent solutions of the non-linear equations describing 
turbulence can be found. Here, as in past work, an assumption specifying 
the role of the non-linear processes will be made. This particular assump- 
tion is a negative one, specifying what the non-linear processes can not do. 
In this way the class of possible turbulent fields is restricted. The deter- 
mination of an extreme member of this class reduces the mathematical 
problem to tractable linear form without the introduction of empirical 
parameters. ‘This absence of empirical parameters permits a rigorous 
appraisal of the validity of the assumption. Indeed, if the quantitative 
solutions to the problem are to compare favourably with the observations, 
the assumption must embody all the significant restraints imposed by the 
equations of motion. 

The particular physical situation to be discussed is the fully turbulent 
flow in channels. ‘This flow has been studied experimentally for many 
years. Hence it is an excellent testing ground for theoretical models. 
The model evolved here has its origin in conventional laminar stability 
theory. ‘Therefore a brief restatement of the conclusions of this theory 
is of value. 

The stability of a time-independent state of a fluid field can be established 
by study of the growth or decay of superimposed infinitesimal disturbances. 
If any disturbance releases a larger amount of potential energy from the 
pressure field than it dissipates through the action of viscosity, it will grow. 
The mathematical formulation of this kind of problem usually leads to a 
linear characteristic value equation. ‘The external conditions which delimit 
the range of stability of the initial time-independent state of the fluid are 
determined by the minimum characteristic value. In simple shearing 
flow a disturbance can abstract energy from the pressure field only through 
the agency of viscous stress. Hence, in this case, viscosity plays the double 
role of permitting the growth of disturbances and also dissipating them. 

When a disturbance achieves finite amplitude, its non-linear interaction 
with the initial motion alters both this motion and the disturbance. 
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First-order stability theory tells one nothing of this alteration, while inclusion 
of the non-linear interaction in the problem leads to such intractable mathe- 
matics that even second-order effects have seldom been investigated. 
However, the physics of the non-linear effects is fairly clear. ‘The disturb- 
ance must grow until its alteration of the initial motion reduces the rate of 
release of potential energy to the rate of dissipation of this energy. If the 
disturbance is an aperiodic turbulent motion rather than another steady 
laminar motion, the equalization of rates of energy release and dissipation 
must still occur on the average. In either the laminar or turbulent case, 
alteration of the initial motion is accomplished by the transport of heat or 
momentum, whichever was responsible for the instability. Hence the most 
important aspects of the non-linear interaction of disturbance and initial 
motion are those advective terms which readjust the energy sources of the 
disturbance. 

In the new ‘stability’ that is achieved, the finite-amplitude disturbance 
can be looked upon as the stabilizing sink of energy for the destabilizing 
energy source in the pressure field. It will be assumed that this is so on 
the average even in the fully developed steady-state turbulence, i.e. that 
the mean non-linear momentum transport terms in the equations of motion 
are entirely stabilizing, and that the only energy source for those disturbances 
which enter into the mean momentum transport is the pressure field of the 
mean motion. 

In the first section, a limiting ‘statistical stability’ condition on the 
spatial structure of the mean momentum transport is found as a consequence 
of this assumption. ‘The momentum transport is then resolved into a 
spatial spectrum whose smallest scale of motion is the scale of the smallest 
statistically unstable motion. ‘The determination of the particular spectrum 
satisfying these stability constraints and the boundary conditions, and also 
releasing the largest amount of potential energy, leads to an unusual formal 
problem. In this problem the non-constant cofficients of a linear character- 
istic value equation are varied under constraints to determine a minimum 
least characteristic value, also the corresponding minimum least character- 
istic function and the optimum form of the non-constant coefficients. 

The second section outlines a variational approach designed to find these 
optimal functions and constants; i.e. to find the momentum transport 
spectrum and mean velocity profile which lead to a marginally stable mean 
field of maximum dissipation rate. Considerable use is made of the 
knowledge acquired in the past thirty years about solutions of the laminar 
stability equation for shearing flow. A first approximation to the spectrum 
is found, and a comparison is made between this theoretical extreme and 
the observations. 

The concluding section describes other more detailed tests of this work 
by comparison with recent channel flow data, e.g. the observation of a 
smallest scale of motion. An attempt to rationalize the basic assumption 
of the theory raises questions concerning the uniqueness of solutions to 
the equations of motion, while the apparent existence of extreme states in 
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turbulent fields suggests an explicit statistical mechanical origin. Finally, 
the application of this approach to quite different turbulent phenomena 
is outlined. 


1. FORMULATION OF A TRACTABLE PROBLEM 

This is a study of the steady-state fully evolved three-dimensional 
turbulence observed between fixed parallel surfaces. Here ‘fully evolved’ 
will mean both that the Reynolds number (based on the channel half-width 
and the average velocity) is far above its value for marginal stability, and that 
there are no variations of moments of the fluctuating field in the downstream 
direction. In figure 1 the coordinate frame is chosen so that the x direction 
is the direction of the mean velocity U, which in turn is some unknown 
function of the cross-stream direction z. ‘The ‘mean’ will be defined 
here as an average in the y direction (perpendicular to the plane of figure 1). 
It is assumed that the effect of the distant y boundaries on real flows does 
not influence the nature of the motion far from those boundaries. Hence 
the remainder of this work deals with parallel surfaces of infinite extent. 
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Figure 1. Coordinates chosen to describe flow between parallel surfaces 
It is also assumed in general that compressibility, heating effects and 
non-linear viscous effects are unimportant in the experimental situation. 
Therefore the equations of motion and continuity are 
Dv l 
Pa - VP+V*v and V.v = 0, (1.1) 
Dt p 
where D/ Dt is the substantial derivative, P the pressure, p the fixed density, 
v the constant kinematic viscosity, and v the velocity vector. 
A first concern will be the determination of statistical stability constraints 
imposed on the steady-state flow by these equations. 
As implied by the name, statistical stability means the constancy in 
time of some space (or ensemble) average of the field of motion. Here an 
investigation will be made of the time dependence of the quantity 


1 ft 
¥(x,2,t) = E | vdy | =v-v’. (1.2) 


u 


Future use of the bar will indicate the mean with respect to y, as in (1.2). 


Expanding (1.1) with the aid of (1.2), one may write 
ev. ov ae ae : 1 ie seis 
> +a + v.Vviv.Vv'+v .Vviv’.Vw = ——VP+ 1V2v+ Vv’, 
at oO p 


and V.v+V.v' = 0, (1.3) 
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The mean value “ 1.3) is then 


Hak 


+v.Vvtv'.Vv' = _ $y Fats 
p 


Y 
ow 


and V,.v = 0, (1.4) 
where V, indicates the two-dimensional (x, z) del operator. Subtracting 
(1.4) from (1.3), one obtains 


4. Vv'+v'.Vv)+(v’.Vv'-—v v’.Vv’) = = — “(TP V, P) + vV2v’, 


and V.v’ = 0. (1.5) 
Now if the mean flow is to be stable and independent of x, i.e. if 
v = U(s)i, then if follows from (1.4) that 


A277 > > 
| -* se i+|-' = ji+[ 5 = |k- = -[ wv ji-| ow k, (1.6) 
os” | pox pos 


where i and k are unit vectors in the x and z directions respectively, and 
u’ and w’ are the components of v’ in the x and z directions respectively. 
Equation (1.6) is, of course, the usual equation for the mean flow. _ Inte- 
grating the k equation, one finds that Pp+w’? = P,/p, where Py is the 
pressure at either fixed surface; and since @w’*/éx = 0 (no variation of 
mean flow moments in the x direction), then ¢P/éx = dP, ex. Integrating 
the i equation, one finds that, since there must be no torque on the fluid 





as a whole, 


vB+w'u' = 7 =T%)-, (1.7) 
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where 8 = —cU/ez, and 7, = -(- 0 P,/cx }&%> is the stress per unit mass at 
pP 


the fixed surface. Equation (1.7) specifies the spatial relation of the mean 
gradient and the mean momentum transport. ‘Thus it is basic to the 
following work. 

However, the stability of any particular mean flow is not assured unless 
an arbitrary mean disturbance disappears in a finite time. If we take 

v = U(z)i+ev,(x, 2, t)+ ev,(x My Bb) + sy 

and P = P(x, 2) +¢€P,(x, 2, t) + P,(x, 2, t)+... (1.8) 
where ¢ is a constant parameter, then V,, P,, V2,,P., etc. must all decay in 
time for an arbitrary choice of «. From (1.4), the zeroth-order terms in 
e give (1.6). The first-order € terms in (1.4) are 


Il 


ov OV eer a ia 

=< + St (SS Ji oH + SP, = —[Av’.Vv’] 
ct Ox Oo ii |) 

and V,.¥, = 0, (1.9) 

where [Av’. Vv’] signifies any first-order change that occurs in the mean 

non-linear stress due to the complex interaction in (1.5) of v, and P, with 
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the field of finite fluctuations. If arbitrary v, and P, are to decay in time, 
equation (1.9) places restrictions on the form and amplitude of U(z), and 
hence, through (1.7), on the form and amplitude of w’u’. 

At this point in the analysis, the mechanistic assumption that v’. Vv’ 
plays only a stabilizing role is introduced. The conditions placed on 
stable forms of U(z) and w’u’ will be at least as severe as the actual case if 
the assumed stabilizing effect of [Av’.Vv’] is neglected. ‘Then, with 
[Av’. Yv’] set equal to zero, equation (1.9) assumes the familiar linear form 
of the equation for the disturbance of a time-independent shearing flow. 
Here, however, the problem is to find the entire class of U, for given 
boundary conditions, which lead to marginally stable v and P. 

‘Taking the curl of (1.9) to remove the pressure term, and assuming that 


¥,.k = (zs) expi(x/z)[v—(cU,,)f], (1.10) 


one may write 


; . ‘i 
a A ee ee a4), (1.11) 
aR \e(s 3)" O(3/%)° 
where %/3, is the downstream wave number of the disturbance of the mean, 
cU,, is its phase velocity, U,,, isthe s-averaged mean velocity, and R = 3)U,,,/v 
is the Reynolds number. Equation (1.11) is called the Orr-Sommerfeld 
equation (Orr 1906-7, Sommerfeld 1908), and has been studied for many 
years to determine initial instability criterion for a given laminar U. Its 
appearance in a problem of fully turbulent flow is due in part to the present 
special assumption and in part to the vanishing of explicit interaction 
terms between v and v’ in (1.4). 
In his study of equation (1.11), Heisenberg (1924) has shown that for 
a given R some smallest scale of motion (largest ~) is unstable. (This is 
also established in equation (2.41) of this paper.) In addition if w’w’ is 
entirely stabilizing, it cannot produce motions requiring greater destabilizing 
force than that available in the mean flow. Hence a second consequence 
of the assumption is that the spectrum of w’w’ cannot contain a scale of 
motion smaller than twice the scale of the smallest marginally stable motion 
according to (1.11). (We specify ‘twice’ since w’u’ is a product of the 
Huctuation velocities.) as 
This restriction on the spectrum of u’w’ is, through (1.7), a restriction 
on the possible spatial structure of U. Hence the class of possible U is 
doubly restricted by the marginal solutions of the stability equation (1.11). 
Within this class, for a given U’,, one function leads to the release of more 
potential energy per unit time than any other. This function can be 
thought of as the ‘most stable’ in the sense that any alteration of it will 
increase the potential energy available to disturbances. Ignorance of the 
detailed mechanisms of momentum transport permits one to hope that this 
extreme state is approache the observed turbulence. Hence in the 
remainder of this sectionat: — sle formal problem for the extreme is posed. 
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In the steady state, the rate of release of potential energy equals the rate 
at which mechanical work is done by the shearing stresses, and this equals 
the rate at which energy is dissipated into heat by the action of viscosity. 


Hence 
ov\2 ov\2 dv\2 
(Br),, = | (x) + (=) aL (=) | , (1.12) 


where (év/ex)? = (év/cx) .(ev/ex), etc., the subscript m again indicates an 
average value in the s direction, and 


(Br), = T9(Bs So)m = To Bo (1.13) 





from (1.7). Also 


. : ~ | Sake ae Oe 
“Oo « Zz, “0 a, 


aSy . So “39 


Bo ia Is. | B 





Hence the total rate of dissipation per unit mass, 7) U,,,, will be a maximum 
for a given U,, when 7, is a maximum. In a non-dimensional form, a 
maximum for (%9/v)(7)/U,,,) is sought for a given (z9/v)U,, = R. However, 
it is convenient to invert this problem and seek the minimum Reynolds 
number associated with a marginally stable solution to (1.11), holding 
(z»/v)(7>/U,,) constant. This latter quantity is proportional to the ratio 
of the actual dissipation rate and the dissipation rate which would result 
from a laminar parabolic flow with the same U,,. Indeed, it is readily 
shown that 

(To )aminas = 3 vp. (1.15) 


To 1 /2 TT B(z) 
Shani” 3(5) (UL) = i wr 


will be held constant during the minimization of R. ‘The other variables 
in (1.11) for which optimum values will be sought are «, c and the mean 
velocity U. ‘There are, however, certain geometric and boundary con- 
straints on U which first must be determined. 

One notes that since u =v = w= 0) and ow/dz = 0 (from V.v = 0) at 
the boundary surfaces z) and — 2», then 


Hence the ratio 


at ts, and —2,. Hence, from (1.7), we have 

















1 s —— al ' T; a T, 
B=-(7——-u'w')= —-~—; thus A(z) =—, B(-%) = —-3 
Vv Zo os Vv Vv 
op T 1 ou'w’ ol op op 7 
I 0 7 | i 0 
Re en nee oe ae thus AS (Zp) = es Zo) a Spat 
C« Vso V OS COs” Oz Cs V&o 
e228 1 32u’w’ au 28 a28 = 
eg Be Rs eg ae Ot ae thus 22 (20) om 5,3 (— 20) an; (1.17) 


and also U = 0 at 3 and — 2,. 
Yet another constraint on U arises from certain general conclusions 
regarding the stability of solutions of (1.11). It has been found (Heisenberg 
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1924, Lin 1945-46) that a necessary, but not sufficient, requirement for 
stable solutions at large Reynolds numbers is that the curvature of the 
mean flow be always of the same sign. A type of instability called 
‘inflexional’ can result if this is not the case. Hence one further requires of 
U that ney 

— >> = Const. x d?, (1.18) 


os 





where d is a real quantity. It is convenient to define the dimensionless 
quantity 


l 
U = Tw (1.19) 
and choose the constant in (1.18) so that 
o2// 12 Ty }2 Kd (1 20) 
— = —- —-—— @ —Kd*, 2 
d¢* 7° (To)1, ) 


where ¢ = 37.(1—2/z)) (see figure 1). In terms of ¢, the boundary 
conditions on d?, for equations (1.17), may be expressed as follows: 


ts ' (Ss a @: 
(4) 0b} or 7 - 


(b) (H), =~ {MM = 1; 
i; 1.21 
(c) -| @dp=1; ( ) 
2/7 "+20 
(d) =| pd? dd l, or sé dz = 0. 


‘The definitional constraint 7! | //dé = 1 may be written 
~ O 
oe 4 
5 | ¢°d? dp = 1- —; (1.22) 


or? 
/ 0 a7 


Now (d) of (1.21) is just the requirement that d? be symmetric ; accordingly, 
equation (1.22) may be written 


ne ¢?\ 1, 
oe (s- Jar a, (1.23) 


— 
~ 0 a 


From (1.7) we have 





uw’ 2 +? 


| (d?—1) dé. (1.24) 


Hence the spatial spectrum of the momentum transport is determined 
by d*. We have also 


= Ki ind—d | d2 db+ [ da? dh. (1.25) 


The final constraint on U is due to the requirement that the field of 
transporting motions, u’w’, has some smallest scale, eventually to be 
determined by the stability equation. This constraint is introduced by 
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the requirement that d has some smallest scale mo, that is 
d= Y/Y, 4,@. (1.26) 


where ¢,,(¢) is some real orthogonal set of functions satisfying the appro- 
priate boundary and symmetry conditions, and the Y,,, are arbitrary reai 
amplitudes. 

The formal linear problem may now be stated: First, for any U satisfying 
—(U/dz?) > 0, find from the Orr-Sommerfeld equation (1.11) the 
smallest scale of motion and the relation between R, x and c for marginal 
stability. Second, expressing U in terms of the d from (1.18) and taking 
some convenient set y,,, vary Y,,, subject to the constraints (1.21) and (1.23), 
to find the minimum R for constant K. This analysis is undertaken in the 
next section. 

2. A VARIATIONAL SOLUTION TO THE STABILITY PROBLEM 

The stability equation (1.11) has been studied since 1908. A recent 
interesting inquiry (Lin 1945-46) contains, as a partial list, seventy 
references. A first purpose in all these papers was to establish the 
Reynolds number for the initial instability of laminar flow. Only an out- 
line of this theory will be given here. 

For a given velocity profile U, the Orr-Sommerfeld equation determines 
a relation c = c(x, R) for each solution ®. For R and « real, and c complex, 
this single relation can be separated into a real and imaginary part. For 
marginal stability the imaginary part of cis zero and this determines R = R(), 
the curve of marginal stability. For each point on this curve the real part 
of the relation c = c(x, R) determines the corresponding value of c. No 
complete solution has yet been found. However, asymptotic solutions of 
equation (1.11) for large values of xR have proved adequate in the study of 
initial laminar instability (and should be much better for the turbulence 
problem since «R is orders of magnitude larger). Recent quantitative 
work rests on the relations* : 

(®, / =) = $, F(n), (2.1) 
/ op ¢=0 
where ®, is the highly oscillatory part of the asymptotic ®, and @, is that 
value of ¢ at which U =c; 











7 = [C1 +A) —y)'9] GPK) 8c(2R), (2.2) 
where 
£i*, A 1 (4/2 (¢ ., d - 
yas | ‘ d* dd, ris | : (- | é d? ds) ?; (2.3) 
and 
fo def eutde He) (3 02)*2] ; 
F(n) =* = 1 —- ——-: 
-7 | é 2déH\)) [2 (ie)? 2] J (n) 


* Lin’s equations (6.28), (6.29), and (12.32) correspond to (2.1), (2.2), and (2.4) in 


the notation used here. 
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The function F() was first numerically calculated by Tietyens. Lin has 
improved and extended this calculation, of which relevant results appear 
in table 1 below. 





| » | 26 | 28 


i) 


3-4 3-6 3°8 


Ww 


( 91 | 0-22520 | 0-27193 | 0-30705 | 0-32130 
0- sis 72 | 0-58082 | 0-56401 | 0-51074 | 0-43560 
0:-46456 | 0-41947 | 0-36110 | 0-28802 | 0-20352 


| 
| 
| 


: Eee 

Fy) | 0-07203 | 0-12220 
Jn) | 0-31558 | 0-46043 
Fn) | 052773 | 049952 

















Table 1 


A first relation imposed by the boundary conditions on the several 
parts of the asymptotic ® is written 


pee Moz og? 
J () = - Ia?K(1 -20| Taa7ar | 


= 2¢,d?[(1+A)(1—y)]} = v(c), (2.4) 
where the primary contribution to z(c) is due to the pole of equation (1.11) 
at U = c, and where A, as will be justified @ posteriori, is assumed to be very 
small. ‘The value of ¢, is determined from (2.4) at the first interception 
of v(c) and J;(y). Equations (2.1), (2.2) and (2.4) can then be used to 
determine «R as a function of c; hence they determine (xR)min and the 
corresponding c. An important characteristic of the R(x) curve, used 
since Heisenberg first developed this approach, is that the value of R at 
(zR)min is very close to the minimum value of R. This property of R(x) 
will also be used here to find Rmin. The real task of recent work has been 
the evaluation of x from the complicated second relation imposed by the 
boundary conditions on the asymptotic ®. This will be avoided here 
through the use of an integral condition for stability discussed later in 
this section. 

‘Turning now to the turbulence problem, one notes that the scale of 
motion in @ is determined by the imaginary part of the derivative of ®. 
The largest value of this derivative will occur in the sharp shear zone at the 
boundary. Hence 


Then, from (2.1), one may write 


l 
eS. = ae ke 2.6 
°F) ) 
where €, = m)¢. Equation (2.4) is now written 
(d?/no), = 3[((1 +A) —y) PA.) Fi). (2.7) 


To determine (xR)min, Heisenberg first observed that for c(2R)/ec=0, one 
has from (2.2) and (2.4) 
J (y) = cv'(c), (2.8) 


and (aR) min a Gd mK)? (“) 1 (1 4 As(1 ee y), (2.9) 





| (1 2s) | , 
~-G3)L.. 2s 
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where 7 is that » which satisfies (2.8), and where the small spatial variation 
of (1+A)3(1—y) has been neglected. Hence, if the Y,, in the expansion 
of d?, equation (1.26), are constrained to satisfy (2.7) and (2.8), the resulting 
profile will be both marginally stable and continuously at (2R)min, no matter 
what further variations of Y, are made. 

Now a convenient orthonormal set ¢,(¢) which automatically satisfies 


(a) and (d) of (1.21) is 


ws,(6) = 1/2 cos 2nd, Wo(?) = 1, (2.10) 
where n is any integer. Hence, to satisfy (b) and (c) of (1.21), we must have 
¥ Y2 =1, (2.11) 


and 


Equation (1.23) may be written 


a. Me 2.43 
K sy 2 Bun ; ; . ( ’ ) 


where 
db? 


B = | (¢ - = [cos 2(n + m)d + cos 2(n — m)d| dd. 


Also from (1.24) and (2.10) the momentum transport spectrum 1s 


u'w’ 


‘O } 


n 
> T,, sin 2nd, 
1 
where 
T 7 2 Mier Ms iit (2.14) 

It is now possible to write down a complete set of equations optimizing 
Y, for a minimum («R)min. However, not only is this a most complicated 
problem, but its solution will be as approximate as the stability criteria. 
Thus, to establish the more significant physical consequences of the extreme 
state, the spectrum }, will be approximated to roughly the same extent 
that ® is in equation (2.1). ‘This approximation of the optimum Y,,, will 
be done in two ways: First, only the consequences of Y,, at large my will 
be considered, so that in expansions in powers of 1/m,) only the first term 
will be kept. Second, rather than optimize each Y’,, separately, optimization 
will be made between sets of Y,, which satisfy the conditions (2.7), (2.8), 
(2.10) and (2.11). 

Although other trial functions may converge more rapidly to the 
optimum dissipation rate, a simple power series has many mathematical 


advantages. Hence )/, is expanded as 


i 4 ' c a ae Se = 7 
"= ToL” (uel2)” * (ml2e” * (alae *y (2.15) 


‘To have even one free parameter for a variational optimization, Y,, must 











532 W. V. R. Malkus 


have five free parameters. Therefore Y,, in (2.15) will be terminated at 
the n* term, leaving A, m, C,, C, and C, to be determined. From (2.15) 
and (1.26) one finds that 


A_[f, , sin(m+ 1b) Cr f (my, 1\sin(t + 1)o _ 
(2n,)!? | sin d f a (No 2) ‘(3 , ) sing 


_ 1 (1 —cos(m + 1)) C. (3 , sin(m, + 1) " 


























2 sin? d (ng/2)?|\2 sing 
(n) +3) cos(m)+2)6 1sin(m,+3)¢ 
7 Zz sin’d = 2s sinh 7 
s{(# i t I)p 3 (ny + +) onto I)d cos _ 
Ny|2 2° sind ’ 3 z sind 
n, + 1\sin(n, + 1)¢(1+ cos*d)  1(1+2cos*d) 
" 4 2 ye sin’ "4 sintd 
1 cos(m, + 1)¢ cos 6(5 + cos*) | 
~ 8 sintd |, (2.46) 


jn 

where the various sums, >} ncos2nd, etc., have been performed by 
the calculus of finite differences. In keeping with the requirement of 
large mp, it is convenient to divide d? into a function describing its complicated 
transition region near the boundary and a function which is important away 
from the boundary. ‘These are 


(d?ny), _ = : 


try 
tw 
ro Nh 
a 
re eee 
Bel ae 
2 | 
3.| 
Im 
im 
be cca 
te 
a 
nN 
a 
~ 








and 
AT  ., (1 — cos €) cos€é 2sin€é 
(d* Ny)4~ } ~ yEe2 Cy sin € + C i ae —>+¢ ? 2 ¢ ea —— a. 
a¢ é L § 
_{,cos€ sing (1—cosé MP os 
+ C343 ——— ae he en | ~ oe L(€), (2.18) 


where Cy = 1+C,+C,+C;, and €=n,¢. Hence a (2.7) becomes 
AL(E,) = (1+ A(1 =) 2E4,(40)/ Fito)? (2.19) 


Moreover, if one neglects the variations of [(1 + A)(1 — y)]®? with €,, equation 
(2.8) becomes 


L'(é ; loFs (0) 
eel _. 2we. Mo, Co) 2.2 
Te) ~ Hew) |1+ re | aii 
Equation (2.12) can be summed to give 

2= —C,-—3C,-3C;, (2.21) 


to the zeroth order in 1/m), while (2.11) is summed as 
2= A*X1+ C, e 3(2¢ ‘2 +¢ ) i) HC, C2 +C3)+4 
+ 2(2C, C3 + C3) + $C, C34+4C3}. (2.22) 
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Then equations (2.6), (2.19), (2.20), (2.21) and (2.22) are five relations 
between six independent variables A, C,, Cy, C3, m9 and €... A choice for €,, 
say, determines the other five. Actually, A and y are small quadratic 
functions of these independent variables, from equation (2.3). However, 
to a zeroth approximation they may be assumed to be zero, their values then 
computed for each €, and these values used for a first approximation of 
A, C,, Cy, C, at each €.. 

Before resolving these algebraic relations the remainder of the problem 
will be discussed. From (2.9) we have 


nv | 
(*R) min = nk(22) (1—y)’ (2.23) 


i.e. C(1+A) = 7K €,. Hence the dependence of m)/K on 
f= E(A, C, C2, Cs, No) 

must be found in order to determine that €.(say €.) which makes («R)min 
a minimum. Following this, to determine the relation between my and R, 
% must be determined. Now «, the down-stream wave number, bears 
some (probably fixed) relation to }zm,, the cross-stream wave number, say 

lin, = ra. (2.24) 
(If r = 1, the smallest ‘eddy’ would be circular; if r = 10 the smallest 
‘eddy’ would be a 10 to 1 ellipse with major axis parallel to the boundary.) 
Then from (2.6) and (2.23) it follows that 


- } 3 
Rmin = (1 — (2) (32) mo (2.25) 





Hence a determination of the ratio r establishes the relation R = R(n)). 
‘To interpret the various formulae in conventional experimental terms, 
it is convenient to define a ‘boundary Reynolds number’: 
sl, 

Rg (2.26) 
where s is a distance proportional to 3 /m, say s = 39/%, and U, is the 
velocity that would be due to the initial gradient at that distance from the 
boundary. So defined, equation (2.26) may be written as 








a*RK % | n e.1u 
fe ae Se ees ‘ie — u. + = foe = 27 
‘ge ° Sly © gage 
where L_ v7, is the conventional ‘friction velocity’, in terms of which 
most data have been described. From (2.27) and (2.25) we have 
at = no\3 
p> ees =]. 2.28 
i aeita sini 


Finally, to find m)/A and the velocity profile as functions of 
£. = &.(A, CL Cy Cs), 
equations (1.23) and (1.26) may be integrated using (2.17) and (2.18). From 
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(1.26), (2.17), and (2.18), it follows that 





7 = | ( ‘ *) Pick < ” G ; ~ )(d?n) _, dO + 
haces’ | &2 \ /d? 
J) E-seMeda® = O 
where 6 = }7—4¢, and 1/n,< € < 1. For large m, 
my|/K = $A?(hCi(In Png + 1) +7?7/12+C], (2.30) 
where I’ = 1-78* is Euler’s constant, and 
Cn EE 4 ed - )2@C,—3 +4C,) (2.31) 


after some condensation with the help of (2.21). From (1.25) and (2.18), 
the velocity field near the boundary may be determined in terms of sine and 
cosine integrals: 


n a ee * (d* 
— ll = I\n€—€ | (<) ds+ | (¢) ds. 
K 0\"Mo]/ 6~1 Ny -0 \Xo d~ln 


After rising linearly, as }7é, from the boundary, , ///K turns sharply 
through a ‘transition’ region to its asymptotic value 


my U/K = $A*{3CR(In 20 + 1) 4+7?/12+C]. (2.32) 
From (1.23) and (2.17), direct integration gives 
22 (Minas HM) = 1.42[162 —10C?Incos8]. (2.33) 


Now combining (2.30), (2.32) and (2.33) with (2.27), one may write 
a theoretical form for the gross features of channel flow bevond the 
‘transition’ region: 


l 3) (sg—3)U,. / aI 
a \ Ry 4 In puna (1 + In a : 








Rp 

2 (7 
i essen, Seen ay Bu 1 ? 36 
als c) | (2.36) 


Those readers familiar with the experimental literature will have recognized 
a ‘velocity defect law’ in equation (2.35) and a ‘logarithmic boundary law’ 
in (2.34). The usual empirical constants in these experimental laws 
depend here on the determination of € = €,(A, C,, Cy, C3, m9) and the ratio r. 

Returning to the algebraic problem of finding €, = €,(A, C,, Cy, C3, no) 


from (2.6), (2.19), (2.20), (2.21) and (2.22), one notes that, for a particular 
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—, 7m is determined from (2.6) and table 1. Resolving the other four 


equations, it is found that 
2 

Y ~ €3\ 7 : ’ y y ’ ‘ 

C,= -2-(3 + #)e4 Q; C, = OG,-G,; C3 = O+e, Cy, 


where ; 2 g. p12 
. A 41 4) 41 














° oe Zz 
q, = -—¢ G, e,Gr- 420 a ee” 
a ee ee ee 
qo = Jg Tg Gy 71 U9 15 e’ 
Y . Ze 
3 = —¢,G,-e,G3-3+4; 
q3 & 2 ey 2 3 e? 
also e,—1 €1€,—e 
G, = , ‘ G, = ++, 
Cs — 4 ae 
| (LEA =P ECF)? 
, (d;—3d;) 
ey (2d, — dy) a (3d, —d,) 
a” (ds, 5d,)’ _ (d; - 5d) 
_f 2 7 ll , - 7 11 
«= (ist imot a4) = (ign t 2H0%) 
1 2 pies 5 F(y)[1 + oF (M0) (ds — 341) 
(6 = te = — és J, i ; ; 
> 2a SI (o)(a; — $d)) 
i. (2d; — dj) si (3d; —ds) 
2" G-iM)y *~ @-i4)’ 
and 


LIE.) — (dy +C,d,+Cyd,+Cyds), L'(é,) = (d+ C, d; + Cyd, + Cyd;). 
For each chosen €,, the quantity (m/K)(/€,.) of equation (2.23) may be 
found from these relations and (2.30). However, the restriction that the 
first interception of J,(7) and 7(c), equation (2.4), determines &, sets an 
upper limit on the €, that may be chosen. ‘Table 2 lists the computed 
values for (m/K)(j/&)), the corresponding values of y and A, and the 


maximum value of 





I = AL(E)E'*{[1 + A(E)] . [1 — 1(E)]} 7° */[4 (m0) ]"? E< So 

at several €, on either side of this upper limit. When ‘ = 1, J,(n) marginally 

| | | 
é. (2 2) | (72 7) | y A I(€) é 
I\K ° &,J/n,=102 | K * &¢J/ng=10* | 
6-00 14-87 30-14 | 0-095 | 0-045 | 0-846 | 0-937 
| 350° — 6-1086 14-40 28-78 0-121 | 0-065 | 1-042 | 0-977 
27 13-84 27-80 0-183 | 0-114 | 1-375 | 1-057 
t 
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intercepts v(c) at the smaller €. Hence, froma graph of table 2, the minimum 
value of (m)/K)(/€.) occurs at €, = €,= 6-081. Substitution of this €, 
into the relations (2.37) gives _ 

C, = </, O = —1-630, A*® = 11-43; 


and Ho — 2 . ; 
2 = 04875, CZ = 2308, 9 C+ T= - 1-908, (2.38) 
for y 0-114, A = 0-059. 

The ratio r can be determined from the second relation imposed by the 
boundary conditions on the asymptotic ®. (Equation (2.4) was the first 
relation.) However, this is a difficult task involving sums of series of 
integrals of functions of the optimum velocity profile. Instead, an integral 
condition for stability due to Synge (1938) will be used. Synge observed 
that multiplying the Orr-Sommerfeld equation (1.11) by ®* dd, integrating 
from 0) to 7 and adding the resulting equation to its own complex conjugate 
leads to J?24. 292]? + 472 = —aRM'(0)I,—2Re(1? + 22), (2.39) 

- - 
where [> - Z ‘beds, I? = | O'*'dd, /? = @’)*’ dd, 
~ O 7 ~ 0 
bi | ( /1'/11'(0))(@'D* — O*") dd, 
and the agers denotes differentiation with respect tod. When ® is stable, 
we have ¢, ; hence 


2 is vO) wy s@ @p*’ 
1/03 | - ent - Fs): 


2 (bye 1, \2 
4 =(Far) |oo* dd +( 3 =] > 0. (2.40) 
(2 


A sufficient condition for stability is that the integral in (2.40) be greater 
than zero. If then ®’/® = 4+7B, this integral is positive definite when 
Ke 424 B)bb*dd i) ‘ddd 
['(0)RP = 83M, M=— . (2.41) 
LF H’/ 71 (0)) Booed} 
Synge noted that, since [B%bd* db |d* ds > [| BOd* dg)’, M in (2.41) 
must be greater than one. Hence, [//'(0)R|?=8<" sets a lower limit 
on x for stable ©. However, a closer condition can be set by using 
the asymptotic ® toestimate MW. Nowatthe boundary, ¢ = 0, A(0) > F(%) 
and B(O) = Fy). From (2.38) and table 1, F',(y)) = 0-465 and 
F.(m)) = 0-174. 

Also from table 1 one sees that F',(y)/F,(y9) increases for smaller 7 (i.e. 
larger ¢) through the range where //’///'(0) is large. Hence, in a conservative 
estimate, we have 











ar 4017 Pa 
A2+ BYbO* dd > 114] B21 | Bowe dé, 
i- ) 6 L Eo JJo $ 
/0-465\2 
and M [1+ (saz), | =8. 


Since the value of M will enter final formulae only as M1, this represents 
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only a slight improvement on Synge’s condition. An extension of the 
preceding analysis shows that J,/a?/J, and J,/«J, approach zero for large x 
and larger R; hence using (2.41), we have 


o 2\77) - 
C; > -[#(1-3) | > 0. 


Equation (2.41) then becomes as accurate an estimate of the marginally 
stable x as one could hope to attain by the tedious method mentioned earlier. 
Finally, equating the R in (2.41) and (2.23), one obtains 


¢ 
r=(1—y)! »(=)(smy 6 — 3-94 (2.42) 
‘lo 


for M = 8, m/§ = 0-4875 and y = 0-114 from (2.38). Inserting the 
numerical results of (2.38) and (2.42) into (2.28) and the constants of (2.34) 


and (2.35), one obtains 


os Ff No\3 

=o 45 (¢) =127 2.43 
Ry rt eee (2) 71, (2.43) 

AO x 
‘ Ry (— =3-014; ev. =3-04 041, (2.44) 

71" 2 [x 

+n see ) + cea +0 fe + POR: Vv. =1:8+0- 2.45) 
(1 n= ) Alp c) + 1-001; e.v.=1-8+0-1, (2.45) 
where e.v. means the experimental values recently determined by Laufer 


(1950). 

The constant of equation (2.44) (inversely proportional to von Karmin’s 
constant) is the most important in these results since it enters into every 
other relation. The ‘logarithmic intercept’ constant, equation (2.45), 
involves the relatively sensitive quantity C (equation (2.31)), perhaps 
accounting for the error in this determination. An observation of r will 
be discussed inthe conclusion. ‘The dissipation rate is inversely proportional 
to the square of U,,,/U_, equation (2.36); hence it is larger than observed 
due to the error in C. However, at large R the percentage error due to 
C is small, and the dissipation rate approaches the observed curve. 

It should be emphasized that an extension of the technique outlined 
here involves considerable computational labour. In addition, further 
variational freedom for the Y,, spectrum may necessitate the use of Rmin 
rather than the approximate («R) min. 

To summarize this section: All known gross qualitative features of 
turbulent flow in smooth channels result from an optimization of the 
dissipation rate subject to the mean stability constraints ; a decisive quanti- 
tative comparison with the observations is not achieved primarily due to 
the mathematical difficulties encountered in finding a general solution 
to the Orr-Sommerfeld equation. 


3. CONCLUSION 
From a mathematical viewpoint the foregoing analysis is most laborious. 
The asymptotic solution to the Orr-Sommerfeld equation for a quite 
general U is described. Only then is U restricted to an optimal form, 
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However, the stability problem might be considerably simplified if the 
major restrictions on U imposed by the optimization could be made first. 
Some slight progress has been made along this line by performing con- 
strained variations of U in the integral stability equations (2.39). This 
approach offers hope of a more complete understanding of the problem 
than the maze of formulae employed here permits. 

On the physical side this work strongly supports the basic mechanistic 
assumption. Its validity can be tested by measuring the smallest scales of 
motion in the transport and energy spectra. Laufer (1950) has made 
time-spectra studies in channel flow from the boundaries to channel centre. 
The smallest scale of motion he detected in the energy spectrum was in 
the middle of the ‘laminar’ boundary region, and was twice as large as the 
observed thickness of this region. Since the energy spectrum has twice 
the maximum wave number of the velocity spectrum, this observation 
establishes that the minimum ellipticity of the smallest scale of motion is 
approximately four to one. Equation (2.42) gives a theoretical optimum 
value of 3-94. A more successful investigation of the detailed stability 
restrictions on the spectral transfer of energy may one day remove the 
assumption entirely. 

Before leaving the topic of spectra, it should be noted that the spatial 
spectrum for the transport given by (2.14) differs in kind from the temporal 
spectra obtained by the experimentalist. In homogeneous turbulence, 
spatial and temporal spectra would be identical; but in this inhomogeneous 
and non-isotropic shear flow, the observed spectra are functions of position. 
The work presented here does not determine the temporal spectra and their 
spatial dependence, but does provide a framework for such an inquiry. 

The apparent achievement of a statistical extreme state by the turbulent 
fluid raises questions of a different sort than those raised by the mechanistic 
assumption. It suggests that one observes the most probable of a large 
set of possible motions: that turbulent solutions to the equations of motion 
may be highly degenerate. If this is so, what determines the observed 
spatial and spectral structure of statistical observables, such as y’?, not 
studied here? One might think that degeneracy implied the equal a priori 
probability of each solution, but the infinite series of higher order small 
corrections to the usual equations of motion remove degeneracies and 
couple the macroscopic field of motion to the underlying microscopic 
particle motion. Hence it is proposed that the probability of occurrence of 
a degenerate solution is directly proportional to its degree of disorganization, 
or entropy. ‘The search for a useful measure of the disorganization of 
irreversible fluid processes involves sufficiently different formal techniques 
to necessitate separate discussion. However, even the usual thermostatic 
statistical mechanics suggests that a steady-state system with a fixed number 
of particles would have maximum ‘entropy’ when it has a maximum 
‘temperature’ and when all the dynamic variables have, as nearly as possible, 
normal Gaussian distributions about their mean values. In steady-state 
shear flow, the mean temperature is a maximum relative to the boundary 
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temperature for a maximum rate of dissipation of potential energy into heat. 
Yet to achieve this first extreme the mean momentum transport must be 
highly organized, as studied in this paper. Hence a completely random 
distribution of either macroscopic or microscopic dynamic variables is not 
possible. A study of possible secondary extremes for partially constrained 
macroscopic observables of turbulent shear flow will be presented in another 
paper. 

Finally, it is of interest that the techniques used here are of value in the 
study of other forms of turbulence. Many of the laminar stability equations 
developed for geophysical and astrophysical problems can serve as the 
foundation for determining the gross statistical features of their corre- 
sponding turbulent fields. Certain of these stability problems are much 
easier than the Orr-Sommerfeld problem. In particular this is true of 
turbulence due to thermal convection. With heat flux replacing momentum 
transport, mean temperature profile replacing the mean velocity profile, and 
Rayleigh’s equation replacing the Orr-Sommerfeld equation, turbulent 
convection may be treated just as shear flow was here. ‘Two earlier papers 
(Malkus 1954) report the initial results of such a study. 


‘The author wishes to thank Dr G. E. R. Deacon, Director of the National 
Institute of Oceanography, Surrey, and Professor P. A. Sheppard of the 
Department of Meteorology, Imperial College, London, for their kindness 
to him during the year in England in which this work was undertaken. 
He is particularly indebted to Professor C. C. Lin of the Massachusetts 
Institute of Technology and to Dr G. K. Batchelor of Trinity College, 
Cambridge, for their comments on the stability problem. ‘This work was 
performed under the auspices of the Office of Naval Research, and is 
contribution No. 804 from the Woods Hole Oceanographic Institution. 
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SUMMARY 

Small spheroidal particles suspended in a sheared viscous 
liquid are sometimes observed to take up slowly preferred orienta- 
tions, relative to the motion of the undisturbed liquid, which are 
independent of the initial conditions of release. ‘These observa- 
tions cannot be accounted for by the solution, obtained by Jeffery 
(1922), of the linearized Navier-Stokes equations. It is shown in 
this paper that the effect of the inertia of the liquid is to alter 
slowly the orbit of the particle in accordance with Jeffery’s hypo- 
thesis that the particle ultimately moves in such a way that the 
dissipation of energy is a minimum, but that this effect is orders of 
magnitude too small to account for any of the experimental 
observations. 

It is suggested that non-Newtonian properties of the liquid 
account for the observations. It is shown that the rate of orienta- 
tion of a particle would then be independent of its size, and this 
prediction is verified experimentally. Other experimental evi- 
dence in support of this suggestion is also described. 

Some remarks are also made about the possible effect of 
collisions between the particles when more than one particle is 


present. 


1. INTRODUCTION 
‘The slow motion of a small ellipsoidal particle in a uniformly sheared vis- 
cous fluid was examined by Jeffery (1922). The usual approximation was 
made that the acceleration (commonly called inertia) terms p(¢u/ot + u.Vu) 
in the hydrodynamical equations of motion may be neglected when the 
Reynolds number is sufficiently small. ‘The approximate equations of 
motion 


pV2u—Vp = 0, div u = 0, (1) 
were then solved with the boundary condition that u = Qar at the surface 
of the ellipsoid. Here, u, p, « and p denote respectively the velocity, 
pressure, viscosity and density of the liquid, Q is the angular velocity of the 
ellipsoid (to be found) and the origin of r is the centre of the ellipsoid. The 
angular velocity was determined by the condition that the total force and 
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couple on the body should vanish; this condition is valid since, as can be 
shown, the inertia of the body may be neglected to the same order in the 
Reynolds number as the approximation involved in (1). 

The equations that determine the motion or orbit of the ellipsoid can 
be integrated when the ellipsoid is one of revolution. Let @ be the angle 
between the axis of symmetry of the spheroid and the direction of the vortex 
lines of the undisturbed shear, and let ¢ be the angle between the plane that 
contains the axis of symmetry and the direction of the vortex lines and a 
fixed plane containing the latter direction, that is, 6 is the azimuthal angle 
of the axis of symmetry about the direction of the vortex lines. ‘Then the 
orbit of the spheroid is given by the equations 


(a2 +b?)@ = x(a?—b?)sin@cos@sindcosd, | 

(2) 
(a? +b?)b = x(a* cos*¢ + b sin*4), 

which have the solution 


tand = - tan —; 


ab 





moe k(a? cos? + b?sin2)!? ’ 

where 2a is the length of the axis of symmetry, 20 is the length of the other 
principal axes, « is the magnitude of the undisturbed shear, and &, the 
orbital parameter, is an arbitrary constant of integration that depends upon 
the initial conditions of release. 

These equations show that the motion is periodic with the axis of sym- 
metry describing a (non-circular) cone about the direction of the vortex 
lines, so that the motion is kinematically similar to that of a top. It also 
follows that the value of k, which specifies the particular orbit, is deter- 
mined by the initial conditions, and there is no tendency for the spheroid to 
set itself so that its axes point in a particular direction or to choose a parti- 
cular orbit, that is, there are no preferred orientations relative to the motion 
of the surrounding fluid that the spheroid will eventually take up irrespective 
of the initial conditions. 

Jeffery did not consider the existence of this indeterminateness to be 
entirely satisfactory and expressed the opinion that a more complete investi- 
gation would reveal the existence of preferred orbits. In effect, the sug- 
gestion made was that over a short period of time the orbit is given accurately 
by equations (2) and (3), but that effects neglected in the analysis have a 
cumulative effect and slowly alter the orbit, so that eventually the particle 
is moving in a preferred orbit which is independent of the initial conditions. 
Jeffery suggested further that these preferred orbits would correspond to 
those solutions of the approximate equations of motion (1) for which the 
rate at which energy is dissipated by viscosity is a minimum. It was then 
shown that, according to this hypothesis, prolate spheroidal particles would 
tend to set themselves with the axis of symmetry parallel to the vortex lines 
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and would rotate about this axis with constant angular velocity, so that the 
whole motion is steady; and oblate spheroidal particles would tend to set 
themselves with an equatorial axis parallel to the vortex lines and would 
rotate about this axis with a variable angular velocity, so that the whole 
motion is periodic but not steady. 

Several experiments have been performed in order to test the accuracy 
of equations (2) and (3) and to investigate whether preferred orbits do 
indeed exist. ‘The results of these experiments, which will be described 
briefly in $2, are not uniform. In some cases, preferred orbits do not 
appear to exist ; in other cases, they do exist but are orbits of either minimum 
or maximum energy dissipation. 

This question of preferred orbits is not without practical interest; for 
example, it is closely related to the viscosity of suspensions of non-spherical 
particles. It is the purpose of this paper to discuss some of the effects 
neglected by Jeffery in order to see whether there are any theoretical reasons 
why particles should take up preferred orbits in accordance with Jeffrey’s 
hypothesis, and to suggest a possible explanation of the experimental results. 

In § 3 we shall discuss the problem generally and show that many effects, 
such as those due to the inertia of the particle or the presence of rigid walls, 
need not be considered; and that the two most likely effects are the inertia 
of the liquid, which gives rise to the acceleration terms in the equations of 
motion, and non-Newtonian properties of the liquid. The rate at which 
the orbit of the particle is altered by the inertia of the liquid was calculated 
and the results are given in §4. ‘These results provide theoretical support 
for Jeffery’s hypothesis, but are unable to account for any of the experimental 
observations since the predicted rates of orientation are orders of magnitude 
toosmall. It is suggested in § 5 that non-Newtonian properties of the liquid 
are responsible for the observed phenomena, and experimental evidence in 
support of this suggestion is described. 

When more than one particle is present, as is the case for a suspension, 
there is a further effect to be considered, namely that of collisions and inter- 
actions between the particles (see Manley, Arlov & Mason 1955). In §6, 
we shall give the results (mainly negative) of calculating the mechanics of 
the interaction between two particles whose separation is always large 
compared with their size, and discuss briefly the significance of these results 
and some intuitive ideas about the problem. 


2. RELEVANT EXPERIMENTAL RESULTS 

Experiments using prolate and oblate spheroids were carried out by 
Taylor (1923). Binder (1939), Trevelyan & Mason (1951), and Manley & 
Mason (1956) describe experiments in which cylindrical rods of large 
aspect ratio were used (the aspect ratio is the ratio of the length to the 
diameter). In all these experiments, the particles were suspended in 
liquids contained between rotating concentric cylinders. ‘This motion 
differs from the uniform shear of infinite extent envisaged by Jeffrey in the 
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existence of Coriolis forces, wall and end effects, and the fact that the shear 
in the neighbourhood of the particle is only approximately uniform. How- 
ever, these effects can easily be shown to be negligible to the order considered 
by Jeffery, although their possible influence on the existence of preferred 
orbits has to be kept in mind (see §$3 and 4 below). 

Taylor used aluminium spheroids in water-glass. He observed that, 
as far as he could tell, the spheroids moved as predicted by equations (2) and 
(3), and that, after about 20 minutes for the prolate spheroids and five 
minutes for the oblate ones (these figures correspond roughly to 180 com- 
plete revolutions for the prolate spheroids and 40 for the oblate ones), the 
spheroids took up orbits in accordance with Jeffery’s hypothesis. In these 
experiments, « = 1-65 sec™!,e = 0-5 when a = 0-13 cm, and € == — 0-75 when 
a-=0-05cm. Here ¢€ denotes the ratio (a—b)/a, sometimes called the 
‘ellipticity’ of the spheroid. The viscosity of water-glass is very roughly 
1,000 gm/cm sec. 

Binder used particles of human hair suspended in glycerine. ‘The 
orbits of the particles are described as being in qualitative agreement with 
equations (2) and (3), although there must be some discrepancy since these 
equations describe the motion of spheroids and not cylindrical rods. When 
the aspect ratio of the rods was less than 15, the rods orientated themselves 
along the vortex lines, in accordance with Jeffery’s hypothesis, but when 
the aspect ratio was greater than 15, the rods set themselves perpendicular 
to the vortex lines in an orbit of maximum energy dissipation. 

Trevelyan & Mason and Manley & Mason used cylindrical rods of 
glass, dacron, nylon and human hair incornsyrup. They verified equations 
(2) and (3) to a reasonable degree of accuracy considering that rods and not 
spheroids were used, but their observations of the existence of preferred 
orbits for a single particle are inconclusive. ‘They sometimes observed a 
slow tendency for a particle to align itself along the vortex lines or to take up 
an orbit in which it was perpendicular to the vortex lines, but in other cases 
the orbit varied erratically with time. The manner in which the orbit 
varied appeared to be random and independent of the particle size, the 
aspect ratio, and the magnitude of the undisturbed shear. 

Manley & Mason carried out also experiments with suspensions of 
particles of large aspect ratio. Suspensions in which the particles were 
initially distributed isotropically were sheared, and a statistical tendency 
was found for those orbits with large values of k to be favoured, that is, 
there was a tendency for particles to orientate themselves along the vortex 
lines. 

3. "THE NON-LINEARITY OF THE PROBLEM 

We shall now consider whether a theoretical treatment which uses the 
approximate equations (1) can explain these experimental observations or 
predict the existence of preferred orbits. In other words, when the per- 
turbations due to the effect of rigid walls in the experimental apparatus, or 
the possibility that the shear is only approximately uniform, or any similar 
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alterations in the idealized external conditions envisaged by Jeffery, are 
taken into account using the approximate equations (1), does the orbit of 
the particle remain constant in time or alter slowly until a preferred orbit is 
reached ? 

Now equations (1) are linear and hence all pressures and stresses are 
proportional to the magnitude of the undisturbed shear and are therefore 
linear in x. It follows that the equations determining the orbit are linear 
in x, i.e. 6 and ¢ are given by the products of « and functions of 6 and ¢. 
‘The extra effects can be regarded as giving rise to perturbations of the orbit 
given by (3), and any variations of the orbit with time are then best repre- 
sented by changes in the mean value of the orbital parameter k, where k can 
vary between 0 and «, being zero when the axis of symmetry is perpendi- 
cular to, and infinite when parallel with, the vortex lines, and the mean is 
taken over a complete revolution of the particle. On putting 6, as given by 
(3), in the expressions obtained for # and ¢@ when the extra etfects are taken 
into account, and averaging over a complete revolution, we then obtain 


_ = « fl), (4) 


where f(k) depends on the shape but not the size of the particle and (as in all 
subsequent equations of this kind) & is to be understood as representing the 
average value over a complete revolution. Hence, if the shear is reversed, 
by, for example, the cylinders being rotated in opposite directions, then 
dk dt will change sign; and if a preferred orbit exists in one case, a different 
preferred orbit will arise when the shear is reversed. 

It follows immediately that, if the motion is such that the physical 
situation is unaltered by a reversal of the shear, then f(k) must be zero and 
the particle rotates in a constant orbit. In other words, the use of equations 
(1) will predict that the orbit remains fixed unless there is a certain degree 
of asymmetry in the external conditions. (Another way of seeing this 
result is to note that, according to equations (1), the motion is reversible.) 
This result was indeed verified directly for the case of a spheroid near a 
rigid plane wall. ‘The wall effect was calculated by a method due to Lorentz 
(1907), and the necessary modifications to (2) and (3) were found. It was 
found, as expected, that the mean value of & over one revolution remains 
constant. 

Further, this argument also shows that the inertia of the particle may be 
neglected since the equations for the orbit, with the inertia of the particle 
taken into account, are linear in kx. 

It is clear now that the question of preferred orbits is to some extent 
dependent on the design of the experimental apparatus. However, the 
Couette flow between rotating cylinders that was used in the experiments 
described in the previous section and with which we are primarily concerned 
possesses the necessary symmetry for the above result to be valid. Further, 
it is to be expected that flows for which the symmetry condition is not valid 
will be the exception rather than the rule; and even in these cases, the result 
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may still hold since the linearity of equations (1) and the consequent 
superposability of solutions enables the motion to be split into two parts, 
one of which is reversible and the other not, and the part that is not reversible 
may not produce a couple onthe particle. For example, if account is taken of 
the relative velocity of particle and liquid owing to their different densities, 
this relative velocity does not change sign when the shear is reversed but 
the couple it produces on the particle is zero. 

We must therefore look for non-linear terms that are neglected in (1). 
The first such terms that come to mind are the acceleration (inertia) terms 
and in the next section we shall examine their effect. We shall see, however, 
that the experimental phenomena cannot be explained in this way and some 
other non-linear effect must be sought. _Non-Newtonian properties of the 
liquid give rise to non-linear effects and these will be discussed in ‘§5 to- 
gether with some relevant experimental results. 

There is, of course, the possibility that external effects of a random nature 
influence the orbits of the particles. It is difficult to understand the experi- 
ments with single particles by Manley & Mason (1956) in any other way. 
Indeed, they say that slow convection currents were present and regard 
these as the explanation of the randomness of their observations, and the 
opinion is expressed that the orbit would remain constant in the absence of 
such random effects. However, experiments that are capable of being 
repeated, such as those of Taylor or those described in § 5, cannot be explained 
in this way. 

4. 'THE INERTIA OF THE LIQUID 
The full equations of motion of an incompressible liquid are 
Du 

P Dt 

To calculate the effect of the acceleration terms by obtaining the complete 
solution of (5) is, in general, a problem of great difficulty, but we can by 
simple dimensional considerations obtain the order of magnitude of their 
effect. For p(Du/Dt)/uV?u is of order R, where R denotes the Reynolds 
number xa?/v (v = u/p), and on making the plausible assumption that the 
extra stresses introduced when the acceleration terms are taken into account 
are smaller by a factor of order R than the stresses calculated neglecting 
these terms, we obtain for the rate of change of the orbit due to the inertia 


1 dk are (k 
= 7 s( ), (6) 


—Vp+pV2u, divu = 0. (5) 


of the liquid 





k dt v °\a 
where g(k/a) depends on the shape of the particle, and the factor «(= (a—b)/a) 
has been introduced to emphasize the fact that dk/dt vanishes from consider- 
ations of symmetry when the particle is a sphere. 
Now the sign of g(k/a) is of interest in connection with Jeffery’s hypo- 
thesis, and it is also not without value to have an estimate of its numerical 
It was therefore decided to calculate this function by taking 


magnitude. 
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account of the acceleration terms in (5) by means of an iterative procedure. 
‘That is, denoting the solution of (5) with pDu/Dt neglected by u™ and p”, 
where u = Qar at the surface of the spheroid and u tends to the 


undisturbed velocity as r —- © (this solution is, of course, the one found by 
Jeffery), we consider the equations 

Du : ™ 

uV2u—Vp = Pr divu = 0, (5’) 


and try to find a solution which satisfies the boundary conditions u = Qar 
on the spheroid and u— u'? —- 0 as r—> a. 

A most serious difficulty now arises since it can be shown that no solution 
ot (5’) can be constructed which satisfies the boundary condition at infinity : 
this is sometimes known as Whitehead’s paradox and was first encountered 
by Whitehead (1893) when considering the uniform translation of a sphere. 
However, there are reasons for believing that this difficulty can be overcome 
in the following way: there is evidence that if a suitable boundary condition 
at infinity is chosen for u, then the solution of (5’) that satisfies this boundary 
condition will represent accurately in the neighbourhood of the particle the 
second approximation to the solution of (5). This boundary condition at 
infinity is deduced from the solution of the Oseen-type equations for the 
problem, but in the present problem these equations are difficult and have 
so far not yet been solved. However, by analogy with simpler cases, it 
appears that the requisite boundary condition is |u—u‘ —ARu™| = 0 (r) as 
r—» «x, where A is a numerical constant (possibly zero) of order one. Equa- 
tions (5’) can then be solved, the second approximation to the couple on the 
particle can be calculated in terms of 0, « and A, and the orbit of the particle 
can be determined. It is clear that the effect of A is equivalent to a small 
change in the magnitude of the undisturbed shear and so, to the order con- 
sidered, A will not enter the expression for dk/dt and its actual value is not 
needed for the purpose of this paper. It must be emphasized, however, 
that this boundary condition may well be wrong and this doubt is to be kept 
in mind whilst interpreting the results obtained below.* 


Result of the calculation 

Although it seems possible, in principle, to carry out the iteration for a 
spheroid of arbitrary shape, the analysis that appears to be required is such 
that no practical method could be found. ‘The analysis was therefore con- 
fined to the comparatively simple case of a nearly spherical spheroid, so 
that powers of e higher than the first could be neglected. The solution was 
then obtained by expanding the velocity and pressure in series of solid 
harmonics multiplied by powers of r. Since the algebra is tedious and 


* I am deeply indebted to Dr I. Proudman for information and advice concerning 
the validity of iterated solutions of the Navier-Stokes equations and about the 
boundary condition to be satisfied at infinity. It is regretted that a description of 
the reasons for suggesting the above boundary condition would be quite lengthy 
and out of place in the present context. 
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involved, it would serve no useful purpose to reproduce it here and only the 
final result will be given. 

It was found that when the acceleration terms are taken account of in the 
above manner, 
dk a*x*€ ? 
= 0-24 — (7) 
This result is equivalent to g(k/a) having the value 0-24 + O(e). 

For a prolate spheroid, « ~ 0 and so dk/dt > 0, i.e. the long axis tends to 
set itself in the direction of the vortex lines; for an oblate spheroid, « < 0 
and so dk/dt > 0, i.e. the short axis tends to set itself perpendicular to the 
vortex lines. Thus, the effect of the inertia of the liquid is in accordance 
with Jeffery’s hypothesis and provides theoretical support for this hypo- 
thesis. 





al 


Comparison with experiment 

Let us now compare the observed rates of orientation with those pre- 
dicted by equations (6) or (7). ‘The relevant experiments are those of 
Taylor (1923) (these experiments were repeated by myself and similar 
results obtained). It is true that € is not sufficiently small for (7) to be 
accurate, and that Coriolis forces of the same order of magnitude as the 
acceleration terms should be taken into account, but (7) should give the order 
of magnitude of the effect of the inertia of the liquid without serious error. 
(The ratio of Coriolis to acceleration terms is (2’/« or (2’/«)*, according as 
to which is the larger, where {2’ is the angular velocity of the particle around 
the axis of the cylinders ; this ratio is either of order one or zero in the experi- 
ments described in this paper.) On putting the data given by Taylor into 
equation (7), or using equation (6) with the plausible assumption that 
g(k/a) is of order one, we find that (1/k)(dk/dt) is at most 10 sec"'. Nowa 
change in the orbit should be observable when k changes by a factor e, and 
the time for this to happen is not less than 10° seconds == 28 hours! Since 
the observed times were not greater than 20 minutes, there can be little 
doubt that the effects of the inertia of the fluid are too small to produce the 
observed phenomena. 

When the aspect ratio of the particles is large, equation (7) is likely to be 
in serious error, and a direct comparison with the experimental observations 
on cylindrical rods cannot be made. However, a general examination of 
the equations for the case of large aspect ratio shows that (1/k)(dk/dt) is then 
probably proportional to b?x?(loga/b)?/v. If this be so, it is easily shown 
that the times needed for a change of orbit to be clearly observable are orders 
of magnitude larger than those observed by Binder (1939), and are also long 
compared with the duration of the other experiments on cylindrical rods 
quoted previously. 


The motion of small particles in a tube 
Of further interest is another phenomenon predicted by Jeffery. It 
can easily be shown that small particles suspended in the liquid flowing 
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through a capillary tube will, according to equation (1), remain at a fixed 
distance from the axis of the tube; however, it is a consequence of Jeffery’s 
hypothesis of minimum energy dissipation that they should move slowly 
towards the axis since the shear there is least, so the extra energy dissipation 
owing to the presence of the particles is then reduced. (‘This phenomenon 
appears to be quite well known, and is observed, for example, by physio- 
logists studying the motion of blood corpuscles in the capillaries of the body. 
I am not aware, however, of any accurate quantitative measurements that 
have been published.) 

‘The etfect of the inertia of the fluid on the motion of a small sphere of 
radius a in a flow with a parabolic velocity profile vu = «z+ Bs? was worked 
out by iteration as before and it was found that the inertia terms give rise to 
a sideways velocity w, perpendicular to the directions of the flow and the 
vortex lines, given by 


i 
w= —0-43 —. (8) 


This velocity moves the sphere to a region where the magnitude of the shear 
is least, and so the effect of the inertia of the fluid is again in accordance with 
Jettery’s hypothesis. Dr M. J. Taylor has, however, communicated to me 
some rough experimental data of his own, from which it appears that the 
velocities predicted by (8) are too small by a factor of order 10. 

‘lo sum up, if the inertia of the liquid is the dominant non-linear effect, 
then there is theoretical support for Jeffery’s hypothesis. It appears, how- 
ever, that in the experiments described in this paper, the inertia terms are 
too small to account for the observations, and we shall now go on to discuss 
the possible effects of non-Newtonian properties of the liquid. 


5. NON-NEWTONIAN PROPERTIES OF THE LIQUID 
It is assumed in the derivation of the Navier-Stokes equations (5) that 
the stress tensor p,; is linearly related to the rate-of-strain tensor 


by the equation Pijy = — pd, + he;;- (9) 


A liquid for which (9) does not hold is said to be non-Newtonian, and this 
definition obviously covers a wide range of liquids with various properties. 
However, we are concerned here with incompressible liquids which are 
almost Newtonian, that is, liquids in which anomalous effects caused by 
deviations from equation (9) are small. For such liquids, we may take the 
stress-strain relation as 


Pi; = — pd; + me; + aO,;, (10) 
in which « is constant and oQ,;, where o is a constant with the dimensions 
of viscosity divided by time which is some measure of the non-Newtonian 
properties of the liquid, is a quadratic functional of the rate-of-strain tensor 
which may involve both the values of e,; at earlier times and higher spatial 
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derivatives of e,;. It is also supposed that the motion is sufficiently slow 
for powers of e,; higher than the second to be neglected. 

The nature of the quadratic term in (10) is by no means clear. If it be 
supposed that the liquid is isotropic when at rest, that the stress depends 
only on the instantaneous rate-of-strain, and that the stress-strain relation 
is analytic, then it can be shown that oQ,; may be taken as p’e,, €,; + ue, 5,;, 
where »’ and yp” are constants (for a full review of the subject, see Truesdell 
(1952)). But in the present problem, where the motion is not steady, this 
form of the quadratic term may well not be valid. Apart from some evidence 
that .’ is positive if this form of the quadratic term is valid, there appears to 
be no non-Newtonian liquid for which the form or magnitude of the quad- 
ratic term has been established. 

The equations of motion of the liquid are ¢p,;/0x; = 0, since the inertia 
of the liquid may be neglected. Because the non-Newtonian terms are 
supposed to be small, a solution of these equations could be attempted by 
iteration, provided that a suitable form for O,; were postulated. (It is to 
be noted that the difficulty at infinity which arose when the acceleration 
terms were taken into account by iteration would not arise here.) It can 
be seen that, after solving the equations of motion and equating the couple 
on the particle to zero, we should finally obtain 

1 2g 
He = — hh a), (11) 
where /(k/a) depends on the shape of the particle and on the form of Q,;.. ‘The 
right-hand side of (11) might be identically zero, but there is no a priori 
reason why this should be so and it seems most unlikely that this should be 
the case. 

The rate of orientation of the particle is known when h(k/a) is known 
but no attempt was made to determine h(k/a) because the form of Q,; is not 
known beyond reasonable doubt. Even if oQ,; were known or postulated 
with a reasonable degree of certainty, there appears to be no practical way 
of solving the equations for a spheroid of arbitrary shape, and to carry out 
the analysis for a nearly spherical spheroid would require very heavy algebra 
and would probably not be of much use since oQ,; would probably contain 
sufficient unknown constants for agreement with experiment to be reached 
without difficulty. However, it follows from (11) that the rate of orien- 
tation of a particle of given shape in a given liquid is independent of its size ; 
whereas the rate of orientation due to the inertia of the liquid is, for a given 
shape, proportional to the surface area of the particle. ‘The following simple 
experiment was therefore carried out in order to see whether the rate of 
orientation is independent of the size; it was verified roughly that this is 
indeed the case. 


Experimental procedure 
The apparatus consisted of two concentric cylinders: the diameter of 
the outer one was 8-4cm and that of the inner one was 2‘8cm. Prolate and 
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oblate spheroids, similar in shape but of different sizes were made on a lathe 
from aluminium wire. ‘The longest and shortest axes of the smallest prolate 
spheroid were 0-25 cm and 0-10cm respectively, and of the smallest oblate 
spheroid 0:17cm and 0-10cm. ‘The other spheroids were twice and three 
times as large. ‘The outer cylinder rotated one revolution in 4:3 seconds 
and the inner one was kept at rest. ‘The particles were placed midway 
between the cylinders so that « was about 0-82sec"!. ‘T'wo liquids were 
used: a highly viscous form of ammonium polymethyl acrilate, and water- 
glass. ‘The densities of these liquids were measured in the usual way and 
the viscosities by timing the rate of fall of a small lead sphere of known size. 
The density of the ammonium polymethyl acrilate was close to that of 
water, that of the water-glass was about 1-5 times that of water; the kine- 
matic viscosities were, as it happened, approximately the same, being about 
1000 cm? sec™!. 

A common property of a non-Newtonian liquid is that, when confined 
between rotating cylinders, it will climb up the inner cylinder, provided the 
upper surface is free (see, for example, Weissenberg 1947). It was noticed 
that the ammonium polymethyl acrilate climbed up the inner cylinder by as 
much as one centimetre, whilst the water-glass did not climb up at all. This 
indicated that the non-Newtonian properties of the water glass were small 
compared with those of the other liquid. 

The experimental results for the ammonium folymethyl acrilate were 
as follows. ‘The prolate spheroids were placed with their long axis hori- 
zontal, and after about five minutes the long axis was vertical. The oblate 
spheroids were placed with their short axis vertical, and after about six 
minutes the short axis was horizontal. As far as I could tell from these 
crude observations, the rate of orientation was independent of the size of 
the particles. When water-glass was used, it was found that the orbits of 
the spheroids remained constant over times up to an hour. 

These results are consistent with the suggestion that the explanation of 
the observed phenomena lies in non-Newtonian properties of the liquid. 
There is an apparent discrepancy, however, between these results for water- 
glass and the observations of Taylor (1923). But Sir Geoffrey Taylor has 
informed me that, in his experiment, the water-glass rose some distance up 
the inner cylinder, and this was also the case when I repeated Taylor’s 
experiment with a different supply of water-glass and at probably a higher 
temperature than the experiments described above. Bearing in mind that 
water-glass is by no means a standard product, and that non-Newtonian 
properties are often critically dependent upon temperature (although the 
temperature dependence here is in the opposite direction from what would 
be intuitively expected), this apparent discrepancy appears to yield further 
evidence in favour of the above suggestion. Further, Dr S. G. Mason has 
informed me that the corn syrup used in the experiments in which he parti- 
cipated (where no preferred orbits were observed with single particles) was 
found to be Newtonian over a wide range of velocity gradients. 

There thus seems to be strong evidence that the existence of preferred 
orbits is related to the effect of non-Newtonian viscosity and that some of 
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the experimental observations of spheroidal particles suspended in the 
Couette flow between concentric cylinders will be explained when the action 
and effect of non-Newtonian properties can be calculated. 

Further, o is a measure of the elastic properties of the liquid, since it is 
easily shown that the quadratic term in (10) does no work when the liquid 
is deformed and returned to its original state, so that o/u is a measure of the 
relative magnitudes of the elastic and viscous properties of the liquid.* 
(c/ 1s probably proportional to x/y, where y is the elastic modulus of rigidity 
of the liquid.) ‘To account for the experiments described above, o/u for 
ammonium polymethyl acrilate would have to be of order 10-* sec. 

An experiment to investigate the possible relationship between the elastic 
properties of a variety of non-Newtonian liquids and the rate at which pre- 
ferred orbits are reached would be of great interest in this connection, as 
would also be an accurate experiment to check that the rate of orientation 
is proportional to the square of the undisturbed shear. 


6. SUSPENSIONS OF MANY PARTICLES 

When considering a suspension of many particles, we are primarily 
interested in the proportion of particles whose orbital parameters lie in a 
given range, since this is needed for a calculation of the bulk properties of 
the suspension and is also the quantity that is actually measured. ‘That is, 
we wish to determine the values of the function p(k) where p(k) dk denotes 
the proportion of particles with orbital parameters in the range (k, k + dk). 
If we were to follow a definite particle, we would find that its orbit is con- 
tinually being altered in apparently discontinuous jumps by collisions with 
other particles (see Manley & Mason 1956) and that there might also be a 
slow change due to the effects described in the previous sections. In this 
section, however, we shall not discuss these latter effects but confine our- 
selves to a discussion of the effect of collisions. (It is worth mentioning, 
in passing, that a suspension of non-spherical particles is anisotropic and 
this implies that the bulk properties of the suspension may be to some 
extent non-Newtonian. This in turn may effect the external motion—for 
example, a suspension contained between rotating cylinders may climb up 
the inner cylinder, even though the pure liquid does not—and so affect the 
motion of an individual particle and give rise to preferred orbits. ‘This 
effect will not be discussed here.) 

It does not necessarily follow, of course, that p(k) is altered by collisions, 
for example, it may be that the effect of a collision is to interchange the 
orbital parameters of the particles that collide; however, there is no a priort 
evidence that this is the case and it seems most likely that the effect of colli- 
sions is to alter the orbits until a statistically steady state is reached. It 
might at first be thought that, if the interactions were calculated using 
equations (1), the linearity of the equations and the consequent reversibility 
of the motion would imply that p(k) is not altered by collisions, since if 


* This was pointed out to me by Dr K. Weissenberg. 





552 P. G. Saffman 


dp(k)/dt is not zero, its sign is reversed by reversing the external shear. 
This is a well-known paradox of Statistical Mechanics, however, and it can 
in fact be shown that the reversibility of the motion does not imply that 
dp(k)/dt = 0, see ‘Yolman (1935) where this problem is discussed in con- 
nection with the Boltzmann H-Theorem. In other words, it is sufficient 
to calculate the collisions using the approximate equations (1). 

The possibility exists that the steady state could be calculated without 
a detailed knowledge of the mechanics of an interaction by means of the 
ideas of Statistical Mechanics. ‘This can be done if it is known whether 
there is any general property of the orbits that is conserved during a collision, 
something which corresponds, for example, to the conservation of energy 
for collisions between the molecules of a gas. With this idea in mind, the 
mechanics of an interaction in which d, the minimum distance of separation, 
is always much greater than a, a typical dimension of the particles, was 
worked out. Since if there isa quantity that is conserved for a general inter- 
action, it will be conserved for an interaction in which a/d is small, and it is 
possible to carry out the analysis for this case by neglecting powers of a/d 
higher than the third. 

The result of this lengthy calculation was that the arbitrary parameters 
that give the relative configuration of the particles before the interaction 
commences can be so varied that the orbits of the particles can be altered 
independently by the interaction. Further, it was found that the orbits of 
the particles are changed by the interaction in such a way that the statistical 
properties of the suspension are unaltered. Hence, there is no conservation 
principle, so that the ideas of Statistical Mechanics cannot be applied, and 
it is necessary to consider the detailed mechanics of a close interaction or 
collision before information about the effect of collisions can be obtained; 
unfortunately, there appears to be no way of doing this. 

Finally, there are two intuitive hypotheses about the etfect of collisions 
that are perhaps worth mentioning. ‘The first of these is that those orbits 
are favoured for which the energy dissipation is small, and the second is that 
the orbits favoured are those for which the volume swept out by the particle 
is small. ‘The first hypothesis is reminiscent of Jeffery’s, but the situation 
here is different from that of §4 since we are concerned with solutions of 
equations (1) and not with the etfcct of terms neglected in these equations. 
The second one, which was suggested to me by Dr J. W. Glen, is based on 
the intuitive idea that the smaller the volume swept out by a particle, the 
smaller the frequency of collisions made by the particle, with the result that 
the particle stays in its orbit longer. It follows from the first hypothesis 
that there is a tendency for particles to align themselves with the long axis 
along the direction of the vortex lines, so that p(k) is large for large k if the 
particles are prolate and for small k if they are oblate; and it follows from 
the second hypothesis that the tendency is for the axis of symmetry to lie 
along the direction of the vortex lines, so that p(R) is large for large k for both 
prolate and oblate particles. Thus, for prolate particles, the predictions 
are the same and in accordance with the observations of Manley & Mason 
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(1956), but the predictions for oblate particles are in disagreement. An 
experiment using oblate particles would be of some interest in this 


connection. 


I wish to thank Sir Geoffrey Taylor for suggesting this problem to me, 
and I am grateful to both him and Dr G. K. Batchelor for the kind interest 
that they have shown in this work. I also wish to thank Dr K. Weissenberg 
for the supply of some ammonium polymethyl acrilate. 


REFERENCES 
BInpER, A. 1939 7. Appl. Phys. 10, 711. 
1922 Proc. Roy. Soc. A, 102, 161. 


Jerrery, G. B. 
1907 Abhandlungen uber Theoretische Physik, p. 23. Leipzig. 


LORENTZ, H. A. 
MANLEY, R. St. J., ARLov, A. P. & Mason, S. G. 1955 Nature, Lond. 175, 682. 
MANLEY, R. St. J. & Mason, S. G. 1956 To be published. 

1923 Proc. Roy. Soc. A, 103, 58. 


TAYLOR, G. I. 
1938 The Principles of Statistical Mechanics. Oxford University 


TOLMAN, R. C. 
Press. 
TREVELYAN, B. J. & Mason, S. G. 1951 #. Coll. Sct. 6, 354. 
TRUESDELL, C. 1952 ¥. Rat. Mech. Anal. 1, 125 
WEISSENBERG, K. 1947 Nature, Lond. 159, 310. 
WuitEHEAD, A. N. 1893 Quart. J. Math. 23, 143 


F.M. 








REVIEWS 


The Structure of Turbulent Shear Flow, by A. A. Townsrnp. Cam- 
bridge University Press, 1956. 315 pp. 40s. or $7.50. 


Any steady flow of a fluid is unstable when the Reynolds number (product 
of a typical flow velocity with a typical length scale, divided by the kinematic 
viscosity) is great enough, although the required magnitude varies extremely 
widely with the type of flow and the size of such disturbances as are present. 
At Reynolds numbers just large enough to produce instability, the disturbed 
flow may settle into a regular eddy pattern (such as the Karman trail or 
Bénard cells), but at higher values an irregular, ‘turbulent’ flow is produced, 
in which the velocity at any point varies in a statistically random manner. 
Fourier analysis shows that fluctuations in a range of frequencies are present, 
the width of the range increasing with Reynolds number. 

The study of turbulent flow is a part of science not founded on the 
‘principle of uniformity’ as expounded by students of scientific method. 
The velocity record at a fixed point obtained by a hot-wire anemometer 
recording on to an oscilloscope is quite different every time a given flow 
(under fixed conditions) takes place. ‘The aim, throughout the seventy 
years since Reynolds initiated the study of turbulence, has been to recover 
some partial degree of uniformity by statistical studies, that is, studies of 
averages of various quantities. In theory, a fluid flow which is known 
precisely at one instant can have its later development predicted uniquely. 
Hence, if a complete statistical knowledge of the flows generated in a certain 
way were available for one instant of time, such knowledge could in theory 
be deduced for later instants. However, it is of the essence of the problem 
that the statistical distribution of the initial random disturbances that 
lead to turbulence is not known in practice. Hence, no definite results 
can be expected unless there is a tendency for the development with time to 
obliterate special effects of this initial distribution. Actually, some such 
tendency is observed, but it is only partial. 

Other branches of statistical physics have been more successful than 
turbulence theory. ‘lhe obliteration tendency appears most strongly in 
the kinetic theory of gases, but the ‘molecular chaos’, which implies a 
statistical independence of velocity distributions at different places however 
close, makes any analogy with turbulence tenuous indeed. ‘The existence 
of spatial correlation brings turbulence closer to that other rather unsuccessful 
subject, the kinetic theory of liquids. 

\gain, there is an obvious difference from linear problems where response 
is frequency-conserving. A better analogy is provided by forced mechanical 
vibrations whose amplitude is governed by non-linear forces which transfer 
energy between modes, such as could be observed in the film of the Tacoma 
Bridge disaster. However, turbulence at high Reynolds number is different 
in that it possesses a far greater range of wave-numbers, due to weaker 


internal damping. 
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This presence of modes of fluctuation with a wide range of wave-numbers, 
transferring energy among themselves as a result of the non-linear character 
of the inertial forces, is the basis of the Kolmogoroff—-Onsager theory of 
the ‘small eddies’, that is, of the components of flow which in a spatial 
Fourier analysis have large wave-number. ‘This theory is based on the 
idea that transfer of energy takes place in the direction of increasing wave- 
number, and in fairly short steps (in terms of the ratio of wave-number of 
a larger eddy to that of the smaller eddy which draws off some of its energy). 
Our knowledge of the basic instability of the steady flow already indicates 
that energy passes from that minimum-wave-number motion most readily 
to a motion with wave-number greater by half to one order of magnitude, 
and if this is true right up the wave-number scale the idea follows that the 
small eddies may have been reached by a statistical combination of many 
elementary processes which could obliterate the pecularities of the original 
large-scale disturbances. Hence they may be in an equilibrium which 
depends only on the viscosity which dissipates their energy and on the total 
rate at which energy is passed down to them, or (in the case of flows where 
this rate is changing with time) in a moving equilibrium which may depend 
also on the rate of change of that rate. 

‘This theory has given a fairly satisfactory statistical picture of the small 
eddies, close toisotropy and equilibrium, which is useful in many applications. 
However, because of the increase in rate of energy loss to smaller eddies 
and in rate of dissipation as wave-number increases, the actual energy in 
the small eddies at any instant is a small fraction of the whole. ‘The 
processes, if any, that govern the character of the ‘energy-containing 
eddies’ (those larger ones which contain the bulk of the turbulent energy) 
have remained a mystery, except that some progress has been made in the 
case of homogeneous turbulence (with a uniform mean flow), which was 
reviewed in 1953 by Dr G. K. Batchelor in the Cambridge series of 
Monographs on Mechanics and Applied Mathematics (this book contains 
also a good account of the small eddies). 

Now at last, in the same series, a book has come out which, largely 
eschewing the small eddies as in the main already understood, takes as its 
principal subject the problem of how, in general turbulent flows, the character 
of the energy-containing eddies is controlled. It is called The Structure 
of Turbulent Shear Flow, presumably to disavow any réstriction to the 
cases with uniform mean flow already described by Dr Batchelor, and 
perhaps also to exclude turbulent flows due to thermal convection, but it 
might have been called simply Turbulence. Its author, Dr Alan 
Townsend, has a combination of experimental and theoretical ability 
which no-one entering the subject has possessed since Sir Geottrey ‘Vaylor 
(in whose Cambridge laboratory he works, and whose encouragement he 
acknowledges in a preface), and which were almost certainly essential 
requisites for his task. He has been developing his ideas in a number of 
exciting articles over the last nine years, but it has not been easy to follow 
from these what the conclusions would be, as in many cases the views 
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expressed in one article were retracted in a succeeding one, and the author 
was plainly feeling his way. But now we have an opportunity to take in 
as a whole the theory and the experimental evidence on which it is based, 
both being very plainly set out in the book under review. 

A number of basic ideas dominate Dr ‘lownsend’s account of an 
extensive range of different flows. ‘The first three are simple general ideas 
which have been adopted fairly widely before, although not always, as 
Dr ‘Townsend points out, with a full understanding of which results follow 
from them and which from more specialized assumptions about structure. 
‘They say that at sufficiently high Reynolds number the energy-containing 
eddies (and, generally, all except the smallest-scale features of the flow) 
(i) are independent of Reynolds number (because the viscosity affects only 
the small eddies) and in time reach a state in which (11) their velocity and 
length-scales may be changing but other quantities of the same dimensions 
change in proportion, and which (iii) is independent of the details of such 
boundary conditions as are basically responsible for the turbulence, depending 
only on an overall resultant of them. ‘The last principle is close to that of 
St Venant in elasticity, and as in that subject it is the resultant force (together 
with, sometimes, a length scale) which matters. ‘hus, turbulence in a jet 
depends on its net thrust, in a wake on the drag of the body producing it, 
and in pipe flow on the wall resistance. ‘he independence of Reynolds 
number (principle (i)) and of details of geometry is postulated only for 
given values of these forces, which on the other hand for a given geometry 
may vary with Reynolds number owing to the kinematic boundary condition 
of no slip being applied in the ‘ viscous layer’ near the wall where the variation 
of mean velocity has a very small length-scale. 

Dr ‘Townsend begins his book with two chapters in which he sets forth 
the principal experimental and mathematical tools which are available for 
the study of turbulence; here, he departs from previous practice principally 
by stressing the idea of an eddy as something finite in extent (and even 
showing a preference for it being round, as indeed so many eddies observed 
in nature are!) rather than as being simply a Fourier component. In 
chapter 3 he sets out the evidence for the three principles just outlined in 
the case of turbulence with uniform mean flow, for which particularly 
detailed evidence is available. If the principles were completely false the 
subject of the book could hardly be said to exist, in that practically no 
uniformity could be imposed on the data. In fact they are broadly true, 
although only to a rough approximation. 

He then develops a number of further principles, and devotes the rest 
of the book to showing how they illuminate the main features of turbulent 
flow in wakes, jets, mixing zones, pipes, channels, boundary layers, and 
between rotating cylinders, and to showing how well the various principles 
are borne out by experimental observation. 

One principle that could be expressed better than in the book itself is 
that when turbulence is subject to a constant mean rate of strain the 
anisotropy of its energy tensor tends to a limiting form. If 


S,, = 4(U,/ex; + 3U;/ax,) 
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represents the mean rate of strain and u,u; the turbulent energy tensor, this 
anisotropy is represented by the behaviour of u,u;/q2 (where q?2 = u,u,). 
Chis must depend on S,; in a manner invariant under change of axes, and 
hence it can only be of the form 
Bh + Bg + C84 Ba, (1) 
where A, B, C are functions of the invariants S;,S,, and A (the determinant 
of S;). A form which fits the experimental observations cited by 
Dr ‘Townsend is 
Ui; U, 0-28 
— = 40,5 — so (2) 
» ot Y v 2 be —_— 
(7 S37) 
‘This seems more satisfactory than Dr ‘Townsend’s form, which is not 
invariant. "lhe reason why C in (1) is put equal to zero is that in plane 


strain (in the x, = 0 plane) the value of Dr ‘Townsend’s K, = 3uj/q?—1 
is observed to vanish. ‘The best-known case of equation (2) is the tendency 
of the ‘shear coefficient’ (—uz/u?) to be about 0-4 in simple shearing flow 


with dU/dy > 0. 


The physical explanation of the principle just described is based on the 
one hand on the known tendency for strain in one direction to increase 
vorticity in that direction and hence to favour velocity components 1n the 
direction at right angles (at the expense of components in that direction). 
To balance this tendency we no longer have the once-credited * general 
tendency to isotropy’ of the energy-containing eddies, which this book 
(together with a number of other contemporary investigations) finally 
explodes. If, in fact, the mean rate-of-strain is removed, the anisotropy 1s 
found to remain practically unchanged (until some new mean rate-of-strain 
is applied). In place of this, we have the suggestion that eddies too much 
stretched in the direction of mean strain decay more rapidly than less elon- 
gated ones. 

A second important principle inferred by Dr ‘lownsend from the 
data takes its simplest form in free turbulence. ‘This consists of a region 
of ‘uniformly turbulent’ fluid (containing all the vorticity and practically 
all the turbulent energy, uniformly distributed, but with anisotropy as 
governed by the former principle) with a bounding surface of extremely 
irregular shape, different at different instants. ‘The uniformity of the 
turbulent fluid is explained as due to diffusion. ‘The spread of the turbulent 
region is explained as due to contortion of the bounding surface by ‘large 
eddies’ and, on a smaller scale, by all the eddies, which taken together 
produce an extremely large boundary area over which viscous forces can 
diffuse vorticity. ‘That this process operates very effectively at high 
Reynolds numbers is clear from the fact that the observed rate of growth 
does not decrease as the Reynolds number increases (that is, as the time 
scale for operation of viscous forces increases relative to the basic kinematic 
time scales). 

The ‘large’ eddies (larger even that those which contain the bulk of the 
turbulent energy control the ultimate speed of spread of the turbulence. 
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‘They derive energy from the primary instability of the mean flow and pass 
energy down to the energy-containing eddies. Here Dr ‘Townsend assumes 
(with some experimental support) that one of the simplest kinds of eddy 
capable of extracting energy from the mean flow predominates in this range 
of sizes. ‘Vhis ends up as a round eddy (with axis in the direction of flow) 
which just fills up the region of substantial mean rate-of-strain of any one 
sign. Before the eddy is sheared into this form it is bent back at an angle 
and during this shearing process it gains energy (as the area of its closed 
streamlines decreases). 

Dr ‘Townsend derives the energy of the main turbulence by applying 
a condition that it extracts energy from a ‘large eddy’ of this form at the 
same rate at which the large eddy extracts energy from the mean flow. 
Thus, he is applying arguments similar to the ‘equilibrium’ arguments of 
Kolmogorott at the small-wave-number end of the scale. ‘The rate of 
dissipation of energy of the large eddy by turbulence is computed by the 
use of an ‘eddy viscosity’, and it is the magnitude of this eddy viscosity 
(determined empirically in the older theories) which emerges from the 
Townsend theory as a predicted quantity in satisfactory agreement with 
experiment. 

‘The mere existence of an effective eddy viscosity is not, however, taken 
as an a priori assumption. Dr ‘Townsend shows how it follows from his 
theory of the effect of mean rate-of-strain on anisotropy. One of his 
‘large eddies’ effectively re-orients the local rate-of-strain environment of 
the turbulence, and this alters the local anisotropy of the energy-containing 
eddies and so produces an addition to the Reynolds stresses. ‘The derivation 
is simpler from equation (2) above that from ‘lownsend’s form, and so 
may as well be given. A change 4S,; in the mean rate of strain produces 
a change 
5S, a S88. , 
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in the Reynolds stresses. In the case of a plane shearing motion only the 
term S,, = Sj, = }dU/dy is present in S,,, and then (3) becomes 
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which, for an eddy without an especially large 5S,. component of rate-of- 
strain, represents closely a viscous resistance by an effective viscosity equal to 
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which is the ordinary ‘eddy viscosity’ of the old theories. 
Throughout the book Dr ‘Townsend walks a very thin tight-rope in his 
attitude to eddy viscosity. He expects the reader to believe that it is 
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‘accidental’ (p.161) that the assumption of an eddy viscosity in the 
ordinary sense, which at any cross-section of the turbulent layer is constant 
within the ‘turbulent fluid’ and zero outside it, gives dead-on agreement with 
the observed mean velocity for all free turbulent flows. Earlier (p. 44), he 
admits the success with which Heisenberg and his followers applied the 
eddy viscosity idea to predict the whole spectrum of the energy-containing 
eddies in isotropic turbulence, but presumably he regards this as accidental, 
too! The reviewer would agree that Dr ‘Townsend’s data show that there 
is no simple single reason why eddy viscosity seems to work, but rather a 
number of conspiring circumstances—for example, there is Dr ‘Townsend’s 
point that where dl//dy is nearly zero, the eddies consist partly of those 
which have been positively sheared and partly of those which have been 
negatively sheared, and hence #7 also will be nearly zero. However, the 
rather sweeping and undeveloped arguments on which the older ideas are 
rejected in this book (for example, on p.94) do not greatly reinforce the 
simple objection (p. 128) that they contain undetermined constants. 

Chapters 4 to 8 are devoted to the development of the principles described 
above and their application to flow in wakes and jets. In these applications, 
it is clear that something has been sacrificed by the decision to study only 
flows which have reached the self-preserving state (principle (i1)), because 
the time taken by turbulence to change from one type to another is in fact 
long. In particular, the strong ‘mixing-region’ turbulence takes a long 
time to settle down into the relatively weaker wake turbulence or jet turbu- 
lence. On the other hand, the success of this book in throwing light on 
turbulence completely vindicates Dr ‘lownsend’s decision to confine his 
observations at this stage to the regions of self-preserving flow, far downstream, 
where this settling-down has taken place, and where there was more hope 
of uncovering the mechanics of the motion. 

In chapter 9, on flow in pipes and channels, the complications due to the 
presence of a wall are first introduced. Correlation and other measurements 
in the region of early constant stress near the wall are shown to be consistent 
with the idea that cylindrical eddies like the ‘large eddies’ of free turbulent 
flow sit on the wall and extract energy from the mean flow at the same rate 
as they pass energy down to not-so-large eddies. ‘This leads to the right 
ratio of eddy viscosity to distance from the wall. 

In chapters 10 and 11 the ideas appropriate to free turbulent flow and 
to flow near a wall are combined in a discussion of the turbulent boundary 
layer. ‘The theory of equilibrium of the large eddies is still remarkably 
effective and the evidence from correlation measurements for their presence 
and general character is good. In chapter 11 the ideas of the book (includ- 
ing, it must be noted, an eddy viscosity!) are used effectively on the interesting 
self-preserving flows with adverse pressure gradient recently discovered 
by Dr Clauser. In both chapters a discussion of three-dimensional 
effects including secondary flows is given. Finally, in chapter 12, the flows 
between rotating cylinders are fitted to some extent into the general frame- 
work of the book. 
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A large measure of sound understanding, with some notable gaps, was 
already present in the world literature on turbulent shear flow, so far as the 
behaviour of the mean flow and the small eddies was concerned. Dr ‘Town- 
send’s book, while not omitting this material, will always be remembered 
principally for having first made a reasonably convincing bridge between 
these extreme features of the flow and given us a pattern for the intermediate 
parts of the spectrum and a way of thinking quantitatively about the larger 
eddies that is sufficiently consistent with both the theoretical and the 
experimental evidence. 

M. J. LIGHTHILL 








